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I. INTRODUCTION 

In 1998 the Reviews of Modern Physics pubhshed a 
summary and critical analysis of the nuclear reaction 
cross sections important to solar burning. That effort, 
Adelberger et al.\ ( 1998 ) and denoted here as Solar Fu- 
sion I, began with a meeting hosted by the Institute for 
Nuclear Theory, University of Washington, 17-20 Febru- 
ary 1997. A group of international experts in the nuclear 
physics and astrophysics of hydrogen-burning stars met 
to begin critical discussions of the existing data on rele- 
vant nuclear reactions, with the aim of determining "best 
values" and uncertainties for the contributing low-energy 
S-factors. The group also considered opportunities for 
further improvements in both measurements and theory. 

Such data and related nuclear theory have been cru- 
cial to the standard solar model (SSM) and the neutrino 
fluxes it predicts. Indeed, measurements of nuclear re- 
actions gave the field its start. In 1958 jHolmgren and| 
Johnston|(|l958l|l959|) showed that the ratefor^He+^He 



r- 'Be -I-7 was ~ 1000 times larger than expected, and 
thus that the pp chain for *He synthesis would have addi- 
tional terminations beyond '^He-l-'^IIe ^He -I- 2p. This 
result led Davis to recognize that his chlorine detector 
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might be able to see the higher energy neutrinos from 
these other terminations, and spurred BahcaU and oth- 
ers to develop a quantitative model of the Sun capable 



the dual discoveries of Super-Kamiokande ( Fukuda et al. 



2001 ) and SNO (Ahmad et al. 2001 1 - two distinct neu- 



of predicting those fluxes (Bahcall and Davis Jr. 1982 1. 



At the time of the 1997 meeting, three decades of ef- 
fort in solar neutrino physics had produced four measure- 
ments that were at variance with the SSM and the stan- 
dard model of electroweak interactions. The measure- 
ments came from the pioneering work of Ray Davis, Jr. 
(Davis Jr. 1994 Davis Jr. et al. 1968 1; the observation 



of neutrinos in the Kamiokande water Cerenkov de- 



tector (Fukuda et al. 1996 


); and the GALLEX ( 


Kirsten 


et al. 


2003) and SAGE ( 


Gavrin et al. 


2003 


) radiochemi- 



The resulting pattern of fluxes that emerged from these 
experiments was difficult to reconcile with any plausible 
variation in the SSM, requiring a much sharper reduction 
in the "^Be neutrino flux than in the ®B flux, despite the 
greater sensitivity of the latter to changes in the solar 
core temperature. 

For this reason it was argued in Solar Fusion I that 
the measurements provided evidence for new physics be- 
yond the standard model. New solar neutrino exper- 
iments that promised much more precise data - the 
50-kiloton successor to Kamiokande, Super-Kamiokande, 
and the heavy-water-based Sudbury Neutrino Observa- 
tory (SNO), with sensitivity to both electron and heavy- 
flavor neutrinos - were then underway. The authors of 
Solar Fusion I, recognizing that the impact of these new 
experiments would depend in part on the quality of the 
nuclear microphysics input to the SSM, thus undertook 
an extended study of the key reaction rates for the pp 
chain and CNO bi-cycle. The effort appears to have been 
of some value to the community, as Solar Fusion I has 
become one of the most heavily cited papers in nuclear 
astrophysics. 



A. Solar Fusion II: the 2009/10 effort 

Ten years after publication of Solar Fusion I a proposal 
was made to the INT to revisit this process, in order to 
produce a new evaluation that would reflect the consid- 
erable progress made in the past decade, as well as new 
motivations for further constraining the SSM. Examples 
of advances in the nuclear physics include the LUNA II 



program at Gran Sasso ( jCostantini et al.\ 20091, which 
has provided remarkable low-energy measurements of key 
reactions such as "^116(0!, 7)^Be and ^'*N(p,7)^^0; several 
high-precision measurements addressing the key pp-chain 
uncertainty identified in Solar Fusion I, ''Be(p,7)^B; the 
application of new theoretical techniques to the p-t-p and 
hep neutrino reactions; and the resolution of several unre- 
solved questions about screening corrections in plasmas. 

The context for these measurements has also changed. 
In 1997 the field's central concern was, in some sense, a 
qualitative one, the origin of the solar neutrino problem. 
This question was answered in spectacular fashion by 



trino oscillations responsible for the missing atmospheric 
and solar neutrinos, largely determining the pattern of 
the light neutrino masses. But issues remain, and most 
of these require precision. There is intense interest in ex- 
tending direct measurements to the low-energy portion 
of the solar neutrino spectrum (< 2 MeV), where exper- 
iments with good energy resolution can determine the 
separate contributions of pep, CNO, ^Be, and pp neutri- 
nos. There is the potential to further constrain the solar 
neutrino mixing angle 012'. the solar luminosity deter- 
mines the pp flux to high accuracy, and the low-energy 
spectrum lies in the vacuum region of the MSW trian- 
gle, in contrast to the high-energy *B neutrinos, where 
matter effects are significant. Thus precise low-energy 
measurements have considerable "leverage" to test 9i2 
and the consistency of the conclusions we have drawn 
from SNO, Super-Kamiokande, and the KamLAND re- 
actor neutrino experiment. Borexino, now entering its 
calibration phase, is the first effort in this program of 
high-precision spectroscopy of low-energy solar neutrinos. 

But the resolution of the solar neutrino problem has 
also returned the field to its roots: Davis built the chlo- 
rine detector to probe the interior of the Sun and thereby 
test directly the theory of stellar evolution and nuclear 
energy generation (Bahcall and Davis Jr. 1982 1. Davis 



was diverted from that goal by the missing solar neutri- 
nos. But as the weak interaction effects responsible for 
that anomaly are now reasonably well understood, solar 
neutrinos again have become a quantitative tool for as- 
tronomy. Indeed, the program carried out by SNO and 
Super-Kamiokande has already yielded one remarkable 
constraint on the Sun, a direct determination of the core 
temperature to high precision, through measurement of 
the ^B neutrino flux (0(^B) cx T^'^). The 8.6% precision 



of the SNO NCD-phase results ( [Aharmim et al.\ |2008| ), 
(/)(8B) = {5. 54^0 fit oil) X lOVcmVs, implies a sensitiv- 
ity to core temperature of ^ 0.5%. 

New questions have arisen about the Sun that neu- 
trinos could potentially address, provided the associated 
laboratory astrophysics has been done. One important 
success of the SSM in the 1990s was in predicting the 
local sound speed c(r). Comparisons between c(r) de- 
duced from helioseismology and the predictions of the 
SSM yielded agreement at ~ 0.2% throughout much of 



the Sun. Bahcall and others argued (Bahcall et al. 2001 ) 
that helioseismology is a more severe and detailed test of 
the SSM than neutrino production, so that SSM success 
in reproducing c(r) made a particle-physics resolution of 
the solar neutrino problem more likely. 

The sound speed is a function of the Sun's interior 
pressure and density profiles, which in turn reflect ther- 
mal transport properties that depend on the Sun's metal 
content, through the opacity. Thus the comparison be- 
tween helioseismology and the SSM tests a key assump- 
tion of the SSM, that the metals are distributed uni- 
formly throughout the Sun, apart from small corrections 
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due to diffusion. This assumption allows one to equate 
SSM interior metal abundances to convective-zone abun- 
dances deduced from analyses of photospheric absorp- 
tion lines. Such analyses had been based on ID models 
of the photosphere. Recently ab initio 3D analyses have 
been developed, yielding significant improvements in pre- 
dicted line shapes and in the consistency of metal abun- 
dance determinations from various atomic and molecular 
lines. However, this work also reduced metallicity es- 
timates from Z ~ 0.0169 to ~ 0.0122 (lAsplund et al. 



2005), destroying the once excellent agreement between 



hehoseismology and the SSM. 

It has been suggested that this difficulty may re- 
flect, contrary to the SSM, differences in solar core and 
convective-zone metallicities that could have arisen from 
the late-stage evolution of the solar disk: as a great deal 
of metal was scoured out of the disk by the formation of 
the giant planets, the last few percent of gas deposited 
onto the Sun could have been depleted of metals ( ^Hax-, 
ton and Serenelli 2008). Indeed, recent studies of "solar 



twins" show abundance trends that correlate with the ex- 



istence of planets ( Israelian et al. 2009 Ramirez et al. 



2009 ) . Haxton and Serenelli ( 2008 ) argued that a direct 



measurement of solar core metallicity could be made by 
observing CNO solar neutrinos. 

In both of the above examples - using neutrinos to 
determine the solar core temperature and metallicity - 
nuclear physics uncertainties remain one of the limiting 
factors in the analyses. 

The proposal to revisit in 2009 the deliberations of 
1997 thus had several motivations: 

• providing a set of standard S-factors and uncertain- 
ties that reflect the progress made in laboratory 
and theoretical nuclear astrophysics over the last 
decade; 

• enabling more precise analyses of solar neutrino ex- 
periments designed to constrain neutrino oscilla- 
tions and other new physics, e.g., future pp and 
pep neutrino experiments that exploit these well 
understood fluxes; and 

• enabling analyses in which solar neutrinos are used 
as a probe of the solar core. 

The 2009 INT workshop^ was modeled after that of 1997, 
with invitations extended to and accepted by representa- 



tives from most of the experimental groups active in the 
nuclear physics of hydrogen burning stars. There was also 
active involvement of theorists, reflecting the progress 
that has been made in ab initio calculations. The work- 
shop participants are the authors of this manuscript. As 
in 1997, early organizing included the selection of working 
group leaders who identified key papers, which were then 
entered in a database for review, prior to the start of the 
workshop. These materials were then summarized and 
discussed during the workshop, as the various working 
groups considered the state of the data and outlined any 
additional work that would be needed for this review. 
The process of critically analyzing both new and older 
data and working toward a consensus on best-value cross 
sections and uncertainties continued throughout 2009. A 
few new topics not considered in 1997 but now recog- 
nized to be quite important, such as the shape of the 
B neutrino spectrum, were addressed. (The neu- 
trino spectrum is one of the inputs to SNO and Super- 
Kamiokande analyses.) The workshop included working 
groups on indirect techniques for constraining cross sec- 
tions, to summarize the progress that has been made in 
validating such approaches, and on new facilities and in- 
strumentation, in view of the facility investments that 
are being considered in laboratory nuclear astrophysics 
(above and below ground). 



B. Contents of this review 

The review begins in Section II with a description of 
hydrogen burning by the pp chain and CNO bi-cycle, and 
the neutrino byproducts of these reaction chains. The 
role of S-factors and the associated questions of screen- 
ing and of extrapolating data to the solar Gamow peak 
are discussed. We provide a fairly complete overview of 
progress in theory, which in some cases provides our only 
estimate of S-factors, and in other cases determines the 
forms of the functions that are needed for data extrapo- 
lations. 

Discussions of individual reactions are organized by 
chapter: Sees. III-IX discuss the pp chain reactions p-|-p 
A+e++Ve] d-hp -J> ^He-h7; ^He-h^He "^He-l-p-Hp; 
^He-h^He ^Be-h7; ^He-Kp ^ ^Yle+e++Ve\ '^Be, pp, 
and CNO nuclei electron capture; and ''Be-l-p — *B-|-7. 
Sec. X discusses the spectrum of *B neutrinos produced 
in the (3 decay to a broad resonance in ^Be. Sec. XI dis- 



cusses N-|-p 
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0+7 and other reactions contribut- 



The workshop was proposed in a letter to the Institute for 
Nuclear Theory's National Advisory Committee (NAC) and ap- 
proved by the NAC and INT Director at the time of the NAC's 
August 2008 annual meeting. Wick Haxton (lead), Eric Adel- 
berger, Heide Costantini, Peter Parker, R. G. Hamish Robertson, 
Kurt Snover, Frank Strieder, and Michael Wiescher formed the 
organizing committee and served as co-editors of this paper. Ad- 
ditional community members joined this group to act as working 
group heads: Jiunn-Wei Chen, Barry Davids, Stuart Freedman, 
Alejandro Garcia, Uwe Greife, Michael Hass, Gianluca Imbri- 
ani, Kuniharu Kubodera, Daniela Leitner, Laura Marcucci, 



Filomena Nunes, Tae-Sun Park, Paolo Prati, Hanns-Peter 
Trautvetter, and Stefan Typel. The working group heads were 
responsible for organizing discussions, creating section drafts, 
and responding to subsequent criticisms of the drafts. Organiz- 
ing committee members, in their capacity as co-editors, were 
responsible for creating from the drafts a coherent document, 
and for addressing any issues unresolved by the working groups. 
Workshop presentations are archived on the INT's web site, 
^http:/ /www. int. Washington. edu/PROGRAMS/solar_fusion.html| 
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ing to the CNO cycles. Sec. XII describes the progress 
that has been made in developing and validating indirect 
methods, while Sec. XIII describes future facilities and 
instrumentation that could further advance the field. 

The conclusions of this review, in some cases, required 
the working groups to make some judgments. There are 
discrepant data sets, and there are cases where data ex- 
trapolations have some dependence on models. We have 
tried to treat such questions as consistently as possible, 
aware that excessively optimistic treatments of uncer- 
tainties could be misleading, while excessively conser- 
vative treatments would degrade the value of the best 
experiments done in the field. In most cases our working 
groups were able to reach consensus. In cases where sig- 
nificant differences remained among the experts, we have 
tried to identify the source of the disagreement, so that 
"consumers" will be aware that full consensus may have 
to await future measurements. 

Table U summarizes the conclusions of this review. 



II. NUCLEAR REACTIONS IN HYDROGEN-BURNING 
STARS 

Observations of stars reveal a wide variety of stellar 
conditions, with luminosities relative to solar spanning 
a range L ~ 10"'* to 10^ Lq and surface temperatures 
Ts -2000-50000 K. The simplest relation one could pro- 
pose between luminosity L and Tg is 

L = 47ri?VsB T* ^ L/Lq = [R/Rq? {Ts/Tq)\ (1) 

where ctsb is the Stefan-Boltzmann constant, and Lq, 
Tq, and i?0 are the solar values. This relation suggests 
that stars of a similar structure might lie along a one- 
parameter path (in this simplified example, defined by a 
fimction of the blackbody radii, (i?/i?o)^) in the luminos- 
ity (or magnitude) vs. temperature (or color) plane. In 
fact, there is a dominant path in the Hertzsprung-Russell 
color-magnitude diagram along which roughly 80% of the 
stars reside. This is the main sequence, those stars sup- 
porting themselves by hydrogen burning through the pp 
chain, 

4p->%e + 2e+-f-2i/c, (2) 

or CNO cycles. The laboratory nuclear astrophysics of 
hydrostatic hydrogen burning is the focus of this review. 

As one such star, the Sun is an important test of our 
theory of main sequence stellar evolution: its properties - 
age, mass, surface composition, luminosity, and helioseis- 
mology - are by far the most accurately known among the 
stars. The SSM traces the evolution of the Sun over the 
past 4.6 Gyr of main sequence burning, thereby predict- 
ing the present-day temperature and composition pro- 
files, the relative strengths of competing nuclear reac- 
tion chains, and the neutrino fiuxes resulting from those 
chains. The SSM makes four basic assumptions: 



• The Sun evolves in hydrostatic equilibrium, main- 
taining a local balance between the gravitational 
force and the pressure gradient. Knowledge of the 
equation of state as a function of temperature, den- 
sity, and composition allows one to implement this 
condition in the SSM. 

• Energy is transported by radiation and convection. 
The solar envelope, about 2.6% of the Sun by mass, 
is convective. Radiative transport dominates in 
the interior, r < O.72i?0, and thus in the core 
region where thermonuclear reactions take place. 
The opacity is sensitive to composition. 

• The Sun generates energy through hydrogen burn- 
ing, Eq. ([2]). Figure [1] shows the competition be- 
tween the pp chain and CNO cycles as a function 
of temperature: the relatively cool temperatures of 
the solar core favor the pp chain, which in the SSM 
produces — 99% of the Sun's energy. The reactions 
contributing to the pp chain and CNO bi-cycle are 
shown in Fig. [2] The SSM requires as input rates 
for each of the contributing reactions, which are 
customarily provided as S-factors, defined below. 
Typically cross sections are measured at somewhat 
higher energies, where rates are larger, then extrap- 
olated to the solar energies of interest. Corrections 
also must be made for the differences in the screen- 
ing environments of terrestrial targets and the solar 
plasma. 

• The model is constrained to produce today's solar 
radius, mass, and luminosity. The primordial Sun's 
metal abundances are generally determined from a 
combination of photospheric and meteoritic abun- 
dances, while the initial ^He/H ratio is adjusted to 
reproduce, after 4.6 Gyr of evolution, the modern 
Sun's luminosity. 

The SSM predicts that, as the Sun evolves, the core He 
abundance increases, the opacity and core temperature 
rise, and the luminosity increases (by a total of — 44% 
over 4.6 Gyr). The details of this evolution depend on a 
variety of model input parameters and their uncertain- 
ties: the photon luminosity Lq, the mean radiative opac- 
ity, the solar age, the diffusion coefficients describing the 
gravitational settling of He and metals, the abundances 
of the key metals, and the rates of the nuclear reactions. 

If the various nuclear rates are precisely known, the 
competition between burning paths can be used as a sen- 
sitive diagnostic of the central temperature of the Sun. 
Neutrinos probe this competition, as the relative rates of 
the ppl, ppll, and ppIH cycles comprising the pp chain 
can be determined from the fiuxes of the pp/pep, ^Be, 
and neutrinos. This is one of the reasons that labora- 
tory astrophysics efforts to provide precise nuclear cross 
section data have been so closely connected with solar 
neutrino detection. 

Helioseismology provides a second way to probe the 
solar interior, and thus the physics of the radiative zone 
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TABLE I The Solar Fusion II recommended values for S(0), its derivatives, and related quantities, and for the resulting 
uncertainties on S{E) in the region of the solar Gamow peak - the most probable reaction energy - defined for a temperature 
of 1.55 X lO^K characteristic of the Sun's center. See the text for detailed discussions of the range of validity for each S{E). 
Also see Sec. |VIII| for recommended values of CNO electron capture rates. Sec. |XI.B| for other CNO S-factors, and Sec. |X]for 
the *B neutrino spectral shape. Quoted uncertainties are la. 



Reaction 


Section 


S(0) 


S'(0) 


S"(0) 


Gamow peak 






(keV-b) 


(b) 


(b/keV) 


uncertainty (%) 



p(p,e+i^e)d 


III 


d(p,7)'He 


IV 


^He(^He,2p)*He 


Pi 


^He(^He,7)''Be 




^He(p,e+!/e)*He 


VII 


■^Be(e-,//e)^Li 


V^III 


p(pe",!^e)d 


vm 


^Be(p,7)«B 


IX 


"N(p,7)i50 : 


■CI.A 



(4.01 ± 0.04) X 10" 
(2.14 



+0.17N 



xlO" 



(5.21 ± 0.27) X 10^ 

0.56 ± 0.03 
(8.6 ± 2.6) X 10"^° 
See Eq. Jioi 



(4.49 ± 0.05) X 10"^ 
(5.56l°i«)xl0-« 
-4.9 ± 3.2 
(-3.6 ± 0.2)xl0~'' 



(9.3+^;^) X 10-^ 
(2.2 ± 1.7) X 10"^ 
(0.151 ± 0.008)xl0"^ " 



See Eq. ||46l 
(2.08 ± 0.16)xl0 
1.66 ± 0.12 



(-3.1 ± Q.3)xl0-^ 
(-3.3 ± 0.2)xl0"^ 



(2.3 ± 0.8) xlO"'^ 
(4.4 ± 0.3)xl0"^ " 



± 0.7 
± 7.1 " 
=t 4.3 " 
± 5.1 
± 30 
± 2.0 
± 1.0 
± 7.5 
± 7.2 



"Error from phenomcnological quadratic fit. See text. 
''S'(0)/S(0) taken from theory; error is that due to S(0). See text. 
'^S"(0)/S(0) taken from theory; error is that due to S (0). See text. 
''Estimated error in the pep/pp rate ratio. See Eq. |46| 
'^Error dominated by theory. 



I ' I ' I ' I ' I ' I ' I ' I 



-4 ^ 



CN cycle 




pp-chain 



T (10 K) 

FIG. 1 The stellar energy production as a function of temper- 
ature for the pp chain and CN cycle, showing the dominance 
of the former at solar temperatures. Solar metallicity has 
been assumed. The dot denotes conditions in the solar core: 
the Sun is powered dominantly by the pp chain. 



that the SSM was designed to describe. The sound speed 
profile c(r) has been determined rather precisely over the 
outer 90% of the Sun and, as previously discussed, is now 
in conflict with the SSM, when recent abundance deter- 
minations from 3D photospheric absorption line analyses 
are used. 



A. Rates and S-factors 

The SSM requires a quantitative description of relevant 
nuclear reactions. Both careful laboratory measurements 
constraining rates at near-solar energies and a supporting 
theory of sub-barrier fusion reactions are needed. 

At the temperatures and densities in the solar inte- 
rior (e.g., Tc ~ 15.5 X 10^ K and pc ~ 153 g/cm^ at 
the Sun's center), interacting nuclei reach a Maxwellian 
equilibrium distribution in a time that is infinitesimal 
compared to nuclear reaction time scales. Therefore, the 
reaction rate between two nuclei can be written (Bur- 



bidge et al. 1957 Clayton 19681 



^12 



ni n2 



1 



'12 



(3) 



Here the Kronecker delta prevents double counting in 
the case of identical particles, ni and n2 are the number 
densities of nuclei of type 1 and type 2 (with atomic 
numbers Zi and Z2, and mass numbers Ai and A2), and 
{crv)-^2 denotes the product of the reaction cross section 
a and the relative velocity v of the interacting nuclei, 
averaged over the collisions in the stellar gas. 



{f^V)i2 



(7{v) V ^{v) dv. 



(4) 



Under solar conditions nuclear velocities are very well 
approximated by a Maxwell-Boltzmann distribution. It 
follows that the relative velocity distribution is also a 
Maxwell-Boltzmann, governed by the reduced mass ^ of 
the colliding nuclei. 



^^^] Any' 
2kT 1 



dv. 



(5) 
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FIG. 2 The left frame shows the three principal cycles comprising the pp chain (ppl, ppll, and ppIII), with branching percentages 
indicated, each of which is "tagged" by a distinctive neutrino. Also shown is the minor branch ^He+p — >■ '*He+e^+!/e, which 
burns only ~ 10"*^ of "^He, but produces the most energetic neutrinos. The right frame shows the CNO bi-cycle. The CN cycle, 
marked I, produces about 1% of solar energy and significant fluxes of solar neutrinos. 



Therefore, 



yields (Bahcall 1989) 
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{av)i2 = 



7rfj,{kT) 



(6) 

where E is the relative kinetic energy and k is the Boltz- 
mann constant. In order to evaluate (crv) 12 the energy 
dependence of the reaction cross section must be deter- 
mined. 

Almost all of the nuclear reactions relevant to solar 
energy generation are nonresonant and charged-particle 
induced. For such reactions it is helpful to remove 
much of the rapid energy dependence associated with the 
Coulomb barrier, by evaluating the probability of s-wave 
scattering off a point charge. The nuclear physics (in- 
cluding effects of finite nuclear size, higher partial waves, 
antisymmetrization, and any atomic screening effects not 
otherwise explicitly treated) is then isolated in the S- 
f actor, defined by 



a{E) 



SiE) 
E 



exp [-27r77(£')] 



(7) 



with the Sommerfeld parameter "q^E) = Z1Z2 cx/v, where 
V = ^2E I \x is the relative velocity and a the fine struc- 
ture constant (?i = c = 1). Because the S-factor is slowly 
varying, one can extrapolate S(i?) more reliably from the 
range of energies spanned by data to the lower energies 
characterizing the Gamow peak. 

A substitution of Eq. ([t]) into Eq. ^ followed by a Tay- 
lor expansion of the argument of the exponentials then 



{(yv)i2 



2 A^Jo 
likT kT 



fo S,s exp [-3Eo/{kT)] 



1.301 X 10"^^ cm7s 



Z1Z2 
A 



1/3 



^" T^9~'^'exp[-3£;o/(fcT)], (8) 



where 



Eo ^ 
kT 

AEo 
kT 

and 

Seff = S(0) ( 1 



i7:Z,Z2a/V2)'/^ [ti/{kT)f\ 



Eq 
3kT'- 



5kT 
36Eo 



A 



A1A2 
A1 + A2' 



S\0)Eo 



35kT 
■MEo 



+ -S"{0)E^ 1 + 



89kT\ 
36^j 



Eq, the Gamow peak energy where the integrand of 
Eq. ^ takes on its maximum value, is the most prob- 
able energy of reacting nuclei. A.Eq corresponds to the 
full width of the integrand at 1/e of its maximum value, 
when approximated as a Gaussian. Equation Q includes 
a factor fo, discussed below, to correct for the effects of 
electronic screening on nuclear reactions occurring in the 
solar plasma. 
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Rates in an astrophysical plasma can be calculated 
given S{E) which by virtue of its slow energy dependence, 
in the case of non-resonant reactions, can be approxi- 
mated by its zero-energy value S(0) and possible correc- 
tions determined by its first and second derivatives, S'(0) 
and S"(0). It is these quantities that we need to deter- 
mine by fitting laboratory data, or in cases where such 
data cannot be obtained, through theory. For most of 
the reactions contributing to the pp chain and CNO bi- 
cycle, data have been obtained only for energies in regions 
above the Gamow peak, e.g., typically E > 100 keV, so 
that extrapolations to lower energies depend on the qual- 
ity of the fit to higher energy data. Ideally one desires a 
fitting function that is well motivated theoretically and 
tightly constrained by the existing, higher-energy data. 
The purpose of this review is to provide current best val- 
ues and uncertainties for S(0) and, if feasible, its deriva- 
tives. 

S-factor uncertainties, when folded into SSM calcula- 
tions, then limit the extent to which that model can pre- 
dict observables, such as the depth of the convective zone, 
the sound speed profile, and the neutrino fluxes. It has 
become customary in the SSM to parameterize the con- 
sequences of input uncertainties on observables through 
logarithmic partial derivatives, determined by calculat- 
ing the SSM response to variations in individual input 
parameters. SSM compilations of the logarithmic par- 
tial derivatives provide, for example, a way to assess the 
importance of each S-factor uncertainty on neutrino flux 
predictions. 

The partial derivatives a{i,j) for each neutrino fiux 0^ 
and SSM input parameter /3j are defined by 



01n[^,M(0)] 
91n[/3,//?,(0)] 



(9) 



where 0i(O) and f3j{0) denote the SSM best values. The 



a{i,j) for 19 SSM input parameters f3j are given by Peha- 
Garay and SerenelUl ([2008) in their 2008 SSM update. 
The I3j include parameters such as the Sun's age and 
luminosity, the abundances of important metals, and S- 
f actors. 

The partial derivatives define the power-law dependen- 
cies of neutrino fluxes with respect to the SSM best- value 
prediction 0^(0), 



N 

.(0)11 



ft- 



-| a{i,j) 



ft(0) 



N 

.(0)11 [1 



<5A 



where the product extends over N SSM input parame- 
ters, and where S/3j = Af3j//3j{0) is the fractional un- 
certainty of input parameter /3j with respect to its SSM 
best value. This expression separates the impact of SSM 
parameter variations on (jji into a solar piece - the in- 
flnitesimal SSM response described by a{i,j) - and a 
laboratory or theory piece - the estimated uncertainty 
6j3j of an input parameter (in our case, that of an S- 
factor). From SSM tabulations of the a{i,j), one can 



estimate the change in a SSM flux prediction (pi, when a 
given SSM parameter Pj is perturbed away from its SSM 
best value by an amount S/3j, without redoing the SSM 
calculation. For example, to assess the impact of an im- 
proved nuclear cross section measurement on (pi , one sets 
df3j to the estimated uncertainty of the corresponding S- 
factor, to obtain the corresponding variation in (pi. In 
this way one can identify nuclear physics improvements 
that will have the most impact on reducing flux uncer- 
tainties. Alternatively, the process can be inverted: a 
fiux measurement could in principle be used to constrain 
an uncertain input parameter. 



For example, Pena-Garay and Serenelli (2008) define 
the dependence of (p{'^B) on the S-factors under discus- 
sion here, 



SS, 



0(«B) oc (1 + 5Sn)-'-''(l + <5S33)-"-""(l + 

x(l + <5Si7)i-"(l + <5Se7)-i-°(l + <5Si 14)-°-"', (11) 

where Sn denotes the S-factor for p+p reaction, etc., and 
SSii = ASii/Sii(0) denotes its fractional uncertainty. 
This review gives the best current values for the needed 
SSs. 



B. Screening of stellar and laboratory reactions 

One must take into account differences in the atomic 
environments to correctly relate screened laboratory and 
solar cross sections, (jI^^{E) and a-^°^^''{E), to each other 
or to the underlying bare cross section a]^{E). Screening 
enhances solar cross sections by reducing the Coulomb 
barrier that reacting ions must overcome. As light nuclei 
in the solar core are almost completely ionized, the solar 
electron screening correction /o. 



foiE) 



-solar 



(E) 



(12) 



can be treated in a weak-screening approximation 



(Salpeter 1954 1. The impact of the modified potential. 



aZiZ2 
V [r) = exp 



RdJ' 



(13) 



on reactions depends on the ratio of the Coulomb poten- 
tial at the Debye radius Rd to the temperature. 



(10) fo 



exp 1 ) = exp (0.I88 Z, Z2 C Po^' T, ^'^ 



(14) 

1 /2 

where C^Rd = [fcT/(47rap)]''", p is the num- 
ber density of nucleons, po is a dimen- 
sionless density measured in g/cm'^, ^ = 

-,1/2 

E X. (^fMO+(/o//o)E ^.(^.M.) , is 

the mass fraction of nuclei of type i, and Tg is the 
dimcnsionless temperature in units of 10^ K. The 
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factor /0//0 ~ 0.92 corrects for the effects of electron 



degeneracy in the solar core ( Salpeter 1954 ) 



The weak-screening approximation requires the aver- 
age interaction energy between particles to be smaller 
than the average particle kinetic energy ( [Baimbetov 



et al. 1995 Kobzev et al.\ 1995 ). This places a con 



straint on the argument of Eg. (|14[), Zi Z^aj (RokT) <C 



1, that is satisfied in the solar core if Z1Z2 < 10 (Gruzi- 



nov 



PP c- 



1998), a condition met by the low-Z reactions ot the 
lain and CNO bi-cycle. However corrections to the 
Salpeter formula are expected at some level. Nonadia- 
batic effects have been suggested as one source, e.g., when 
a high Gamow energy guarantees reacting nuclei having 
velocities significantly higher than the typical ion veloc- 
ity, so that the response of slower plasma ions might be 
suppressed. At the time of Solar Fusion I such dynamic 
corrections were a source of controversy. Dynamic correc- 
tions were first discussed by Mitler ( 19771 and later stud- 
ied by Carraro et al. (1988). Subsequent work showed 



that Salpeter's formula would be valid independent of the 
Gamow energy due to the nearly precise thermodynamic 



equilibrium of the solar plasma ( 


Brown and Sawyer 1997 


Gruzinov 


1998 


Gruzinov and Bahcall , 1998). The argu- 



ments, summarized in Solar Fusion I, were significantly 
extended in 2002 by Bahcall et al. (2002), who pointed 



out a number of contradictions in investigations claim- 
ing larger corrections, and showed that a field theoretic 
approach led to the expectation of only small (~ 4%) 
corrections to the standard formula, for solar conditions. 



However controversies have not entirely died out (Mao 



et al. 2009). 



The Salpeter correction relates the solar and bare cross 
sections, ag°^'^''{E) and ah{E). As the reactions studied 
in the laboratory generally involve target nuclei bound in 
neutral atoms and molecules, not bare ions, a second step 
is needed to extract <j\,{E) from laboratory data. As in 
the Sun, electrons in the laboratory target tend to reduce 
the barrier, so that the screened cross section a'^^^{E) 
will exceed that for bare ions a^iE). The enhancement 



is given by (Assenbaum et al. 1987) 



fuh[E) 



,lab 



[E) 



^hiE) 



exp 



7rv{E)U, 
E 



> 1 for Ue < E, 
(15) 

where Ue is an electron-screening potential energy. This 
energy can be estimated from the difference in atomic 
binding energies between the compound atom and the 
projectile plus target atoms of the entrance channel. Be- 
cause the correction depends on the ratio Ue/E, one ex- 
pects screening corrections to be most important for very 
low projectile energy. 

In contrast with the case of solar screening, a great deal 



baum et al. 


1987 


|Engstler et al. 


1988 


1992 


Greife et al. 


1995 


Prati et al. 


1994 


Rolfs 


2001 Rolfs and Somorjai 



1995) to test our understanding of electron screening in 



terrestrial targets. Studies of reactions involving light nu- 



clei (Engstler et al. 1988 Strieder et al. 2001) revealed 



an u ptur n in cross section at low energies, as pred icted by 
Eq. ( 15). For example, results for ^He(d,p)^He (jAliotta 
et al. 120011) could be represented by Eq. (ITF 



2001 1 could be represented by Eq. ( [T5| for a screen- 
ing potential Ue = 219±15 eV. While this potential is sig- 
nificantly larger than the one obtained from the adiabatic 
approximation, Uad = 119 cV, the analysis requires one 
to assume an energy dependence of the bare cross section 
ah{E). This adds a difficult-to-quantify theoretical un- 
certainty to the extracted potential. It may be possible 
to remove much of this uncertainty through an indirect 
measurement of CTh{E) by the Trojan Horse Method (Lat 



tuada et al.] [200T| [Spitaleri et al.[ [200T| [Strieder et al. 
20011 iTumino et aVi 120031). 



There exist various surrogate environments that have 
been exploited by experimentalists to test our under- 
standing of plasma screening effects. Screening in d(d,p)t 
has been studied for gaseous targets and for deuterated 
metals, insulators, and semiconductors fR aiola et al.\ 
2004). For a summary of the results see iHaxton et al. 



(2006): it is beheved that the quasi- free valence electrons 
in metals create a screening environment quite similar 
to that found in stellar plasmas. Experiments in met- 
als have confirmed important predictions of the Debye 
model, such as the temperature dependence Ue{T) cx 

The tendency of experimentally determined values of 
Ue to exceed theoretical estimates by a factor ^ 2 has 
been noted by [Assenbaum et al] ( [1987[ ); [Rolfs[ ( [2001[ ); 
Rolfs and Somorjai' (1995). Various possible explanations 



have been considered (jBalantekin et al. 


1997 


Fiorentini 


et aL[[2003 


[Flambaum and Zelevinsky||1999 


agino and 


Balantekin 


[2002 


Shoppa et al. 


1993 


). A possible solu- 



tio n of the laborato r y scre ening pr oblem was p r oposed 
in Langanke et al. (1996) and in Bang et al. (19961, 
that the stopping of ions in matter differs at low en- 
ergy from that obtained by extrapolating from stopping 



power tables at higher energies (Andersen and Ziegler 



1977). Smaller stopping powers were indeed verified ex- 



perimentally ( Golser and Semrad 1991 Rolfs 2001) and 
explained theoretically (Bertulani 2004 Bertulani and 



de Paula 2000 ) 



Screening corrections for laboratory reactions are im- 
portant in extracting S-factors in cases where data extend 
to very low energies. In this review two cases of interest 
are '^He-l-'^He — > p-l-p-f-^He, where the lowest data point 
is at i? = 16 keV, and ^^N(p,7)^^0, where measurements 
extend down to 70 keV. 



C. Fitting and extrapolating S-factors 

S((3) (and its derivatives S'(0) and S"(0)) needed in 
Eq. (Is]) could be taken from a polynomial fit to data. 
A quadratic form often provides an excellent representa- 
tion of the data up to a few hundred keV. However, as 
the procedure is purely empirical, it provides no theoret- 
ical justification for extrapolating beyond the last known 
data point. For example, a quadratic fit to the labora- 
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tory data for ^Be(p,7)^B would miss the upturn in the 
S-factor at low energy expected from theory, as this in- 
crease occurs beyond the range of existing data. For this 
reason, we restrict our use of empirical fitting functions 
to cases where the data sets encompass the full range of 
energies relevant to astrophysics. 



1. Theory constraints: model-based methods 

One class of important theoretical constraints makes 
use of the peripheral nature of non-resonant radiative 
capture reactions close to the threshold. If the reaction 
occurs at separations much larger than the sum of the nu- 
clear radii, one can derive the coefficients for the Taylor 
series for S{E) independent of models, as only the asymp- 
totic forms of the bound and scattering initial- and final- 
state wave functions are relevant. This idea has been 
exploited in several ways. 



and Baur ( 2005 1 for explicit expressions of the cross sec- 



Williams and Koonin ( 1981 ) used Bessel function ex- 



pansions of Coulomb wave functions and a hard-sphere 
approximation to derive an expansion of the low-energy 
logarithmic derivative. 



1 dS{E) 
S{E) dE 



a + bE. 



(16) 



This approach was further 



Mukhamedzhanov and Nunes 



developed by 
(|2002|, who consid- 



ered variables such as the remnant Coulomb barrier, 
the initial and final centrifugal barriers, and the bind- 
ing energy (but not the interactions of the colliding 
nuclei in the entrance channel). They found that the 
near-threshold behavior of S{E) could be sensitive to 
such parameters. Baye and collaborators, employing 
zero-energy solutions of the Schrodinger equation and 
their energy derivatives, showed that model-independent 
values for the coefficients in the Taylor expansion for 
S{E) around E — could be extracted from the asymp- 
totic normalization coefficient (ANC) of the bound 
state wave function and the scattering lengths of the 
scattering states, thus including effects from interactions 



in the continuum (Baye 
Brainisl 120001). 



2000 2004 2005 Baye and 



Despite the successful application of the Taylor se- 
ries expansion for S{E), it was noticed that the series 
has a restricted domain of convergence, determined by 
the binding energy Eb of the final state. This is a 
consequence of a pole in the relevant radial integral at 
E = -Eb ( |Bayel [2000{ [Jennings e<~ami998a|b] ). This 
limitation becomes particularly severe for weakly bound 
nuclei: for ^Be(p,7)^B, \Eb\ ~ 138 keV barely reaches 
the domain of experimental data. Thus the alternatives 
of a Laurent expansion of the S-factor in the photon en- 
ergy E-y = E + Eb, an expansion of {E + Eb)S{E), and 
the explicit treatment of the pole have been explored as 
alternatives in the analysis of experimental data ( Cyburt 



and Davids 2008 Cyburt et al. 2004). See also Typel 



tions without the convergence limitation. 

Model-based calculations of fusion cross sections also 
provide a template for fitting and extrapolating experi- 
mental data. Models can be constrained by the known 
properties of the system under study and can be applied 
over a wide range of energies. While they often pre- 
dict the energy dependence of S(i?) accurately, in many 
cases an overall renormalization is needed to give the cor- 
rect magnitude of the S-factor. The need for this scal- 
ing is qualitatively understood, as model calculations of 
interior wave functions are generally done in restricted 
spaces, and thus lack high-momentum (and certain low- 
momentum) components of the true wave function, with 
consequences for the normalization. (The goal of pre- 
dicting both the shape and normalization of S-factors 
is motivating the development of quasi-exact ah initio 
methods, as discussed below.) 

Modeling approaches involve various levels of complex- 
ity. The simplest microscopic reaction theories are the 
potential models, in which the internal structure of the 
colliding nuclei is ignored. The dynamics of the process 
is reduced to a single coordinate, the distance vector be- 
tween the two nuclei. The potential-model Hamiltonian 
is typically a phenomenological one, e.g., a Woods-Saxon 
potential, with parameters that can be determined by 
fitting data, such as the elastic cross section. 

More realism is provided by cluster models like the 
resonating group method (RGM) or the generator- 
coordinate method (GCM), which take into account the 
many-body substructure of the reacting nuclei. These 
models employ fully antisymmetrized many-body wave 
functions of the compound system, though constructed 
in a restricted model space. The full wave function is 
described as a superposition of many-body cluster wave 
functions of fixed internal structure moving against each 
other. The interaction is described by phenomenologi- 
cal nucleon-nucleon potentials with parameters that are 
adjusted for each reaction under consideration. 

Another description of fusion reaction cross sections 
comes from the R-matrix. Space is divided into two re- 
gions, the interior where nuclear forces are important, 
and the exterior where the interaction between the nu- 
clei is assumed to be only Coulombic. The full scatter- 
ing wave function connecting different channels i is ex- 
panded in partial waves with total angular momentum J. 
The Schrodinger equation for the interior Hamiltonian is 
solved, with boundary conditions at the channel radii 
encoding the correct asymptotic behavior. The solutions 
of the Schrodinger equation determine the level energies 
E\ and reduced widths 7Ai that appear in the expression 
for the R-matrix 



Rij{E) 



N 



^ E - Ex' 

A=l ^ 



(17) 



for each J, in the standard approach of |Lane and Thomas] 
( 1958 ) . Simple expressions relate the reaction cross sec- 
tions at energy E to the R-matrix. The cross section 
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should be insensitive to the choice of the channel radii. 
In most applications the R-matrix is viewed as a param- 
eterization of measured reaction cross sections in terms 
of fitted level energies and reduced widths. A connection 
to an underlying reaction model is not required. The 
R-matrix allows one to properly account for penetrabil- 
ity effects, and to adjust the complexity of the fitting in 
response to various practical considerations, such as the 
energy range of interest. 

R-matrix resonance parameters (level energies and re- 
duced widths) are not directly comparable to the exper- 
imental quantities due to level shifts associated with the 
chosen boundary conditions. Generalizing earlier ideas of 
Barker (1971) and Angulo and Descouvemont (2000), an 



combine this potential. 



alternative parametrization of R-matrix theory has been 



developed by Brune 



( 2002 1 where all level shifts vanish 



and the partial widths and level energies are identical to 
the observed parameters. This approach simplifies the 
incorporation of known nuclear properties in the fitting 
procedure and the comparison with experimental reso- 
nance properties. 



2. Theory constraints: ab initio methods 

Ab initio methods - defined here as methods that pro- 
vide a quasi-exact solution to the many-body Schrodinger 
equation, such as the hyperspherical harmonic expan- 
sion (HH) and Green's function Monte Carlo (GFMC) 
methods, or that express observables in terms of a con- 
trolled expansion, such as effective field theory - play 
two critical roles. Two reactions discussed in this re- 



view, p-l-p — > d-t-e"* 



and ^He-l-p ^Re+e^+i^e 



presently beyond the reach of experiment. Thus we are 
entirely dependent on theory for the corresponding S- 
factors. The convincing demonstration that the rate for 
p+p — > d+e++i^e can be calculated to a precision of < 
1% is one of the important achievements of ab initio nu- 
clear theory, as described in Sec. |III[ 

Furthermore, ab initio methods potentially could be 
applied to all other reactions in the pp chain (and, far- 
ther in the future, to the CNO bi-cycle) to provide a 
more reliable basis for extrapolating data. One of the 
impressive examples of progress to date, the agreement 
between and data for d(p,7)'^He and theory (calculations 
employing variational HH wave functions in combination 
with an electromagnetic current operator with both one- 
and two-body components), is discussed in Sec. IV and 
illustrated in Fig. [3j 

Ab initio methods break into two broad categories, 
potential-based calculations and effective field theory ex- 
pansions. The former are distinguished from model- 
based methods discussed in Sec. |II.C.1| in two regards. 
First, they use a realistic interaction that fits two-body 
scattering data in detail, as well as certain bound-state 
properties of the lightest nuclei. Thus the interaction 
has both a rich operator structure and an explicit treat- 
ment of the short-distance repulsive core. Second, they 
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(18) 



with numerical techniques that can accurately treat an 
interaction of such complexity and with such disparate 
spatial scales, producing a quasi-exact solution of the 
many-body Schrodinger equation. The form of the three- 



body potential in Eq. (18), which contributes for A > 3 



but plays a less important role than the dominant two- 
body potential, is typically taken from theory. Once the 
wave functions are obtained, they can be combined with 
electroweak transition operators to produce estimates of 
observables. The transition operators include both one- 
body terms determined from the coupling of the single 
nucleon to the electroweak current, and two-body cor- 
rections, typically derived from one-boson-exchange di- 
agrams. Examples of the potential approach, including 
discussions of the associated issue of transition operators, 
are found in Sees. [111) |TVl and [7111 

The second approach is based on effective field the- 
ory (EFT). EFTs exploit the gap between the long- 
wavelength properties of nuclei that govern nuclear re- 
actions near threshold, and the short-range interactions 
in the NN potential that make an exact solution of the 
Schrodinger equation technically difficult. The calcula- 
tions are restricted to a limited basis describing the long- 
wavelength behavior, and the omitted degrees of freedom 
are absorbed into effective operators that can be orga- 
nized in powers of Q/Acut, where Q is the momentum 
characterizing the physics of interest and Acut is the mo- 
mentum characterizing the omitted physics. If carried 
out completely, no simplification is achieved, because the 
low-momentum EFT Lagrangian has an infinite number 
of such operators. EFT becomes useful when there is 
a significant gap between Q and Acut, so that only a 
small number of the effective operators corresponding to 
the leading powers in Q/Acut must be retained, to repro- 
duce long- wavelength observables to a specified accuracy. 
The coefficients of the leading operators can then be de- 
termined by fitting data: if enough constraints exist to 
fix all of the needed low-energy constants, then accurate 
predictions can be made about new processes. The appli- 
cation of this method to p-|-p — >■ d+e'^+Ve and ^He - I- p 



III 



— > ''He + c+ + j/g is described in some detail in Sees 

and |VII1 respectively. This approach can also be applied 
to d(p,7)3He. 

One of the potential-based methods now being devel- 
oped for reactions should be highlighted because of its 
established success in predicting bound-state properties 
throughout most of the Ip shell. The quantum Monte 
Carlo (QMC) approach combines the variational Monte 
Carlo (VMC) and GFMC methods ( |Pieper[ |2008| . The 
VMC calculation produces an approximate wave func- 
tion by minimizing the energy of a variational wave func- 
tion including elaborate two- and three-body correla- 
tions. The GFMC method is then employed to make the 
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needed small improvements to the VMC result required 
for a true solution to the Schrodinger equation. 

The GFMC method requires a local potential, so its 
use has been restricted t o the Argonne wis NN poten- 
tial (Wiringa et al. 1995), denoted AV18. There is also 



an important three-nucleon interaction, determined by 
fitting 17 bound- and narrow-state energies for A < 8 
dPieper et a"I||2001j ). The high quality of the QMC pre- 
dictions for energies of bound states and sharp resonances 
in nuclei with A < 12, and for charge radii, electromag- 
netic moments, and other observables, has been thor- 



oughly established (Pieper et al. 2001 2002 2004). 



Recent VMC-based calculations of capture cross sec- 



tions using realistic potentials ( 


Marcucci et al. 


2006 


Nol- 


lett 


2001 


Nollett et al. 


2001 


) represent a first step in 



extending the QMC program to reactions. These calcu- 
lations used VMC wave functions for bound states in '^H, 
■^He, ''He, ^Li, ^Li, and ^Be, as well as an exact deuteron. 
Initial states in the reactions d(a, 7)^Li, ^H(a, 7)^Li, and 
''He(a,7)^Be were computed as products of the reactant 
VMC wave functions and a correlation, matched to ex- 
perimental phase shifts, to describe the relative motion 
of the interacting nuclei. Work has focused, in particular, 
on building in the proper long-range clustering of the fi- 
nal states, as this is important in reproducing the proper 
energy dependence of S-factors. Results for ■^H(a,7)^Li 
closely match the measured absolute S-factor. However, 
the prediction for '^He(a,7)^Be lies below the data by 
about a factor of 1.3 to 1.45. 

Better QMC calculations of those and other cross sec- 
tions are possible. VMC wave functions were used partly 
because of the technical difficulty of computing quantities 
off-diagonal in the energy eigenstates using GFMC; this 
problem has now been solved, and electroweak matrix ele- 



ments between discrete levels have been computed ( Mar- 



cucci et al. 2008 Pervin et al. 2007). Scattering wave 



functions are also now being computed directly from the 
NN-I-NNN potential, with successful calculations of low- 
energy neutron-^He scattering wave functions reported 
by Nollett et al. (2007) using particle-in-a-box formula- 



tions of the QMC methods. 

While we have used the QMC approach to illus- 
trate the progress in quasi-exact approaches, there are 
other important efforts underway to compute cross sec- 
tions beyond A=4 from realistic NN potentials. Ex- 
amples include the ab initio no-core shell model both 
alone ( Navratil et al. . 2006a|b ) and in combination with 



the resonating group method (Quaglioni and Navratil 



2009); the Lorentz integral transform method (Efros 



et al.\ 2007); and the unita ry correlation operator method 
( |NefF and Feldmeier 2008 1 . The hypersherical harmonics 
method, which will be discussed in connection with the 
d(p,7)'^He and hep reactions, is also being extended to 
heavier systems. 

We anticipate that quasi-exact methods will soon be 
practical for many scattering and capture processes in 
light nuclei. Calculations based on exact solutions of ac- 
curate interactions will predict not only the energy de- 



pendences of solar fusion reactions but also absolute cross 
sections. Theory may thus provide a firm basis for vali- 
dating and extrapolating data and for resolving system- 
atic differences between measured data sets. 



3. Adopted procedures 

These are the procedures we adopt for fitting and ex- 
trapolating data: 

• In two cases, p-|-p d-|-e+-|-:^e and '^He-l-p — >■ 
*Hc+c++z/e, S-factor estimates depend entirely on 
theory. The goal in such cases should be the appli- 
cation of both potential and EFT or EFT-inspired 
methods, yielding consistent results with quantified 
uncertainties. As detailed in Sec. |III[ one is close 
to achieving this for Sn, with two methods provid- 
ing consistent answers and uncertainties of < 1%, 
and with a third method (EFT) potentially reach- 
ing similar precision, if ancillary measurements can 
better determine the needed low-energy constant. 
In the case of Shop, a less critical cross section, the 
further developments of methods like Green's func- 
tion Monte Carlo will provide an important check 
on the current state-of-the-art, a variational calcu- 
lation in which a correlated hyperspherical harmon- 
ics expansion was used. 

• In cases where data exist through the energy range 
of astrophysical interest, much can be done inde- 
pendent of theory. A polynomial representation of 
S{E), e.g., values for S(0), S'(0), and S"(0), could 
be obtained by directly fitting the data. 

However, as S(£') represents the bare cross sec- 
tion, theory may still be needed to remove the 
effects of screening in the terrestrial target. As 
detailed above, there is some confidence that the- 
ory determines the functional form of the screening 
(Eq. ([Is])), so that such effects can be subtracted 
given sufficient low-energy data to fix the numeri- 
cal value of the screening potential (which theory 
appears to predict less reliably). This issue arises 
in S33. 

• In cases where data exist but are not adequate to 
fully characterize the cross section in the region of 
astrophysical interest, we advocate the use of fitting 
functions motivated by theory to extrapolate data, 
with data determining the normalization. To the 
extent that well-justified models differ in their pre- 
dictions, additional uncertainties must be assigned 
to S(0) and its derivatives. Judgment is required 
in assessing the models and determining how they 
should be applied, e.g., the range in E over which 
a given model is likely to be valid. Each work- 
ing group was asked to consider such issues, and 
to present and justify the procedures it followed to 
assess associated fitting uncertainties. 
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D. Treatment of uncertainties 

The treatment of uncertainties - the statistical and 
systematic errors in data and the impact of imperfect 
theory in fitting and extrapolating data - is discussed 
in some detail in the Appendix. There are cases where 
several high-quality data sets exist, each with errors that 
presumably reflect both the statistical and evaluated sys- 
tematic uncertainties of the experiment, that disagree by 
more than the error bars would indicate. In treating such 
cases, an error-bar "inflation factor" is commonly intro- 
duced, to account for the apparent underestimation of 
systematic errors. We have done so following Particle 



Data Group (PDG) conventions (Amsler et al. 2008), 



with one minor modification described in the Appendix. 
Uncertainties quoted in this review correspond to one 
standard deviation (68% confidence level). 

As discussed in the Appendix, there are alternative 
prescriptions for apportioning the unidentified systemat- 
ics - and thus the inflations - among the experiments 
that disagree. However our group concluded that the 
PDG procedure was the best choice both for technical 
reasons and because the procedure is widely used in the 
physics community. 



III. THE pp REACTION 

The rate for the initial reaction in the pp chain, p-l-p— 
d + e+ + i^e, is too small to be measured in the labora- 
tory. Instead, this cross section must be calculated from 
standard weak interaction theory. 

As in Solar Fusion I, the adopted value and range for 
the logarithmic derivative is taken from Bahcall and May] 
(fT969l), 



S'ii(O) = Sii(O) (11.2 ±0.1) MeV" 



(19) 



This result is in excellent agreement with those obtained 
from linear fits to the modern potential-model calcula- 
tions of Schiavilla et al. (1998), which yield values of 



11.14 MeV"^ and 11.16 MeV^for the full and impulse- 
approximation calculations. As the Gamow peak energy 
is ^ 6 keV for temperatures characteristic of the Sun's 
center, the linear term generates a < 8% correction to 
the E = value. The 1% uncertainty in Eq. ( 19 ) cor- 



responds to a < 0.1% uncertainty in the total reaction 
rate. This is negligible compared to other uncertainties 
described below. Therefore, in the following, we focus on 
Sii(O). 

At zero relative ene rgy Sii(O) can be written (Bahcall 
and Mayl 119681 119691), 



Sii(O) = 6n'^mpa\n2 



a' 



ryO 



Ga 

Gv 



fR 



(ft) 



0+-+0+ 



(20) 



where a is the fine-structure constant; rrip is the proton 
mass; Gy and Ga are the usual Fermi and axial- vector 
weak coupling constants; 7 = (2/ii?^)^/^ = 0.23161 fm~^ 



is the deuteron binding wave number; /i is the proton- 
neutron reduced mass; Bd is the deuteron binding en- 
ergy; is the phase-space factor for the pp reaction 
with radiative corrections; (/t)o+-s.o+ is the ft value for 
superallowed 0+ — > 0^ transitions; and A is proportional 
to the transition matrix element connecting the pp and 
deuteron states. 

Inserting the current best values, we find 



Sii(O) = 4.01 X 10-25 MeV b 



Ga/Gv 
1.2695 




(21) 



We now discuss the best estimates and the uncertainties 



for each of the factors appearing in Eq. (21). 

We take (/i)o+^o+ = (3071.4 ± 0.8) s, the value for 
superallowed 0^ — 0^ transitions that has been deter- 
mined from a comprehensive analysis of experimental 
rates corrected for radiative and Coulomb effects ( Hardy] 



and Towner 20091. This value determines the weak 
mixing matrix eleme nt \Vud\ = 0. 97418(27), the value 
adopted by the PDG ( [Amsler et aH|2008| ). This ft value 



is also consistent with (3073.1±3.1) s used in Solar Fusion 
I. 

For Ga/Gv, we use the PDG value Ga/Gv = 1.2695± 
0.0029 which is consistent with 1.2654 ± 0.0042 used in 
Solar Fusion I. 

For the phase-space factor we have ta ken the value 
without radiative corrections, fpp = 0.142 (Bahcall and 
1969) and increased it by 1.62% to take into ac- 



May 



count radiative corrections to the cross section ( Kurylov 



et al.. 2003). The main source of error is from neglected 



diagrams in which the lepton exchanges a weak boson 
and a photon with different nucleons. These diagrams 
are estimated to modify fpp by ^ 0.1%, based on scaling 
the similar nucleus-dependent correction in superallowed 
(3 decay (Kurylov et al. 2003[ ). It would be useful to 
check this estimate through direct computations. We 
adopt f^p = 0.144(1 ± 0.001), which is consistent with 
0.144(1 ± 0.005) used in Solar Fusion I. 

The dominant uncertainty in Sii(O) comes from the 
normalized Gamow- Teller (GT) matrix element A. A 
great deal of theoretical work since Solar Fusion I has 
focused on reducing this uncertainty. In Solar Fusion I A 
was decomposed into A = A{1 + S), where A represents 
the contribution of the one-body transition operator and 
A6 that from two-body corrections. A thus involves an 
evaluation of the Gamow- Teller operator between the 
initial-state pp wave function and the final-state deuteron 
wave function. A^ = 6.92(1 ± 0.002l°;°J^) was adopted, 
where the first and second uncertainties reflect, respec- 
tively, variations in empirical values of the deuteron and 
low-energy pp scattering parameters, and the model de- 



pendence of the nuclear potential ( Kamionkowski and 



Bahcall 1994). The value and uncertainty of the ex- 
change current contribution, 6 = O.OIIq'qi; "^^-s deter- 
mined from the range of values of published calculations. 



15 



following the conservative recommendation of Bahcall 
and Pinsonneault (1992). 



mate 



Two major steps have contributed to reducing the un- 
certainty on A since Solar Fusion I. The first is a much 
deeper understanding of the correlation between the un- 
certainties in A and SA: the overall uncertainty in A 
can be described by a universal parameter that can be 
fixed by a single measurement. The study of Schiavilla 
et al. ( 1998 ) demonstrated this phenomenologically in the 



context of potential-model approaches, while later analy- 
sis via EFT provided a more formal justification (Butler 



et al 2001 



Park et al. 2003|). The second step is the 

as 



use of the precisely known tritium /3 decay rate FT 



first proposed by Carlson et al. (1991), to fix this uni- 
versal parameter. This has been done in both potential 
models (Schi avilla et al.\ |1998[ ) and in the hybrid EFT 
approach ( Park et al. 2003p . We briefly describe these 
developments. 



A. Progress in potential models 

The most elaborate calculation for the pp fusion pro- 
cess in the potential-model approach (see Sec. II. C. 2) 



was carried out by Schiavilla et al. ( 1998 ) . A comparison 



of the results for five representative modern potentials 
- potentials designed to accurately reproduce nucleon- 
nucleon scattering data - yielded A^ = 6.975 ± 0.010. 
This study demonstrated the importance of using the 
tritium /3 decay rate to constrain the two-body GT tran- 
sition operator. (Both the Fermi and GT operators con- 
tribute to tritium /? decay, but the former can be reliably 
calculated because of the conserved vector current and 
the smallness of isospin breaking effects, ^ 0.06%.) If 
one adjusts the uncertain strength of the exchange cur- 
rent so that the tritium /3 decay rate is reproduced, the 
variation in Sii(O) that otherwise would come from the 
choice of the phenomenological potential is largely re- 
moved. Predictions for five representative high-precision 
phenomenol ogical potentials fall in a narrow interval 7.03 



<A <7.04 (Schiavilla aZ. 



1998) 



We note two other sources of model dependence that 
contribute to the overall uncertainty in A. First, as three- 
body potentials and currents contribute to the tritium 
/3 decay rate, uncertainties in modeling such effects will 
influence the extracted constraint on the two-body cur- 
rents needed for Sii(O). The best estimate of the con- 
sequences of this uncertainty for Sii(O), ~ 0.8%, comes 
from the chiral (or pionful) EFT* approach described 
below. Second, the experimental uncertainties in the ef- 
fective range parameters for nucleon-nucleon scattering 
will propagate to A. We have assigned a 0.5% uncer- 

2 

tainty in A to this source, pending future work in EFT 
to better quantify this uncertainty. By adding in quadra- 
ture these uncertainties of 0.8% and 0.5% and the smaller 
uncertainty associated with the above potential range, 

2 

A = 7.035 ± 0.005, we obtain the potential model esti- 



A = 7.035(1 ±0.009). 



(22) 



B. Progress in effective field theory (EFT) 



II.C.2 to 



The application of EFT, described in Sec. 
the calculation of the pp fusion rate (and several other 
electroweak processes in light nuclei) is one of the notable 
developments since Solar Fusion I. There have been two 
lines of EFT calculations of pp fusion, described below. 



1. Hybrid EFT (EFT*) 



Electroweak nuclear transitions in EFT 



EFT, ^EFT 



(23) 



require initial and final nuclear wave functions and the 
transition operators to be derived from EFT. However, 
this has not yet been achieved in EFT with dynamical 
pions for pp fusion. Instead, a hybrid approach ( [Park 



et al. 2003) called EFT* (or MEEFT) has been devel- 



oped in which transition operators are taken from chi- 
ral perturbation theory (xPT), but sandwiched between 
phenomenological wave functions, \I>p^°" and vJ/P'^''"^ gen- 
erated by a potential model. As discussed below, this 
approach is a substantial improvement over the earlier 



calculation of Park et al. ( 1998 ) 



For the low-energy GT transition that governs pp fu- 
sion, the one-body transition operators Of^"^ are well 
known, while the two-body operators Of:^"^ contain only 
one unknown low-energy constant (LEG). This LEG, de- 
noted by d^, parameterizes the strength of contact-type 
four-n ucleon coupling to the axial current. 



Park et al.\ 



(2003) chose to determine from the tritium /3-decay 
rate Fj. The fact that vj/ph™ jg ^gt exactly an eigenstate 
of the EFT Hamiltonian can in principle be a source of 
concern, but it is plausible that the mismatch affects pri- 
marily the short-distance behavior of the wave function, 
so that the procedure of fixing the relevant LEC(s) to 
data can remove most of the inconsistency: While C^pt 
by construction is valid only well below Aqcd , the use of 
the phenomenological Hamiltonian, Eq. (18), introduces 



high momentum components above Aqcd- To test this 
procedure, one can introduce a cutoff Ann to eliminate 
high-momentum components in the two-nucleon relative 
wave function, fitting the LEG as a function of this pa- 
rameter. One expects, if the fitting of the LEG reason- 
ably accounts for missing or inconsistent short-distance 
physics, little Ann dependence would be found in the cal- 
culated pp fusion rate. The residual dependence on Ann, 
when this cutoff is varied over a physically reasonable 
range, provides a measure of the model independence of 
an EFT* calculation. 
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The Park et al. (2003) calculation included up to 
next-to-next-to-next-to-leading order (N-^LO) terms in 
chiral expansion, and after fitting to Fj, yielded 

2 

A — 7.03(1 ± 0.008). The uncertainty was estimated 

2 

from the changes in A when Ann is varied over an en- 
ergy range typical of vector meson masses, 500 to 800 
MeV. A rough estimate based on higher order chiral 
contributions was also made. Specifically, the contribu- 
tions of the first four chiral orders to A follow the pat- 
tern (l+0.0%+0.1%-|-0.9%), while the fifth-order term is 
estimated to be ~ 0.4%. We assume that the second- 
and third-order terms are accidentally small, while the 
fourth- and fifth-order terms reflect the convergence of 
the expansion in ttIt^ j Kqcd ^1/7. Three-body currents 
contribute in sixth order. We therefore use the size of the 
fifth-order term, 0.4%, as a measure of the uncertainty 
due to neglected higher order contributions (including 
three-body currents). 

Full EFT calculations that use "^^^^ instead of ^-p^^™, 
thus eliminating operator-wave function inconsistencies, 
are an important goal. Progress toward this goal includes 
recent constructions of EFT-based nuclear interactions; 
see, e.g., |Epelbaum (2006) and Gazit et al. (2009). 



The experime ntal goal is to con strain A to < 1. 5% for 
pionless EF T ( |Chen aL[[2005l ) and chiral EFT* ([AUdo 
eiaLl 120021). 



3. Comment on Mosconi et 



Mosconi et al. { 2007 1 have compared v-d reaction cross 
sections for various models that differ in their treatments 
of two-body transition operators, concluding from this 
comparison that the results obtained in potential mod- 
els, EFT*, and pionless EFT have uncertainties as large 
as 2-3%. Although they address only ly-d cross sections, 
a comment is in order here because this process is closely 
related to that for pp fusion. Mosconi et al. (2007) 



reach their conclusions by examining the scatter of un- 
constrained calculations of the v-d cross section. How- 
ever, all state-of-the-art calculations use F^ to reduce 
two-body current and other uncertainties, as we have de- 
tailed here. Once this requirement is imposed, the scatter 
in the calculated value of I'-d cross sections is significantly 
reduced. 



C. Summary 



2. Pionless EFT 

This approach can be applied to processes where the 
characteristic momentum p is much smaller than the pion 



1999 


Chen et al. 


1999 


Ka- 



mass TOtt (Bedaque et al. 
plan et al. 1996[ ), which is the case for solar pp fusion. 



Pions can then be integrated out, so that all nucleon- 
nucleon interactions and two-body currents are described 
by point-like contact interactions with a systematic ex- 
pansion in powers oi p/mj^. The one- and two-body con- 
tributions individually depend on the momentum cut-off 
but the sum does not. Thus, A and AS in pp fusion 
are correlated. In pionless EFT only one two-body cur- 
rent (with coupling Li^a) is needed in the description of 
deuteron weak breakup processes, through next-to-next 
to-leading o rder (NNLO) in the p/m^ expansion (Butler 
et al. 2001). This two-body current is a GT operator 



Other two-body currents are either missing due to conser- 
vation of the vector current, or involve matrix elements 
suppressed due to the pseudo-orthogonality of the initial- 
and final-state wave functions. This means the universal 
number Li^a encodes the two-body contributions for all 
low-energy weak deuteron breakup processes, so that a 
single measurement will fix the rates of all such processes. 
The other approaches discussed above share this feature. 
The computation of A in pionless EFT was carried 



out to the second order by Kong and Ravndal (2001) 



and Ando et al. (20081 and then to the fifth order by 
Butler and Chen (|2001 ). Constraints on Li^a from two 



nucleon systems (Butler et al. 

2 

yield A = 6.99±0.21. The MuSun experiment ( Andreev 



2002} |Chen ^Tal] |2003t 



et al. 2008 ) is taking data on fi capture on deuterium. 



We have seen that the various approaches discussed 

2 

above yield accurate and very consistent values for A . 
The remaining factors in Eq (18) also have uncertainties, 
but these are common to all the calculations. Adding all 
the uncertainties in quadrature, we find that the current 
best estimates for Sii(O) are 

4.01(1 ± 0.009) X 10"^^ MeV b potential models 

4.01(1 ± 0.009) X IQ-^^ MeV b EFT* 

3.99(1 ± 0.030) X IQ-^^ MeV b pionless EFT. (24) 

The larger uncertainty in the pionless EFT result is due 
to the relatively weak constraints on Li a that can be 
imposed within two-nucleon systems but, as mentioned, 
this situation will soon be improved. The agreement of 
the central values obtained in the potential model and 
EFT* indicates the robustness of the results as long as 
the two-body current is constrained by tritium f3 decay. 
Meanwhile, the agreement of the error estimates in the 
two approaches is primarily due to the fact that, as ex- 
plained above, the dominant part of the uncertainty has 
been estimated using the same argument. Based on the 
result obtained in the potential model and EFT*, we 
adopt as the recommended value 



Sii(O) = 4.01(1 ± 0.009) X 10"^^ MeV b. 



(25) 



We adopt the |Bahcall and May| ( |1969[ ) value for S[^{0) 
S[^{0) = Sii(0)(11.2±0.1)MeV'^ , (26) 



Bahcall and May ( 1969| ) also estimated dimensionally 
that (0) would enter at the level of ~ 1%, for temper- 
atures characteristic of the solar center. As this is now 
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comparable to the overall error in Sn, we recommend 
that a modern calculation of S"]^(0) be undertaken. 



IV. THE d(p,7)^He RADIATIVE CAPTURE REACTION 

The radiative capture of protons on deuterium is the 
second reaction occurring in the pp chain. Because this 
reaction is so much faster than the pp weak rate dis- 
cussed in the previous section, it effectively instanta- 
neously converts deuterium to ^He, with no observable 
signature. Thus uncertainties in its rate have no conse- 
quences for solar energy generation. By comparing the 
pp and d(p,7)'^He rates, one finds that the lifetime of a 
deuterium nucleus in the solar core is ~ 1 s, and that 
the equilibrium abundance of deuterium relative to H is 
maintained at ~ 3 x 10^^^. 

However, the d(p,7)'^He reaction plays a more promi- 
nent role in the evolution of protostars. As a cloud of 
interstellar gas collapses on itself, the gas temperature 
rises to the point of d(p,7)'^He ignition, ~ 10^ K. The 
main effect of the onset of deuterium burning is to slow 
down the contraction and, in turn, the heating. As a 
consequence, the lifetime of the proto-star increases and 
its observational properties (surface luminosity and tem- 
perature) are frozen until the original deuterium is fully 
consumed ( jStahler] 1988[ ). Due to the slow evolutionary 
timescale, a large fraction of observed proto-stars are in 
the d-burning phase, while only a few are found in the 
earlier, cooler, rapidly evolving phase. A reliable knowl- 
edge of the rate of d(p,7)'^He down to a few keV (the 
Gamow peak in a proto-star) is of fundamental impor- 
tance for modeling proto-stellar evolution. 

The pd reaction also plays an important role in Big 
Bang nucleosynthesis, which begins when the early uni- 
verse has cooled to a temperature of ~ 100 keV. The un- 
certainty in the pd reaction in the relevant energy win- 
dow (25-120 keV) propagates into uncertainties in the 
deuterium, ■^He and ''Li abundances, scaling roughly as 



^He 
IT 



H 



Of 

oc Ytp^ , 



(27) 



w here Rpd is the value of S12 relative to the fiducial value 



Gyburt (2004). Thus a 10% error in the pd capture 
rate propagates into roughly 3.2%, 3.8% and 5.9% un- 
certainties in the light element primordial abundances, 
d, ^He and ''Li, respectively. 



A. Data sets 

The extensive experimental data sets for pd radiative 
capture include total cross sections and spin polarization 
observables at center-of-mass energies E ranging from 
several tens of MeV to a few keV, covering all the rele- 
vant astrophysical energies. In the regime i? < 2 MeV 
(below the deuteron breakup threshold) , the relevant ex- 
perimental data include Bailey et al. { 1970[ ); Gasella et al. 




E [MeV] 



FIG. 3 (Color online) The astrophysical Si2-factor 



datasets (Casella et al. 


2002 Griffiths et al.\ 1962 |Ma et al. 


1997 Schmid et al. 


19961 are plotted together with theoret- 



ical predictions of Marcucci et al. ( |2005| . The solid line rep- 
resents the "full" theoretical calculation, while the red band 
represents the 68% lower and upper bounds of the adopted 
quadratic best fit to the four experimental datasets (see text 
and Eq. (29 1 for more explanation). In the insert, the S12- 
factor of the ^II(p,7)^IIe reaction in the energy range 0-50 
keV, obtained with the Argonne vis two-nucleon and Urbana 
IX three-nucleon Hamiltonian model in the impulse approx- 
imation (dashed line) and with inclusion of interaction cur- 
rents (solid line), is compared with the experimental results. 



(2002 


); Griffiths et al. ( 


1963 




1962D; Ma et al. 


(19971; 


Schm 


id et al. (|1995 


1996 


). The 


IGrifHths et al. "p963 ) and 



Bailey et al. ( 1970 1 low energy data may be 



high because of the use of incorrect stopping powers (Ma 



et al. 



et al. 



1997; Sc hmid et a^.[|1995[|1996| . Also, the [Schmid 
( |1995( ^1996 1 data sets may have not propagated 



their energy-dependent systematic uncertainties. In Fig. 
|3) the data for S12 used for the best fit in Sec. |IV.G|are 
plotted toge ther with theoretical predictions of [Marcucci] 
\et al. ( 2005 ) . The observed linear dependence of S12 on 



at low energies, as well as the angular distributions of the 
cross section and polarization observables, indicate that 
the d(p,7)'^IIe reaction proceeds predominantly through 
s- and p-wave capture, induced, respectively, by magnetic 
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(Ml) and electric (El) dipole transitions. The Ml tran- 
sitions (proceeding through and ^83/2 pd channels) 
are especially interesting, as the one-body Ml operator 
cannot connect the main s-state components of the pd 
and '^He wave functions at low energies. Because of this 
"pseudo-orthogonality" only the small components of the 
wave functions contribute in the impulse approximation 
(lA). In contrast, as exchange-current operators are not 
similarly hindered, their matrix elements are exception- 
ally large relative to those obtained with the one-body 
Ml operator. The suppression of matrix elements cal- 
culated in the lA and their consequent enhancement by 
exchange-current contributions are a feature common to 
other Ml-induced processes in A=3 and 4 systems, such 
as the nd and n^He radiative captures at thermal neutron 
energies. 



B. Theoretical studies 

The most extensive and recent theoretical studies of 
the d( p,7)^He reaction at low energies have been carried 
out by Marcucci et al. (2005). The calculated S12, shown 
in Fig. 3j is in excellent agreement with data. To describe 
the pd continuum and ^He bound states, these authors 
used variational wave functions built in a correlated- 
hyperspherical-harmonics (CHH) basis for a Hamiltonian 



uncertainties have been underestimated. Serpico et al. 



consisting of the Argonne vis two-nucleon ( Wiringa et al. 
[T995) and the Urbana IX (jPudliner et al.\ |1995|) three- 



nucleon potentials. This Hamiltonian is known to repro- 
duce a variety of three-nucleon bound- and scattering- 
state properties, including binding energies, charge radii, 
elastic and inelastic cross sections, and low-energy polar- 
ization observables, while the accuracy of the CHH varia- 
tional method is comparable to that of other quasi-exact 



methods (Nogga et al. 2003) 



The nuclear electromagnetic current consists of one- 
body terms (the lA currents), originating from the con- 
vection and spin-magnetization currents of individual 
protons and neutrons, and two- and three-body exchange 
currents, constructed from the corresponding potentials 
by a procedure that preserves current conservation (CC). 
The method by which this is achieved has been improved 



over the years (Riska 1984 Schiavilla et al. 19981, and 
its latest impleme ntation is discussed at length by |Mar-| 



cucci et al. (2005). The currents are still model depen- 
dent, of course, as CC places no constraints on their 
transverse components. 

The calculated value for 812(0) including exchange- 
current contributions is 0.219 eV b, in excellent agree- 
ment with the value extrapolated from the LUNA mea- 



surements (0.216± 0.010 eV b), and evaluations by Cy- 



burtU2004) (0.227 ± 0.014 eV b) , [Descouvemont al. 



p004n 0-223± 0.007 e V b) and [Serpico et al.\ (|2Q04p 
(a2T4± 0.007 eV b). In |Descouveinont et aL| ( |2004 D sys- 
tematic and statistical errors are combined before follow- 
ing a standard fitting procedure. However, as this artifi- 
cially reduces the impact of systematic errors, their cited 



(2004) properly separates systematic and statistical er- 
rors in their treatment, but do not cite 68% confidence 
limits, also yielding an error that is too small. The eval- 
uation by Cyburt ( 2004 1 separates systematic and statis- 



tical uncertainties and cites errors consistent with 68% 
confidence limits, yielding realistic uncertainties. 



C. Summary 



evaluate the 


Casella et al. 


(20021, 


Schmid et al 


(1996 


) and 


Ma et al. 



(1997p data, determining Si2{E) as a function of the 



center-of-mass energy by fitting the four data sets by a 
quadratic polynomial in E. We adopt this fitting proce- 
dure, despite our earlier arguments favoring fitting for- 
mulas that are motivated by theory, because the energy 
window of interest is fully covered by the experiments. 
This yields 



Si2(0)=0.214t°:°i^ eV b. 



(28) 



in agreement with previous evaluations. The error is 
larger here, because of the exclusion of the [Bailey et al. 
(1 19701) data. 



We also determined the 68% upper and lower bounds 
for the quadratic parameterizations, valid for E < 1 
MeV, the range spanned by the data we considered. The 
results are (see also Fig. [3]) 



Slower 
1 9 



SiTiE) = 0.2145 + 5.5612 



0.1983 + 5.3636 



2.9647 



— y 

MeV / 



4.6581 



E 



s^pp"(^) 



^MeV 
0.2316 + 5.7381 



6.5846 



— y 

MeV/ 



E 
MeV 

eV b 

E 
MeV 

eV b 

MeV J 
eV b. 



(29) 



The results determine the S-factor and its uncertainty in 
the vicinity of the solar Gamow peak. In particular, for 
a temperature characteristic of the Sun's center, 1.55 x 
10^ K, 

812(^^0 = 6.64 keV) = 0.252 ± 0.018 eV b, (30) 
so that the estimate uncertainty is ~ 7.1%. 

V. THE =^He(^He,2p)*He REACTION 

The '^He('^He,2p)''He reaction is the termination of the 
ppl cycle and thus, as Solar Fusion I describes in more 
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detail, uncertainties in this cross section played a promi- 
nent role in early speculations about a nuclear astro- 
physics solution to the solar neutrino problem. As an 
increase in S^slE) would reduce the branchings to the 
ppll and ppIII cycles ~ thus also reducing the neutrino 
fluxes measured by Davis - the possibility of an undis- 
covered narrow resonance at energies beyond the reach 
of early experiments was raised by |Fetisov and Kopy-| 
( 1972[ ) and |Fowler| ( |1972[ ). This motivated efforts to 
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measure S33(ii^) at lower energies, and particularly stim- 
ulated the efforts of the LUNA collaboration in the 1990s 



to map the cross section in the solar Gamow peak ( 


Arpe- 


sella et al. 1996 Bonetti et al\ 1999 Greife et aLl|1994| 


Junker et al. 


1998 


). The principal result since Solar Fu- 



( 1999 ) , extending measurements to the lower edge of the 
Gamow peak at 16 keV, making S33(i?) the most directly 
constrained S-factor within the pp chain. 

S33 {E) remains of significant importance, as it controls 
the ppl/ppll-|-pplll branching ratio and thus the ratio 
of the pp/pep to ''Be/^B neutrino fluxes. This ratio is 
important to future strategies to better constrain neu- 
trino oscillation parameters and matter effects, through 
comparison of high-energy (matter influenced) and low- 
energy (vacuum) fluxes. The ratio of S33 to S34 enters 
in computing the neutrino energy losses of the Sun, and 
thus influences the connection between the Sun's photon 
luminosity and its total energy production. 



A. Data sets and fitting 



We consider data available at the time of Solar Fusion 



I (Arpesella et al. 


1996 


Bacher and Tombrello 


1965 


Dwarakanath and 


Winkler 1971| Greife et al. 


1994 



Junker et al. 1998 Krauss et al. 1987) as well two new 



data sets: the extreme low energy data of LUNA ( Bonetti 



et al. 19991 and results from the OCEAN experiment 



(Kudomi et al. 2004) at energies shghtly above the so- 
lar Gamow region. In order to follow the recommended 
fltting prescription discussed in the Appendix, one needs 
a detailed discussion of systematic uncertainties, partic- 
ularly common mode systematics. This requirement re- 
duces the datasets considered to just four experiments. 
The earliest of these originates from the Muenster group 



cations 



(Krauss et al. 1987), followed by the two LUNA publi- 



Junker et al. ( 1998 ) (which supersedes Arpesella 



et al.\ (|1996D) and|Bonetti et aL|(|1999[); and the OC EAN 
effort [Kudomi et aL| ( |2004[ ). |Krauss et al.\ ( |1987[ ) and 
Kudomi et al. ( 2004 1 identified a common systematic er- 
ror for their respective data sets while the LUNA group 
provided statistical and systematical errors at each ex- 
perimental energy measured. In order to use a uniform 
treatment we calculated an average systematic error for 
the latter data sets. Larger systematic errors were noted 
only at the lowest energies (due to uncertainties in stop- 
ping power) where the total error is dominated by statis- 
tics. 



Past efforts have fit data to an S-factor including 
screening corrections, with the bare S-factor a polyno- 
mial up to quadratic order. 



S33(i?) 



S^3-^(i?)exp 



E 



(31) 



Sbarc 
33 



{E) = S33(o) + s^3(o)£; + -s;(3(o)i?2. 



Although model calcu l ations of S3f ''^ {E) are available 
(see, e.g., Typel et al. (1991)), a phenomenological rep- 



resentation for the bare S-factor is appropriate because 
the data extend to the Gamow peak. There is no need 
for a theoretical model to guide an extrapolation, apart 
from the functional form of the screening potential. 

The selected data for this review cover the range from 
the solar Gamow peak to 350 keV, providing a limited 
range with which to perform a four parameter fit to 
the S-factor including electron screening (S33(0), S33(0), 
^33(0), Ue). We test the robustness of the fit parame- 
ters, by varying the order of the polynomial for the bare 
S-factor. Our results are in Table ITll 



TABLE II Table of fit parameters and their total errors for 
constant, linear, and quadratic representations of the bare 
S-factor. 



parameter 


constant 


linear 


quadratic 


833(0) (MeV b) 


4.84 ±0.13 


4.95 ±0.15 


5.32 ±0.23 


^'xM (b) 


N.A. 


-1.06 ±0.51 


-6.44 ± 1.29 


S^,'3(0) (MeV-^ b) 


N.A. 


N.A. 


30.7 ± 12.2 


(eV) 


395 ± 50 


360 ± 55 


280 ± 70 




35.4 


34.1 


31.8 


xlot/dof 


0.40 


0.39 


0.37 



Our quadratic fit agrees quite well with the fit derived 



by Krauss et al. ( 1987 ) , adopted in the r eaction rate com- 
pilation of Caughlan and Fowler (1988). However, there 



is a significant spread in fit parameter values for the dif- 
ferent order polynomial fits, with slight decreases in the 
total x^- One can also see this spread in fit results from 



other groups ( 


Bonetti et al. 


Kudomi et al.\ 


2004 


). This su 



1999 Junker et al. 



1998 



have the resolving power to accurately determine all fit 
parameters: there are strong correlations for the choices 
of data and fitting functions made here. Adopting any 
single fit will underestimate the uncertainties due to the 
degeneracy between parameter values. From Bayes's the- 
orem, assuming that the S-factor in this region {E < 350 
keV) can be described without cubic terms, we can de- 
rive constraints on the parameters by weighting each fit 
in Table [n] by its total value. This method takes into 
account the spread from fit-to-fit. We find 



S33(0) 



5.21 ±0.27 MeV b 



(32) 
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FIG. 4 (Color online) The data, the best quadratic+screening 
result for S33(i?), and the deduced best quadratic fit (red 
line) and allowed range (yellow band) for S^t"- See text for 
references. 



S'33(0) = 

S:,'3(0) - 

Ue = 



-4.90 ±3.18 b 
22.4 ± 17.1 MeV" 
305 ± 90 eV. 



The results reveal that existing data cannot strongly 
constrain all of the fitting parameters separately, and in 
particular do not sharply constrain Ue- To improve con- 
straints on the screening potential one will need more pre- 
cise data from near the Gamow peak, as well as new mea- 
surements up to the MeV range (with well documented 
systematics) to better determine the higher-order terms 
in the quadratic fit. New theory efforts in determining 
the shape of this S-factor would also be beneficial, as new 
low energy '^He-'^He elastic scattering data could be used 
as an additional constraint. 

However, our principal concern is the precision with 
which 83!''° (i?) can be determined in the vicinity of the 
Gamow peak, not the separate parameters. From the 
fit's correlation matrix we find 

S!^f *(£;) = 5.21 - 4.90 ( + 11.21 ( MeV b 



McV 



VMeV J 



0.075 - 1.516 ( ) + 14.037 1 1 

VMevy VMevy 



E y 



-15.504 



E 
MeV 



71.640 

VMevy 



1/2 



where 



S',riE)^STiE)±5Sss{E). 



MeV b 



(33) 



Because these results were obtained with a phenomeno- 
logical fitting function, their reliability has been establish 



only for the energy range covered by the data employed 
in the fit. Thus Eq. (33 1 should be used for energies E < 
350 keV. For a temperature 1.55 x 10^ K corresponding 
to the Sun's center, we find at the Gamow peak 

S^3''°{Eq = 21.94 keV) 5.11 ± 0.22 MeV b, (34) 

so that the estimated uncertainty is 4.3%. 



VI. THE 3He(a,7)^Be REACTION 

When Solar Fusion I appeared, the most recent 
•^He(^He,7)^Be measurement was 10 years old. The four 
new measurements that have been published since that 
time, in response to a challenge by John Bahcall, are the 
focus of this section. 

For energies of interest, E < 1 MeV, ■^He(''He,7) ''Be 
is a nonresonant reaction, predominantly external direct 



capture (Christy and Duck 1961 1 by electric dipole emis- 



sion from s- and d-wave initial states to the two bound 
states of ^Be. Reaction measurements have been made 
by detecting the prompt 7-rays, the ''Be activity, and 
the ^Be recoils. Below we discuss the measurements, the 
theory needed to extrapolate the measurements to astro- 
physical energies, and our determination of 834(0). 



A. Experimental measurements 



Groups at the Weizmann Institute ( Singh et al.\ 2004 1 
and at the University of Washington-Seattle (IBrown 



et al.. 2007) carried out cross section measurements in 
the center-of-mass energy range E = 0.42 to 0.95 MeV 
and 0.33 to 1.23 MeV, respectively, using gas cells with 



Ni entrance windows. The LUNA collaboration ( 


Bem- 


merer et al. 2006a Confortola et al. 2007| |Gyiirky et al. 


2007 


) (see also Gostantini et al. 


(20081) carried out low- 



background measurements from E — 0.093 to 0.170 MeV 
at the LUNA facility in the Gran Sasso underground lab- 



oratory, and a European collaboration (Di Leva et al. 



2009) (here called ERNA) made measurements from E 
= 0.65 to 2.51 MeV, both with windowless gas cells. 

An important concern in Solar Fusion I was whether 
•^IIe(^He,7)^Be measurements made by detecting the 
^Be activity might be affected by background ''Be pro- 
duced by contaminant reactions. Possibilities include 
^Li(d,n)^Be or ^°B(p,Q;) ''Be, which could occur given pro- 
ton or deuteron contamination in the '''He beam in com- 
bination with ^Li or ^'^B contamination in the gas cell, for 
example, in the foil or beam stop. Only one of the older 
experiments - that of Osborne - involved measurements 
of both prompt 7s and ^Be activity (see Solar Fusion 
I for older references). While the Osborne experiment 
found agreement between the ■^He(^He,7)^Be cross sec- 
tions determined by the two methods, in general the cross 
section determined from activity-based experiments was 
somewhat larger than that determined from prompt-7 
experiments. 
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In the new experiments, all but the Weizmann group 
measured both prompt 7s and ^Be activity, while ERNA 
also measured ^Be recoils. In each of these experiments, 
the cross sections deduced by the different methods were 
consistent, leading to upper limits on nonradiative cap- 
ture of 2-5% from E = 0.09 to 2.5 McV. This is consistent 
with theoretical calculations that indicate much smaller 
rates expected for EO capture and other electromagnetic 
processes that could pro duce ''Be without accom pany- 
ing energetic prompt 7s (Snover and Kurd 2003). All 



new experiments except that of the Weizmann group em- 
ployed "^He beams and "^He targets, thus minimizing po- 
tential problems with background ^Be production. In the 
new experiments sensitive checks ruled out contaminant 
^Be production at lower levels. Thus we see no reason to 
doubt the new activity measurements. 

^Be activity measurements provide a direct determi- 
nation of the total cross section. In contrast, as prompt 
7-ray yields are anisotropic, one must take into account 
detector geometry and the anisotropy to determine a to- 
tal cross section. [The ^ 30% capture branch to the 429- 
keV first excited state of ^Be has usually been determined 
from the isotropic 429 keV -> ground state yield.] Un- 
fortunately, no angular distribution measurements exist 
at the needed level of precision. The theoretical angular 
distributions of Tombrello and Parker (1963a I (see also 



Kim et al. ( 1981 1) were used to correct the prompt LUNA 
data, while the UW-Seattle data agree better with an as- 
sumed isotropic 70 angular distribution than with theory. 
As the prompt anisotropy corrections can be comparable 
to the overall quoted cross section uncertainty, we de- 
cided to exclude the prompt data from our analysis. We 
do this in part because little additional precision would 
be gained by combining the highly correlated prompt and 
activation data. Hence we base our analysis on activation 
data, plus the ERNA recoil data. 

The ERNA data and the older data of IParker and Ka^ 



vanagh ( 1963 ) extend well above 1 MeV, where measure- 
ments may provide information useful for constraining 
theoretical models of S34(i<^). Of these two data sets, 
only ERNA shows evidence fo r a significant r i se in S 34(ij^) 
above 1.5 MeV (see Fig. 1 of |Di Leva eFall p009| )). 



B. Theory 

Relative (but not absolute) S-factors at energies below 
1 MeV vary by only a few percent among credible mod- 
els, with small differences arising from non-external con- 
tributions and initial-state phase shifts. The two bound 
states of ^Be populated by '^IIe(a,7)^Be direct capture 
have large overlaps with He -I- ''He cluster configurations. 
The Pauli principle requires radial nodes in these over- 
laps, guaranteeing a small (but nonzero) short-range con- 
tribution because of cancellation in the matrix-element 
integral. 

Considerable accuracy below 1 MeV can be achieved 
by a pure external-capture model, with hard-sphere scat- 



tering at a radius chosen to reproduce measured phase 
shifts. In such a model '^He and ^He are treated as 
point particles, and final states are modeled only by 
their long-range asymptotic parts. This is the approach 
of the Tombrello and Parker (1963a) model, used to fit 



S34 in Solar Fusion I. A more realistic treatment of con- 
tributions from 2.8 to 7.0 fm is provided by potential 



models ( 


Buck et al. 


|1985 


Buck and Merchant 




1988 


Dubovichenko and Dzhazairov-Kakhramanov 1995 


Kim 


et al. 


1981 


Mohr 


2009 


Mohr et al. 


1993 


), which gen- 



erate wave functions from a Woods-Saxon or similar po- 
tential, constrained by measured phase shifts. 

Microscopic models take explicit account of nucleon 
short-range correlations. In the resonating-group method 
(RGM) a simplified nucleon-nucleon interaction is tuned 
to observables in the system being investigated (e.g., en- 
ergies of the ^Be bound states), and the phase shifts are 
computed, not fitted. The RGM wave functions are sums 
of states consisting of simple cluster substructure; in most 
^Be calculations, they are antisymmetrized products of 
Gaussians for "'He and '^He, multiplied by a function of 
the coordinate describing cluster separation. 



The RGM cal culations of |Kajino| ( |1986| and the 
potential-model of 'Langanke ( 1986^ (which employed an- 
tisymmetrized many-body wave functions) predicted the 
energy dependence of the '^H(a,7)^Li reaction quite ac- 
curately, prior to the precise measurement of |Brune et al.\ 
(1994). On the other hand, there is some variation of the 
computed ■^He(Q;, 7)^Be S-factors among RGM models 
using different interaction types and different Gaussian 
widths within the clusters. This variation has been shown 
to correlate with measures of the diffuseness of the ''Be 
ground state (Csoto and Langanke 2000:,Kajino 1986). 



Substantial changes in the S-factor and phase shifts also 
occur when ^Li-|-p configurations are added to the RGM 



and Hofmann 1986 ) 



wave functions ( Csoto and Langanke 2000 Mertelmeier 



Calculations using highly accurate nucleon-nucleon po- 
tentials are now possible. In |Nollett, (2001) , both bound 
states were computed using the variational Monte Carlo 
method, while the relative motion of the initial-state nu- 
clei was modeled by one-body wave functions from the 
earlier potential-model studies. This approach should 
provide additional realism to the nuclear wave function 
at short range, and it features initial states that fit the 
measured phase shifts. It produced very nearly the same 



S34(£') energy dependence as Kajino (1986), and an ab- 
solute S34(0) that is lower by about 25^^ 

Through a numerical coincidence, the branching ratio 
for captures to the two final states is very nearly con- 



stant at low energy (Kajino 1986). This circumstance 



and the external-capture nature of the reaction suggest 
that laboratory data can be extrapolated to low energy 
by fitting a single rescaling parameter that multiplies a 
model S34(£^) to match the data. Such a rescaling does 
not have a strong physical justification for microscopic 
models, as they do not have undetermined spectroscopic 
factors. However, rescaled microscopic models should be 
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at least as accurate as potential models and more accu- 
rate than the hard-sphere model. 

A different approach was followed in |Cyburt and 
Davids (2008), where a parameterized function fit was 
made to three of the four modern data sets over a wider 
energy interval than we used to determine our recom- 
mended 834(0) (see below), with the result 834(0) = 
0.580 ± 0.043 keV b. Their fitting function is mo- 
tivated by recent work emphasizing external capture 



and subthreshold poles in low-energy 8-factors ( 


Jennings 


et al. 


1998a|b 


Mukhamedzhanov and Nunes 


20021, and 



it matches expressions for zero pha se shift derived in 
Mukhamedzhanov and Nunes (2002). For 834, the d- 



waves have small phase shifts, and the function describes 
d-wave capture quite well. In the more-important s-wave 
capture, the function does not match detailed models of 
834 (E'), irrespectiv e of fitted parameters; its closeness to 
the expressions of Mukhamedzhanov and Nunes (2002) 



suggests that some other functional form is needed to 
account for nonzero phase shifts. 



1. Model selection for 834(0) determination 

To determine 834(0) from exp erimental capt ure data, 
we use the microscopic models of'Kajino' (19861) and Nol- 
lett ( 2001[ | (Kim A potential), rescaled to fit the data 
below E = 1 MeV (see below). We selected these two 
models based on several factors. 



i) They both accurately reproduce the s-wave phase 
shifts (as given by the phase-shift analysis of 



Tombrello and Parker (1963b)) and the long-range 
asymptotics of the 'Be bound states. The Kajino 
model reproduces the phase shifts without having 
been fitted to them. 



ii) They contain more short-range physics than hard- 
sphere or potential models, which may extend the 
energy range over which they describe the reaction 
correctly. 

iii) They agree well with each other even though they 
were generated by very different computational ap- 
proaches. 

iv) They reproduce the measured energy dependence 
of S34(i?) well, up t o at least E = 1.5 M eV (see 
Fig.[5| also Fig. 3 of |Di Leva eTaLl^m^ ). 



v) They calculate other electromagnetic observables in 
^Li and ^Be, that are in reasonable agreement with 
experiment. 



MeV, suggesting that the calculations are most reliable 
at lower energies, where poorly-constrained short-range 
contributions to 834(2?) are minimized. In Nollett ( 2001| ), 
the contribution of ■^He-^He separations less than 4 fm 
was about 4% of 834(0) and about 8% of 834(1 MeV). 
8incc a uniform 4% at all energies could be absorbed into 
the rescaling, the difference between short-range contri- 
butions at and 1 MeV suggests 4% as a conservative 
estimate of the rescaling error. 



3. Theoretical uncertainty in the 834(0) determination 

We estimate a theoretical uncertainty in the 834(0) de- 
termination by rescaling several models to the capture 
data in the same manner used to determine the rec- 
ommended value of 834(0), and examining the resulting 
spread in 834(0) values. We restrict our consideration 
to microscopic models that reproduce the s-wave phase 



shifts, cho o sing t hose of Walliser et al. ( 1984 ) , Csoto and 
Langanke ( 2000 ) (only those with ■^He -I- ''He clusteri- 
zation ), Nollett ( 2001[ ), and new variants of the Nollett 
(2001 ) calculation possessing phase shifts perturbed from 



the empirical values. 

The full spread among th e chosen set of m odels is 
±0.030 keV b, relative to the|Kajino| (|1986|) and [Nollett 



(2001) (Kim A potential) fits. We somewhat arbitrarily 
recommend two-thirds of this value; i.e., ±0.02 keV b, as 
an approximate 1-cr theoretical error. The scatter among 
these models is not independent of the rescaling uncer- 
tainty estimated above; hence, we have not included an 
explicit rescaling contribution in this estimate. 



4. S-factor derivatives 

The data do not provide a useful constraint on low- 
energy derivatives of 834(2?). Microscopic models that 
reproduce the phase shifts and simpler models that 
focus on wave-function asymptotics produce values of 
834(0)7834(0) in the range -0.55 to -0.79 MeV'^. These 
values depen d on both the model and the method of esti- 
mation. Only Williams and Koonin ( 1981 ), Walliser et al. 



(1983), and Walliser et al. (1984) published enough infor- 



mation to allow one to extract an estimate for 8"a, yield- 



ing 8^(4 (0)7834(0) = 0.26 to 0.43 MeV- 



34 J 

We base our 



recommendations on the Nollett ( 2001[ ) (Kim A) model. 



which yields effectively 8^4(0)7837(0) = -0.64 MeV'^ 
and 8^(4(0)7834(0) = 0.27 MeV'^ from a quadratic fit 
below 0.5 MeV. 



5. Comment on phase shifts 



2. Region of 534(2?) fitting 

We restricted the energy range for fitting to i? < 
1 MeV. The scatter among models (which differ mainly 
at short range) becomes much larger at energies above 1 



As the bound-state ^Be wave functions have known 
asymptotic forms, differences of the low-energy S{E) 
among models arise from differing s-wave phase shifts and 
from short-range contributions. The short-range contri- 
butions, which are difficult to compute convincingly, are 
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0.0 0.2 0.^ 0.6 0.8 1.0 1.2 

E (MeV) 



TABLE III Experimental 834(0) values and 1-a uncertainties 
determined from fits of the scaled NoUett (Kim A potential) 
theory to published data with E < 1.002 MeV. Total errors 
are quoted, including inflation factors, and systematic errors 
of LUNA: ± 2.9%; Weizmann: ± 2.2%; UW-Seattle: ± 3.0%; 
ERNA: ± 5.0%. 



Experiment 


834(0) 
(keV b) 


Error 
(keV b) 


Inflation 
Factor 


LUNA 


0.550 


0.017 


1.06 


Weizmann 


0.538 


0.015 


1.00 


UW-Seattle 


0.598 


0.019 


1.15 


ERNA 


0.582 


0.029 


1.03 


Combined result 


0.560 


0.016 


1.72 



FIG. 5 (Color online) 834 (£^) vs. E. Data points: LUNA - 
green circles; Weizmann - red squares; UW-8eattle - blue dia- 
monds; ERNA - brown triangles. 8olid curve - best fit scaled 
NoUett theory to the data with E < 1.002 MeV. The yellow 
band indicates the zbl-cr error band. Data are shown with 
statistical-plus- varying-systematic errors only; overall system- 
atic errors are not included. 



probed by capture experiments above 1 MeV. With the 
exception of |Mohr et Z] ( |1993[) and |Mohr| ( [2009] ), phase- 
shift fitting for studies of the '^He(a,7)'Be reaction has 
been based almost entirely on the phase-shift analysis of 



Tombrello and Parker (1963b I. While this phase-shift 



analysis provides a useful constraint, it depends mainly 
on a single experiment from the early 1960s, and it does 
not include an error estimation. The modern lMohr et al.\ 



(1993) experiment extended to lower energies, but it has 



no published error estimate or phase-shift analysis. 



C. 834(0) determination 

Figure [5] shows the low energy data with E < 1.23 
MeV, and the fit obtained by scaling the NoUett (Kim A 
potential) theory to best match the data with E < 1.002 
MeV. We used the analytic function 

S,,{E) = 834(0) 

X {l-OAObiE^ + 0.577 E^ ^0.1353E*), (35) 

where E is in units of MeV. Below one MeV this expres- 
sion is valid to better than 0.3%, on average. 

The best-fit curve in Fig. [5] was obtained by fitting 
each data set separately with the scaled theory, and then 
fitting the set of four 834(0) values to determine the mean 
834(0) value and its error. 

As can be seen from Table |III[ the fits to the individ- 
ual data sets are good, indicating consistency with the 
theoretical energy dependence, within the limited energy 
ranges of each set. The fit to the combined set of four 
8(0) values is of marginal quality, indicating a lack of 
good agreement in the absolute normalizations of the dif- 
ferent experiments. The combined fit has ^ jdoj = 2.3 



{dof = 3), corresponding to P(x^do/) = 0.07. All of 
the errors given in Table [ITT| include the infiation factors 
determined from the goodness of fit (see the Appendix, 
8ec. XIII.B I. Fits to these data using the scaled theory of 
Kajino yield slightly smaller values, and reproduce the 
low-energy UW-Seattle data somewhat better; however, 
the mean 834(0), 0.561 keV b, is essentially identical to 
the result obtained with NoUett's theory. 

We have focused here on measurements published since 
8olar Fusion I. We do so because in general they are 
better documented than the older ones, and address is- 
sues such as contaminant ''Be production in a quantita- 
tive manner that lends greater confidence to the results. 
One may judge from the Kajino-fit analysis presented 



Brown et al. (2007), that including older measurements 



would lower the mean 8(0) by at most 0.01 keV b or so. 
Thus including the older measurements would not change 
our result significantly. 

Given the marginal quality of the mean experimental 
834(0) fit, we round off the values given above, and quote 
a "best" result, 

834(0) = 0.56 ± 0.02(expt) ± 0.02(theor) keV b, (36) 

based on activation data and the ERNA recoil data, and 
taking the theoretical error from 8ec. |VI.B.3| 

Our best 834(0) estimate may be compared to the value 
834(0) = 0.53 ± 0.05 keV b given in Solar Fusion I. 

New capture experiments below 1 MeV would be most 
valuable for reducing the experimental uncertainty in 
834 (£^), particularly ones that maximize overlap with the 
existing modern data sets. New scattering and capture 
experiments above 1 MeV, as well as precise angular dis- 
tribution measurements, could be useful for constraining 
future theoretical calculations. ^ 



^ Note added in proof: Recent fermionic molecular dynamics 
(FDM) calculations l |Neff et al.\\M^ of S34(E) are in excellent 
agreement, in both absolute magnitude and energy dependence, 
with the experimental data shown in Fig. |5] and with the high- 
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VII. THE ^He(p,e+!/,)*He REACTION 

The hep reaction 

p + ^He '^He + e+ + I'e 



and hen reactions were performed by Carlson et aL|( 1991 ) 



(37) 



is the source of the pp chain's most energetic neutri- 
nos, with an endpoint energy of 18.8 MeV. The Super- 
Kamiokande and SNO collaborations have placed inter- 
esting limits on the hep neutrino flux by searching for 
these neutrinos in the energy window above the neu- 
trino endpoint, even though the expected flux is very low 
(see Fig. [?]). The hep rate is beyond the reach of current 
experiments: this process is induced by the weak interac- 
tion and further suppressed by a Coulomb barrier and by 
other aspects of the nuclear physics, as explained below. 
Thus theory provides our only estimate of Shop- 

The calculation of Shop is a difficult challenge. The 
leading one-body (IB) Gamow-Teller (GT) transition op- 
erator cannot connect the main s-state components of the 
p-|-^He and ^He initial- and final-state wave functions.^ 
Hence, at the IB level the reaction proceeds through the 
small components of the '^He and ^He wave functions, 
such as d-state components. Consequently, the relative 
importance of other transition operators, such as axial 
meson-exchange currents (MEC), is enhanced, as is the 
contribution from p-wave p-l-'^He capture, normally kine- 
matically suppressed at solar temperatures. The situa- 
tion is further complicated by the fact that the axial IB 
and MEC "corrections" have opposite signs, making s- 
wave hep capture even more suppressed. 



A. hep calculations 

Some of the features mentioned above are shared by 
the hen process (n-l-'^He -^'^Ke+j), in particular the 
strong suppression of IB contributions. The possibility 
of deducing Shop from the known hen cross section was 
explored in early studies: while these reactions are not 
isospin mirrors, there is a close relationship between the 
isovector spin contribution to hen and the GT contribu- 
tion to hep. However the hep S-factors determined in 
these studies differed, in some cases, by orders of magni- 
tude. 

In an attempt to understand the origin of this large un- 
certainty, fully microscopic calculations of both the hep 



energy ERNA data up to 2.5 MeV. The FDM is a nearly ah initio 
microscopic method employing realistic effective interactions. 
^ While the radial wave functions of the four nucleons in ''He can 
all be Is, with the various single-particle states distinguished 
by spin and isospin, this is not the case for the three protons 
in p-|-'^He: the Pauli principle requires that one must be radi- 
ally excited. The GT transition operator does not alter radial 
quantum numbers, only spin and isospin. Thus the GT matrix 
element between p-l-'^He and *He is suppressed due to the s-wave 
orthogonality. 



and Schiavilla et al. (1992), using a realistic Hamiltonian 
with two- and three-nucleon interactions. Among the 
approximations made in the 'Schiavilla et al. ( 1992 ) cal- 
culation were the description of the p-l-'^He initial state 
as s-wave and the omission of the dependence of the weak 
operators on the lepton pair momentum. Corrections to 
the IB GT operator were evaluated, with the largest two- 
body (2B) contributions coming from the excitation of in- 
termediate A-isobars. The A-isobar degrees of freedom 
were explicitly included in the nuclear wave functions, 
using a scaled-down approach to the full -t- A coupled- 
channel problem known as the transition-correlation op- 
erator method. Carlson et al. { 1991 1 and Schiavilla et al. 



f 1992) found that effects such as the different initial-state 
interactions for n-|-^He and p-f'^He were so substantial 
that the known hen cross section was not a useful con- 
straint on hep. Two estimates were given for the hep 
S- factor at zero energy ( ,Schiavilla et al. 19921, 



Shcp(O) = 




X 10" 



keV b. 



(38) 



depending on the method used to fix the weak N — A 
coupling constant, g^NA'- the larger of the results cor- 
responds to the naive quark model prediction for gpNA, 
while in the smaller, g^NA was determined empirically 
from tritium (3 decay. The Solar Fusion I best value for 
Shop is the average of the values in Eq. ( 38 1 . 

This problem was revisited nearly a decade later, fol- 
lowing improvements in the description of bound and 
continu um four-body wave fun ctions. The wave func- 
tions of Marcucci et al. (2000) were obtained with the 



correlated-hyperspherical-harmo nics ( CHH) variational 



method (Viviani et al. 1995 1998), using the Ar- 



gonne wis (AV18) two-nucleon ( Wiringa et al. 19951 and 



Urbana IX (UIX) three-nucleon interactions (Pudliner 



et al. 1995). The method produced binding energies of 



''He and "He and the singlet and triplet p-f ^He scattering 
lengths in excellent agreement with experiment. 



The Marcucci et al. 



(2000) calculation included all s- 
and p-wave capture channels in the p+^He initial state 
and all multipole contributions in the expansion of the 
weak vector and axial-vector transition operators. The 
weak operators corresponding to the space component 
of the IB weak vector current and the time compo- 
nent of the IB axial current, both of order v/c, have 
significant exchange-current corrections of the same or- 
der from pion-exchange. These two-body operators were 
constructed to satisfy (approximately) the constraints of 
current conservation and PCAC (partial conservation of 
the axial- vector current). Corrections to the allowed GT 
operator include both {v/cY IB and exchange-current 
contributions. The treatment of the latter followed ICar 1-1 

( p92] ) in using 



son et al. (1991) and Schiavilla et al 



the transition-correlation operator scheme and in fixing 
gi3NA to the experimental GT strength in tritium /3 de- 
cay. 
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TABLE IV Shep in units of 10"^" keV b, calculated with CHH 
wave functions genera ted from the AV18/UIX Hamiltonian 



(Marcucci et al. 20001 for three p+ He center-of-mass ener- 
gies E. The "One-body" and "Full" labels denote calculations 
with the one-body and full (one- and two-body) nuclear weak 
transition operators. Contributions from the ^Si channel and 
from all s- and p-wave channels are listed separately. 





E = 


keV 


E = 


5 keV 


E = 


10 keV 




^Si 


s-l-p 


^Si 


s-l-p 


^Si 


s-l-p 


One-body 


26.4 


29.0 


25.9 


28.7 


26.2 


29.2 


Full 


6.38 


9.64 


6.20 


9.70 


6.36 


10.1 



Table IV gives the resulting Shep at three center-of- 
mass energies. The energy dependence is rather weak. 
The p waves have a significant effect, accounting for 
about one-third of the total cross section at E=0. De- 
spite the dehcacy of the calculation, [Marcucci et al.\ 
(2000) concluded that the degree of model dependence 
was moderate: the calculations were repeated for the 



older Argonne wi4 (Wiringa et al 



and Urbana VIII (Wiringa 19911 three-nucleon inter- 
actions, but the predicti ons for Shep di f fered only by 



1984) two-nucleon 



6%. The best estimate of [Marcucci et a/.| ( |2000[ ), Shep = 
(10.1 ± 0.6) X 10~^° keV b, is about four times the value 
given in Solar Fusion I. 

A further development came with the use of heavy- 
baryon chiral perturbation theory (HBChPT) to derive 
the n eeded electroweak current operators systematically, 
with Park et al. (20031 carrying out the expansion to 
next-to-next-to-next-to leading order (N'^LO), thereby 
generating all possible operators to this order. These 
operators represent the short-ran ge physics that res ides 
above the scale of the EFT, which [Park et al.\ i2003l de- 



fined via a Gaussian regulator with a cutoff A, a param- 
eter that was varied in the calculations between 500 and 
800 MeV 



(see Sec. 



III). Shep was obtained by calculat- 
ing the matrix elements of these EFT current operators 
with phenomenological wave functions, obtained using 
the AV18/UIX Hamiltonian and the CHH method. (See 
Sec. [HT] for a more extended discussion of such hybrid 
EFT* approaches.) 

To this order, the resulting currents are IB and 2B: 
three-body operators arise at order N^LO. The expansion 
reproduces the one-pion exchange-current corrections to 
the space component of the vector current and charge 
component of the axial current, as dictated by chiral sym- 
metry, while the time component of the vector current 
has no MFC corrections. The MFC contributions to the 
axial GT operator include both a one-pion-exchange term 
and a (non-derivative) two-nucleon contact-term. The 
low-energy constant determining the strength of the con- 
tact term must be determined from an observable. Fol- 



lowing the treatment of 5/3 at a by Marcucci et al. (20001, 
this was done by fitting the GT transition strength ex- 
tracted from tritium /3 decay. 

Table [V] gives the values determined by [Park et al. 



TABLE V The hep GT matrix element Li{q;A) (in fm^^^) 
for the transition from the initial ^Si p-f -'^He state to the fi nal 
''He state, as a function of the cutoff A (Park et al. 2003), at 
E=0. Li{q;A) is evaluated at g = 19.2 MeV, the momentum 
carried out by the lepton pair. Shep (in 10"^" keV b) is also 
given. 



A (MeV) 




500 


600 


800 


Li(g;A): 


IB 


-0.081 


-0.081 


-0.081 


Li{q;Ay. 


2B (no contact term) 


0.093 


0.122 


0.166 


Li{q;A): 


2B (with contact term) 


-0.044 


-0.070 


-0.107 


Li{q;A): 


2B-total 


0.049 


0.052 


0.059 


Shop 




9.95 


9.37 


7.32 



(2003) for Shep(O) and for the GT matrix element be- 
tween the ^Si p+'^He initial and the "^He final states, as 
a function of A. By fixing the strength of the contact 
term to an observable, one hopes in such hybrid EFT* 
approaches to remove most of the calculation's cutoff de- 
pendence. Heuristically, the contact term compensates 
for high-momentum components in the phenomenologi- 
cal wave functions that would not be there had both op- 
erators and wave functions been derived rigorously from 
EFT, with a common cutoff. However, the table shows 
that significant cutoff-dependence remains in the total 
amplitude because of cancellation between the IB and 
2B contributions: the variation in S hep is ~ 15%. This 



is taken as the uncertainty in the Park et al. { 2003 1 esti- 
mate for Shep, Shep(O) = (8.6 ± 1.3) x 10" ^" keV b. The 



result is consistent with that of Marcucci et al. { 2000 ) . 

The predictio n of [Park et al. { 2003 ) was used by Bah- 
call et al.\ ( [2006[ ) and by [Pefia-Garay and Serenelli[ ( [2008^ 
in their latest determinations of the hep neutrino flux, 
(j)^(hep) = (8.22 ± 1.23) x 10^ cm'^ s-\ where the 
error reflects again the 15% uncertainty quoted above. 
The value for (l)^(hep) is in agreement with the Super- 



Kamiokand e (Fukuda et al 



et al. 



^ 2001[ ) and SNO ( |Aharmin 

2006|) upper limits at go'^ronfidence level, 40 x 10'=' 



and 23 x 10 cm s , respectively. 



B. Summary 

Given the two consistent calculations presented above, 
with the internal checks on the sensitivities to input wave 
functions and to cutoffs, and given the compatibility with 
the limits established by Super-Kamiokandc and SNO, 
we recommend 



Shep(O) = (8.6 ± 2.6) X 10^20 ]^gy 



(39) 



where the uncertainty is obtained by doubling the cutoff- 
dependence found in the Park et al. ( [2003 ) calculation. 
One anticipates that the cutoff dependence would be re- 
duced if the operator expansion were carried out beyond 
N'^LO. Thus such a program could increase confidence 
in Eq. ( 39 ) and narrow the uncertainty, even without a 



26 



fully consistent treatment of both operators and wave 
functions. 

Other ancillary calculations that could strengthen con- 
fidence in this S-factor estimate include 

• new studies of the hep reaction in which a broad 
spectrum of Hamiltonian models are explored, as 



was done by Schiavilla et al. ( 1998 ) for the pp re- 
action; 

• study of related electroweak reactions where rates 
are known, such as muon capture, as was done by 



i/^; and 



Marcucci et al.\ ( |2002| ) and |Gazit| ( |2008| ) for /i" 
•^He ^ - 



• further work to understand the relationship be- 
tween the suppressed processes hep and hen. 



VIII. ELECTRON CAPTURE BY ^Be, pp, AND CNO 
NUCLEI 

Electron capture is the source of line features in the 
solar neutrino spectrum, and represents an important 
pathway for energy production in the pp chain. Solar 
electron-capture lifetimes differ substantially from labo- 
ratory values because light nuclei are highly ionized and 
because the continuum electron density is large. 

The relative rates of ''Be electron capture and 
^Be(p,7)®B determine the ppll/pplll branching ratio and 
thus the ratio of the ^Be and ^B neutrino fluxes. The 
electron capture proceeds by the mirror transition to the 
ground state of ''Li (3/2~) and by an allowed transition 
to the first excited state (l/2~, 478 keV). By normal- 
izing the solar rate to the known terrestrial decay rate, 
the nuclear physics dependence of the solar rate can be 
eliminated. The ratio of rates depends on the relative 
electron probability densities averaged over the nucleus. 
This requires a calculation of the atomic probability den- 
sities governing the K and L terrestrial electron capture 
rates, the continuum electron probability densities at the 
nucleus for the solar rate, and corrections to the solar 
rate resulting from incomplete ionization. The solar con- 
tinuum calculation was done by Bahcall ( 1962 ) , and es- 



timates of the bound-electron contributions have been 



made by |Iben, Jr. et al.\ ([1967|, [Bahcall and Moeller 
(1969), and Bahcall ( 1994| . The solar continuum calcu 
lations have typically been done by employing the Debye- 
Hiickel approximation for plasma screening. Electrons 
within the local Debye sphere screen the nuclear poten- 
tial, thus lowering the electron density at the nucleus and 
the electron capture rate, while protons penetrating that 
radius would enhance the rate. 

Our recommended rate is based on the calculation of 
Bahcall and Moeller] ([T969 ) , with updates including the 
currently adopted 'Be half-life of 53.22 ± 0.06 days, a 
total-to-continuum capture ratio of 1.217 ± 0.002 (Bah- 



call 1994), and a terrestrial L/K capture ratio of 0.040 ± 
0.006 (Voytas et al. 2001). We use the original estimate 



of Bahcall ( 1962 ) for the terrestrial K-electron probabil- 
ity at the nucleus. The result, 

i?(^Be + e") = 5.60(1 ± 0.02) X 10"^(p/Aie) 

X T^^''^[l + 0.004(T6 - 16)] s-\ (40) 

vahd for 10 < < 16, is identical to Eq. (26) of Solar 
Fusion I. Here p is the density in units of g/cm'^, Tg is 
the temperature in units of lO^K, and /ig is the mean 
molecular weight per electron. The assigned uncertainty 
of 2% is dominated by possible corrections to the Debye- 
Hiickel approximation for charge fluctuations (reflecting 
the small number of electrons within the Debye sphere) , 
and by breakdowns in the adiabatic approximation, as 
evaluated by Johnson et al. ( [1992 ) in self-consistent ther- 
mal Hartree calculations. The small rate enhancement 
they found, 1.3%, is incorporated into and dominates the 
error in Eq. (40 1. 



Despite the lack of changes since Solar Fusion I, there 
have been developments in two areas, each concerned 
with screening corrections. First, a series of precise mea- 
surements of the terrestrial electron capture rate have 
been carried out to assess the dependence of screening 
on target chemistry, which could alter the L/K ratio (be- 
cause of L-capture sensitivity to changes in the valence 
electrons). Over the pa st decade such changes, first sug- 
gested by Segre (1947), have been explored in a series 
of half-life measurements in which ^Be was implanted in 



metals and insulators, or encapsulated in fuUerene (Das 


and Rayl 20051 1 


Limata et al.\ 20061 |Nir-El et al.\ 2007 


Norman et al.\ 


2001| Ohtsuki et al.\ 2004| Ray et al. 


1999i|2002 2006 


Wang et al. 2006). The pattern of re- 



suits is somewhat confused, with claims of variations up 
to 1.1%, but with other studies limiting effects to level s 
< (0.2-0.4)% (iLimata et 



despite use of 



al. 

lost materials 



2006 



Nir-El et al. 



20071, 



witn substantially ditterent 



electron affinities. Our tentative conclusion is that the 



uncertainty assigned in Eq. ( 40 ) is sufficient to allow for 
likely variations in terrestrial screening corrections. 

Second, questions about the adequacy of solar plasma 
screening corrections, detailed in Solar Fusion I, have not 
died out. IQuarati and Scarfone (2007 2009) reconsid- 



ered the plasma fluctuation contributions to the electron 
capture rate of ^Be, concluding that corrections of 7 
10% are required. The ansatz of Quarati and Scarfone 



( 2007) wa s prev iously considered and rejected by Bah- 
call et al. ( 2002 1 , however. The influence of protons on 
the rate of 'Be electron capture in the Sun was claimed 
to be more significant by Belyaev et al. (2007) than was 
previously thought. Davids et al. (2008), however, reject 



their argument, pointing out that only the previously in- 
vestigated electromagnetic contributions of protons play 
a role, and that the approximations under which a puta- 
tive three-body electromagnetic contribution was calcu- 
lated are invalid. 

The electron captures on p-l-p and on CNO nuclei com- 
pete with the corresponding /? decays, and thus these 
rates have been conventionally normalized to solar de- 
cay rates. As electron capture and /3 decay depend on 
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the same allowed nuclear matrix element, the ratio is in- 
dependent of the nuclear physics. The result from Solar 
Fusion I, from Bahcall and May| ( 1969| ), is 



4.5 



R^'^^ipep) = 1.102(1 ± 0.01) X IQ-^'ip/fi^) 



-1/2 



[l + 0.02(T6-16)]i?^"'=(pp), (41) 



where the superscript "Tree" indicates that the relation- 
ship omits radiative corrections, which are discussed be- 
low. The range of validity is 10 < Tg < 16. 

Radiative corrections were evaluated by |Kurylov et al. 
( 2003 ) for the two pp-chain reactions under discussion, 



p + p + e 

7t 



d + j/g 
'^Li + Ve- 



'Be-Ke" 

The radiative corrections were given as 



(42) 
(43) 



1 Capt 
pTree 
^ Capt 



l+-5Capt(i?e,Q) =C''^^(i?e,g), (44) 



where Fcapt is the total decay width, T^^^^^ is the 
tree-level width without radiative corrections, and 
9CaptiEe,Q) is a calculated factor that depends on both 
the total energy of the captured electron and the Q- 
value of the transition. Figure |6] shows the resulting cor- 
rection factors. 

Because Eq. ( 40 ) corresponds to a ratio of stellar and 
terrestrial electron capture rates, the radiative correc- 
tions should almost exactly cancel: although the initial 
atomic state in the solar plasma differs somewhat from 
that in a terrestrial experiment, the short-range effects 
that dominate the radiative corrections should be simi- 
lar for the two cases. [Indeed, this is the reason the pp 
and ^Be electron corrections shown in Fig. [g] are nearly 
identical.] However the same argument cannot be made 
for the ratio of pep electron capture to pp /3 decay, as 
the electron kinematics for these processes differ. With 
corrections Eq. (41) becomes 



i?(pep) = ^g^^^^^l.l02(l±0.01)xlO-4(p/^,) 



X n 



-1/2, 



l + 0.02(r6-16)]i?(pp). 



(45) 



where the radiative corrections have been averaged 
over reaction kinematics. Kurylov et al. (2003) found 
a 1.62% radiative correction for the /3 decay rate, 
(C'''''^(pp)) ^ 1.016 (see discussion in Sec. Ill), while 
(C'''^'^(pep)) ~1.042. Thus (C"''^'^(pep))/(C"^pp)) - 
1.026, so that our final result becomes 

i?(pep) = 1.130(1 ± 0.01) X 10"''(p/Aie) 



-1/2 



[l-h0.02(r6- 16)]i?(pp). (46) 



While certain improvements could be envisioned in the 
Kurylov et al. (2003) calculation - for example, in the 
matching onto nuclear degrees of freedom at some char- 
acteristic scale ~ GeV - rather large changes would be 
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FIG. 6 (Color online) Calculated radiative corrections for 
p-|-p-|-e~ — >■ d +Ve (dashed line) and '^Be -f — >■ ^Li + Ue 
(dotted line). The solid line is for p+e^ — )■ n + u^. Figure 
from [Kurylov et~al] (|2003[) . 



needed to impact the overall rate at the relevant 1% level. 
For this reason, and because we have no obvious basis for 
estimating the theory uncertainty, we have not included 
an additional theory uncertainty in Eq. (46). However, 



scrutiny of the presently unknown hadronic and nuclear 
effects in gcapt{Ee,Q) would be worthwhile. As one of 
the possible strategies for more tightly constraining the 
neutrino mixing angle 9i2 is a measurement of the pep 
flux, one would like to reduce theory uncertainties as 
much as possible. 

The electron capture decay branches for the CNO iso 
topes "N, 1^0, and ^^F were first estimated by |Bahcall| 



( 1990 ) . In his calculation, only capture from the con- 
tinuum was considered. More recently, |Stonehill et al. 



(2004) have re-evaluated these line spectra by including 
capture from bound states. Between 66% and 82% of the 
electron density at the nucleus is from bound states. Nev- 
ertheless, the electron-capture component is more than 
three orders of magnitude smaller than the /3+ compo- 
nent for these CNO isotopes, and it has no effect on en- 
ergy production. However, the capture lines are in a 
region of the neutrino spectrum otherwise unoccupied 
except for *B neutrinos, and they have an intensity that 
is comparable to the ^B neutrino intensity per MeV (Fig. 
[t]), which may provide a spectroscopically cleaner ap- 
proach to measuring the CNO fluxes than the continuum 
neutrinos do. 

The recommended values for the ratio of line neutrino 
flux to total neutrino flux are listed in Table I VII 

The ratio depends weakly on temperature and density, 
and thus on radius in the Sun. The values given are for 
the SSM and do not depend significantly on the details 
of the model. The branching ratio for ^Be decay to the 
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Neutrino Energy (MeV) 

FIG. 7 (Color online) Solar neutrino fluxes based on the "OP" 
calculations of ,Balicall et al. ^{^2005) , with the addition of the 
new line features from CNO reactions. Line fluxes are in 
cm~^ s~^ and spectral fluxes are in cm~^ s~^ MeV~^. Figure 
adapted from Stonehill et al. ( 2004 \ . 



TABLE VI The ratios of neutrino line intensity to the total 
intensity, after integration over the solar model. 



Source 

p+p 

^He+p 



^Be 
«B 

13n 
15q 



'total 
3 a 

8 b 

■7 c 



2.35 X 10 
4 X 10 
< 7 X 10 
0.8951 
0.1049 " 
2 X lO"'^ 
7.9 X lO""* 
4.0 X lO""* 
5.9 X lO""* 



Ref. 



Bahcall 



Bahcall 
Bahcall 



(19901 



19901 
( ,1990 1 



see text 



Bahcall 



Stonehill et al. 



Stonehill et al. 



Stonehill et al. 



1990 



2004 



2004 



2004 



"includes a 2.6% radiative correction from |Kury lev et aL| l [2003| l 
''to ^He ground state 
'^to ■'He excited state 
''to ^Li ground state 
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first excited state in the laboratory is a weighted average 



of the results from 


Balamuth et al 


( 


1983 


), Donoghue 


et al. (1983), 


Mathews et al. 


(1983 


1,1 


Davids et al. 


(1983 


1, 


Norman et al. 


(jl983a>p, and 


an average of earlier results. 



average, 10.45 ± 0.09% decay to the first excited state, 
is corrected by a factor 1.003 for the average electron 



energy in the solar plasma, 1.2 keV (Bahcall 1994), to 



IX. THE ^Be{p,7)**B REACTION 

The ^Be(p,7)^B reaction at low energies is predomi- 
nantly nonresonant El, s- and d- wave capture into the 
weakly-bound ground state of '^B (Robertson 1973). At 



solar energies the reaction proceeds by external direct 
capture, with matrix-element contributions dominated 
by ^Be-p separations on the order of tens of fermis. 
The energy dependence near the Gamow peak cannot 
be determined from simple extrapolations of higher en- 
ergy data, but must be taken from models. The narrow 
1+ resonance at Ep = 720 keV as well as resonances at 
higher energies are usually treated separately, and have 
little influence on solar rates. 

In Solar Fusion I only one direct ^Be(p,7)®B measure- 
ment was found to be sufficiently well documented to 
allow an independent assessment of the systematic er- 
rors. Consequently the recomm ended 817(0) wa s based 
on a single experiment, that of Filippone et al. (19831. 
Since Solar Fusion 1 new direct 'Be(p,7)*B measure- 
ments have been carried out at Bordeaux/Orsay (Ham- 



mache et al. 



et al. 2003a 



1998, 2001), the Weizmann Institute (Baby 



bl(see also Hass et al. (1999)), Bochum 



( [Strieder et al. 2001) and the Unive r sity o f Wash i ngton - 



Seattle/TRIUMF (jJunghans et al 



2010 2002 20031. 



These modern measurements form the basis for our Solar 
Fusion II Si7(0) recommendation. 

Other new measurements include two performed with 



2009 



Gialanella et al. 



^Be beams (Bardayan et 
2000). Although inverse measurements of this sort are 



much more difficult, they offer the attraction of different 
systematic errors. However, these experiments did not 
reach a precision useful for our purposes and thus play 
no role in our current assessment. 

In addition to direct measurements, Si7(0) has been 
determined indirectly from Coulomb dissociation, as 
summarized below in Sec. |1X.C[ and from peripheral 
heavy-ion transfer and breakup reactions. General as- 



pects of such techniques are discussed in Sec. XII 



yield a recommended branching ratio of 10.49 ± 0.09%. 



A. The direct ^Be(p,7)^B reaction 

All modern ''Be(p,7)*B experiments have employed the 
same basic method of counting /3-delayed as from the 
decay of *B to determine the reaction yield. However, 
different experimental techniques were used, and differ- 
ent levels of precision were achieved in the procedures 
for converting measured yields into cross sections and S- 
factors. Below we discuss the most important issues. 



1. Beam-target overlap 

In a conventional experiment with a beam area smaller 
than the target area, it can be difficult to determine ac- 
curately the overlap of the beam with the target, due 
to non-uniformities in the areal density of typical tar- 
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gets. This is frequently the case for radioactive target 
experiraents, as target designs are often quite compact, 
with cross sections comparable to the beam area, in or- 
der to minimize unused target material. This potential 
problem has been avoided in the most recent ^Be(p,7)®B 
experiments by using small-area targets irradiated by 
uniform beam fluxes. The reaction yield is then pro- 
portional to the product of the beam flux and the total 
number of ''Be atoms. The latter quantity can be deter- 
mined accurately from the ^Be decay radioactivity. As 
the target density may have tails extending to large radii, 
and as the beam density may not be perfectly uniform, 
it is necessary to carry out ancillary measurements to 
demonstrate the accuracy of this technique. Measure- 
ments can include separate determinations of the radial 
dependence of the beam density and the target density, 
and/or the radial dependence of the product of the beam 
and target densities. While the Bochum, Weizmann, and 
UW-Seattle/TRIUMF experiments aU used the small- 
area target/uniform-beam-flux method, only the latter 
two experiments provided sufficient information to per- 
mit an independent assessment of procedures. 



2. B backscattering 

A systematic error in ''Be(p,7)^B measurements that 
was identified after Solar Fusion I is the loss of ^B re- 
action products due to backscattering out of the target 
( [Strieder et~al\ |1998[ [Weissman et al.\ \l9m) . This loss 
is particularly significant for high-Z target backings and 



low proton energy. The Filippone et al. ( 1983 ) and Bor- 



deaux/Orsay experiments used Pt backings, for which 
the backscattering corrections are significant. In the Bor- 
deaux/Orsay experiment, calculated backscattering cor- 
rections were applied to the data, while the |Filippone| 
et al. ( 1983 1 experiment was performed prior to th e iden 
tification of ^B backscattering as a serious concern. Jung 



bans et al. ( 2003 ) estimated that the backscattering cor- 



rection for the Filippone et al. ( 1983 ) data would be be- 
tween -2% and -4% (a factor of two smaller than the es- 



timate given in Weissman et al. (19981). Here we ignore 



this correction because it is well within the overall preci- 
sion claimed in the Filippone et al. ( |1983 1 experiment and 
because it is incomplete, as effects due to target thickness 
nonuniformity (unknown) and surface composition have 
not been included. 

For the other modern experiments, ^B backscattering 
losses are not an issue: the Bochum experiment used a 
low-Z backing, while the UW-Seattle/TRIUMF experi- 
ments used an intermediate-Z backing and demonstrated 
by direct measurement that backscattering losses were 
very small. The Weizmann experiment used implanted 
targets with an intermediate-Z substrate. 



3. Proton energy loss corrections 

Low-energy data must be corrected by energy- 
averaging to account for proton energy loss in the target. 
This requires knowledge of the energy loss profile of the 
target and the target composition, as well as the monitor- 
ing of possible carbon buildup during bombardment. The 
most detailed determination of these quantities was made 
in the UW-Seattle/TRIUMF experiments, where the tar- 
get profile was determined from the narrow (F << 1 keV) 



''Be(a,7)^^C resonance at — 1377 keV. In Junghans 



et al. (2010) a more detailed resonance profile analysis of 
the previously published data was presented, allowing for 
possible depth-dependent target composition. The vary- 
ing systematic errors on the low energy "BE3" thick- 
target data were increased over the original results in 
Junghans et al. (2003) due primarily to larger assumed 



dE/dx uncertaint ies . 

In the [Filippone et al. (1983) experiment, the energy 
loss profile of the target was deduced from the measured 
shape of the 12-keV wide '^Li(p,7) resonance at Ep — 441 
keV, assuming the ^Li and ''Be distributions in the tar- 
get were the same. In the Bordeaux/Orsay experiment, 
Rutherford backscattering and (d,p) measurements were 
used to determine the target composition and proton en- 
ergy loss. In the Bochum and Weizmann experiments, 
the r = 36 keV ''Be(p,7) resonance at Ep = 720 keV was 
used to determine the proton energy loss. The Weizmann 
experiment used implanted targets with known compo- 
sition, verified by direct secondary ion mass spectrom- 
etry measurements. In the Filippone et al. ( 1983 1 and 



Bochum measurements, limits on the composition were 
inferred from the fabrication process. 

Other important factors include determination and 
monitoring of the ''Be target activity, corrections for 
sputtering losses, and determination of the efficiency for 
a detection. For the implanted target of the Weizmann 
experiment, target sputtering losses were shown to be 
negligible. The UW-Seattle/TRIUMF experiments have 
the most extensive error analysis of the modern exper- 
iments. Measurements were made with two targets of 
different thicknesses (labeled BEl and BE3) and with 
two different methods for determining the detection ef- 
ficiency for as. The resulting statistical and systematic 
errors are the smallest yet achieved. 



B. Theory 

Among the many theoretical models that have been 
published, the simplest are those in which the interac- 
tion between the ''Be nucleus and proton are described 



by a Woods-Saxon or similar potential (Aurdal 


1970 


Barker, 1980; Bertulani 1996 Davids and Typel 


2003 


Esbensenj |2004 Kim et al.\ 11987 Krauss et al. 


1993 


Nunes et al. 1997a|b| 1998; Riisager and Jensen 


1993 


Robertson 1973 Tombrello 1965| Typel et al. 


19971 



The main constraints on such moc 



els are the ground 
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state energy, the energies of low-lying resonances, and 
s-wave scattering lengths (Angulo et al. 2003). Charge 
symmetry has been used to obtain potentials from ^Li+n 
scattering lengths and the ^Li(n,7)^Li cross section, but 
persistent difficulties in simultaneously reproducing the 
absolute cross sections for ^Be(p,7)^B and ^Li(n, 7)^Li 
may reflect the greater sensitivity of neutron capture to 
the inner part of the wave function (Barker 1980[ |Es-| 
bensen 2004). Among potential models, only those of 
Nunes et al. ( f 997a|b 19981 include coupling to inelastic 
channels, open above the 430 keV threshold for excita- 
tion of ^Be. No significant effect was found, consistent 
with results of microscopic models. 

Potential models yield a reasonably accurate descrip- 
tion of the external part of the direct capture. The wave 
function at r < 5 fm is not tightly constrained in poten- 
tial models but contributes to the capture at all energies. 



particularly above 500 keV ( Csoto 1997 Jennings et al. 



1998b). However, one requirement is the existence of a 



node in s-wave scattering states, as the scattered wave 
function must be orthogonal to those o f the closed H e 
core assumed in the description of ''Be ( [Aurdal 1970). 
Model spectroscopic factors have been taken from shell- 
model studies, fixed to match transfer-reaction results 
(including the asymptotic normalization coefhcicnts dis- 
cussed in Sec. XII I, or determined by rescaling computed 



S-factors to match capture data. 



R-matrix models of direct capture (Barker 1995 



Barker and Mukhamedzhanov 2000 ) resemble potential 



models in their lack of explicit 'Be substructure, their 
need for fitting constraints, their apparent fidelity at 
large ^Be-p separation, and their relative lack of short- 
range details. Similar data are fitted and similar results 
produced. The R-matrix as applied to direct capture dif- 
fers from the discussion in Sec. [ll] only in its need for 
radiative widths and attention to the long-range tails of 
bound states (Barker 19951. 



"Microscopic" models explicitly containing eight nu- 
cleons can include substructure within ''Be and configu- 
rations not reducible to ^Be+p, calculated from the (ef- 
fective) nucleon-nucleon interaction. The antisymmetry 
between the last or scattering proton and those within 
^Be is maintained. Fully microscopic calculations to date 
generally apply versions of the resonating group method 
(RGM) to significantly simplify the many-body problem 
(|Cs6t6i[l997||Cs6t6 et a l, 199"5;' Descouvemont||2004|[De^ 
scouvemont and Baye |1988| J994 Johnson et al. 1992). 



For Si7 the interaction is usually tuned to reproduce the 
proton separation energy of *B, but may also be adjusted 
to reproduce the scattering length of ''Be-f p in the 5* = 2, 
L = channel that dominates capture at zero energy 
( Descouvemont 2004). RGM models do roughly as well 



as potential models in the external (> 5 fm) region while 
providing a more realistic description of structure in the 
internal region. Nonetheless, RGM results depend on the 
choice of nucleon-nucleon interaction and on the data 
used to fix parameters. RGM predictions of absolute 
cross sections tend to be high relative to measured val- 



ues. Thus RGM results are frequently rescaled, so that 
theory is used only to predict the energy dependence of 
S-factors, in extrapolating higher energy data to the re- 
gion of the Gamow peak. 

Other microscopic approaches have used effective inter- 
actions in combination with the shell model, adapted to 
treat weakly-bound and unbound states of p-shell nuclei 



(Bennaceur ei 1999 Halderson 2006). These studies 
focused on spectroscopic properties of A = 8 nuclei rather 
than the radiative capture. While this approach is not 
as well developed as the RGM method, it has produced 
low-energy S-factors similar to those of the RGM and 
other models. The absolute S-factor of lBennaceur et oil 
(1999) is in good agreement with the data, while that of 
pHaldersonI ( 2006 1 is ~ 40% larger than experiment. 

Ideally microscopic calculations would be carried out 
with realistic nucleon-nucleon interactions, but this is 
challenging due the complexity of the interaction and the 
need for very large spaces. The only published example 
is that of |Navratil et al. ( |2006a|b ), in which the overlap 
integrals between ®B and 'Be-fp were computed from 
seven- and eight-body wave functions obtained with the 
ah initio no-core shell model (NCSM). Due to the finite 
range of the harmonic oscillator basis, the long tails of 
the ^Be-|-p overlaps were corrected by matching their log- 
arithmic derivatives to Whittaker functions at interme- 
diate distances. These overlaps were then used as final 
states, with initial scattering states drawn from previous 
potential-model studies. The resulting Si7(0), 22.1 eV b, 
is close to the experimental value. The calculated Sn{E) 
is relatively insensitive to the choice of initial state for 
E < 100 keV, but more so at higher energies (e.g., with 
variations of 20% at 1.6 MeV). 

The envelope of predicted energy dependences of the- 
oretical models has about a 30% spread over the energy 
range fitted below. While efforts have been made to fit 
Si7(£^) with as little theoretical input as possible, some 



degree of model input appears necessary ( Cyburt et al. 



2004) 



We adopt the RGM calculation of |Descouvemont| 
(2004) as the standard to extrapolate the experimen- 
tal data to energies of astrophysical interest. Among 
available RGM calculations, this one is the most com- 
plete numerically. Of the two NN interactions used 
in Descouvemont (2004), the Minnesota interaction was 



judged to describe light nuclei more accurately. The pre- 
dicted Si7(0) = 24.69 eV b is 19% larger than our rec- 
ommended value, while the calculated shape of Si7(-E) 
provides a marginally better fit to the data, compared to 
other models we considered. Other ^B and *Li properties 
computed in this model also match experiment reason- 
ably well. Nevertheless, the substantial theoretical error 
bar assigned to our end result of Sec. IX. D| - to remove 
much of the dependence on choice of model - dominates 
the overall uncertainty in our value for Si7(0). 

Low-order polynomial representations of Si7(-B) that 
span both the solar Gamow peak and energies where data 
are available have poor convergence due to a pole in the 
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S-factor at —138 k eV ( jJennings et al. 1998a|b Williams 



and Koonin 1981 ). Thus instead we fit the models over 
a more limited energy range important to stellar fusion, 
to 50 keV. A quadratic expansion then provides a good 
representation. This procedure yields Si7(0)/Si7(0) be- 
tween — 1.4/MeV and — 1.83/MeV for the models used in 
our fitting. We recommend as a best value and probable 
range 



S'l7(0) 
Sl7(0) 



(-1.5±0.1)/MeV. 



(47) 



The corresponding values for S"7(0)/Si7(0) vary from 
7.2/MeV^ to 20.4/MeV^; we recommend 



Sl7(0) 



(ll±4)/MeV^ 



(48) 



The ranges are consistent with other published values 
where derivatives were defined by similar procedures 



(Barker 1983 Bennaceur ei aL 1999 Descouvemont and 



Baye 1988 Kolbe et al. 1988 1 . Published values outside 



2000 



1985 



our recommended ranges (Adelberger et al. 1998 Baye 
|Baye and Brainis |2000 Baye and Descouvemont 



Baye et al, 19W{ Jennings et at] 1998bj John- 



son et al. 1992, jWilliams and Koonin , 1981p are either 
mathematical derivatives at = or fits over a wider 



3 



energy interval. For the adopted Descouvemont (2004) 
model with MN potential, the corresponding numbers 
are Si7(0)/Si7(0) = -1.51/MeV and S'/7(0)/Si7(0) = 
13.5/MeV2. 



C. B Coulomb dissociation measurements 

Estimates of direct (p,7) capture cross sections can be 
derived from Coulomb Dissociation (CD) measurements 
(see Sec. XII). Because of the complexity of the associ- 
ated analysis and the absence of convincing benchmarks 
for the CD method, the Solar Fusion I authors concluded 
that it would be premature to use information from the 
CD of ^B in deriving a recommended value for 817(0). 
However, the CD of ®B was identified as a prime test 
case for this method, because this reaction can be stud- 
ied both directly and indirectly, is characterized by a low 
proton binding energy, and is dominated by El tran- 
sitions. Three groups have performed CD experiments 
with radioactive ®B beams of incident energies between 
44 and 254 A MeV. A comparison of their results to those 
from radiative proton capture allows one to assess the 
precision that might be possible with the CD method. 

Exclusive CD measurements were performed at 47 A 



MeV (Iwasa et al. 1996 iMotobayashi et al. 1994) and 



in complete kinematics including 7-rays, but had to cope 
with a large background induced by reactions in the He 
bag between the target and the fragment detectors. The 
MSU experiment suffered from a low detection efficiency, 
particularly at high p-^Be relative energies. The GSI 
experiment eliminated background by reconstruction of 
the break-up vertex and utilized a focusing spectrome- 
ter with large momentum acceptance that provided high 
geometric detection efficiency. These considerations sug- 
gest that the GSI measurement of Schiimann et al. ( 2006 ) 
represents the most complete experimental study of "B 
CD to date. 

The extraction of Si7(i?) from the differential CD cross 
section da/dE, which varies rapidly with energy, is not 
trivial. The poor energy resolution in CD experiments, 
together with the influence of experimental cuts, require 
careful simulations of this distribution using a theoretical 
model. In addition to the dominant single El photon ex- 
change, other potentially important factors are E2 transi- 
tions, nuclear break-up, and higher-order corrections. All 
of these effects are expected to be smaller at the higher 
energy of the GSI experiment than at the lower energies 
of the RIKEN and MSU experiments. However, a proper 
analysis of the GSI experiment requires relativistic mod- 
eling, a step so far taken only in perturba tion theory 
(|Bertulani[ |2005[ [Ogata and Bertulani[ |2009[). 



For the RIKEN case, 



Kikuchi et al. (19971 presented 



differential cross sections da/d9s, where 9s is the scatter- 
ing angle of the excited *B* system reconstructed from 
the ^Be and p momentum vectors, relative to that of the 
incoming ^B. The measured distribution was compared 
to first-order perturbative calculations that included El 
and both nuclear and Coulomb £ = 2 transition ampli- 
tudes. At low relative energies, the authors found good 
agreement of their measured distributions with those 
from a model that assumes only a dipole contribution. 
Later investigations of the same data employed more so- 
phisticated reaction models, stressing the importance of 



all the effects mentioned above ( 


Alt et al. 


2003 


2005 Es- 


bensen et al. 2005 Goldstein et al. 2007 Ogata et al. 


2006 


Summers and Nunes 




2005 


). For example, the 



value of Si7(0) obtained from the continuum-discretized 



coupled-channels (CDCC) analysis of Ogata et al. (2006) 
is 13% larger than that determined in the first-order cal- 
culation of Kikuchi et al. ( 1998 ) . 



At MSU, inclusive measurements were performed to 
test the prediction that interference between El and E2 
transitions in the CD of ^B would produce asymmetries 
in the longitudinal momentum distributions of the emit- 
ted fragments (Esbenscn and Bertsch 19961. Longitudi- 
nal momentum distributions of the 'Be fragments from 
the break-up of *B on Pb and Ag targets at beam ener- 



52 A MeV ( [Kikuchi et aL[|1997 199"8t at RIKEN, at 83 gies of 44 an d 81 A MeV were measured ( [Davids et af. 



A MeV at M SU Pavids et aL[ [2001a|b[), and at 254 A 



MeV at GSI (Iwasa et al. 1999 Schiimann et al. 2003 



2006[ ). For the RIKEN and GSI experiments, the most 
recent publications supersede the previously published 
ones. The RIKEN experiment measured the CD of ^B 



1998 2001b I. Asymmetries in these distributions were 
incontrovertibly observed and were interpreted with both 
first-order perturbative and CDCC calculations. The E2 
strengths deduced from first order perturbation theory 
were found to be somewhat smaller than or consistent 



32 



> 



40 - 



30 - 



20 



10 



T r 



• Schuemann et al 
i< Kikuchi et al. 
a Davids et al. 



T 




500 



1000 1500 



2000 



2500 3000 
E (keV) 



FIG. 8 (Color online) S17 values from CD experiments. 
Full red circles: latest analysis of the GSI CD experiment 
(Schiimann et al. 2006); open blue stars: Kikuchi et al. 
{ 1998 1 a nalyzed in first-order per turbation theory; open blue 
squares: Davids and Typel ( 2003 1. The error bars include sta- 



tistical and estimated systematic errors. The curve is taken 
from the cluster-model theory of |Descouvemont et al.\^200A\ , 
normalized to 817(0) = 20.8 eV b. 



periments as well as the number and precision of the 
CD Si7(i?) data points prevent a precise determination 
of the shape, which therefore has to be taken from the 
radiative-capture measurements. In Fig. [8] we display 
the best-fit curve for the direct (p, 7) data, including the 
dominant El multipolc but not the All contribution (see 
SeclDTDl). 



It is difficult to quantitatively assess the impact of the 
different theories and energy ranges used in analyzing the 
three CD experiments on the derived 817(0) values. The 
resulting values are 21.4±2.0 eV b for the RIKEN exper- 
iment, as reanalyzed by Ogata et al. ( 2006 ); 20.6±1.4 eV 
b for the GSI experiment; and IT.SriyiV b for the MSU 
experiment. Empirically these values are consistent with 
the range Solar Fusion I defined for direct measurements, 
Si7(0) = I9I2 b. Moreover, the good agreement be- 
tween the shapes of the GSI CD and the radiative capture 
data eliminates the concern about systematically differ- 
ent slopes of Si7(i?) derived from the respective methods. 
However, we believe it would be premature to include the 
CD results in our determination of a recommended value 
for Si7(0), as a better understanding of the role of E2 
transitions and higher order effects in breakup at var- 
ious energies is needed. Further discussions can be found 
in Sec. EH 



with all published models of structure. Later, the lon- 
gitudinal momentum distributions of the emitted protons 
were s tudied in the exclusive MSU measurement at 83 A 
MeV ( Davids et al.\ 2001a|b ) and found to be consistent 
with the 'Be distributions observed in the inclusive mea- 
surement. The Sn{E) distribution was extracted from 
dcr/d£' (Davids and Typel 2003 ) with a requirement that 



6'g < 1.8°, corresponding classically to an impact param- 
eter of 30 fm; a small E2 contribution derived from the 
inclusive measurements was taken into account. 



Schiimann et al. (2006) published the most extensive 



set of differential cross sections for the GSI experiment. 
All distributions were gated by 9g, < 1°, corresponding 
to an impact parameter of 18.5 fm. The measured distri- 
butions were compared to theoretical ones filtered by the 
experimental efficiency and resolution using a GEANT-3 
simulation. The event generator employed a simple first- 
order perturbation-theory description of Coulomb break- 
up with only El transitions included. The authors chose 
this simple model for its ease in numerical calculations 
and for its fidelity in reproducing, e.g., the inclusive 6'g 



distribution (Fig. 11 of Schiimann et al. (2006)) and the 
surprisingly symmetric distribution s of the protons 
in the *B* reference system (Fig. 13 in|Schiimann et al. 



(2006)). Consequently, Sn{E) was deduced from this 
model under the assumption that, contrary to theoreti- 
cal expectations, E2 transitions could be ignored. The 
data points resulting from all three CD experiments are 
shown in Fig. [sj (Note that the RIKEN data points were 
taken from the first-order perturbation-theory analysis 
by [Kikuchi et aL| ( |1998[ ).) 

The different assumptions made in analyzing the ex- 



D. Direct '^Be(p,7)®B analysis and Si7(0) determination 

Figure [o] shows the modern ^Be(p,7)*B data with 
center-of-mass energy E < 1250 keV. We analyzed the 
Filippone et al. (19831 data using the ''Li(d,p) cross sec- 



tion given in Solar Fusion I. Total errors, including sys- 
tematic errors, are shown on each data point, to facilitate 
a meaningful comparison of different data sets. All data 
sets exhibit a similar Si7(£^) energy dependence, indicat- 
ing that they differ mainly in absolute normalization. 

Following the discussion in Sec. IX. B[ we determine our 
best estimate of 817 (0) by extrapolating t he data using 
the scaled theory of [Descouvemont ( 2004| ) (MN calcula- 
tion). We performed two sets of fits, one to data below 
the resonance, with E < 475 keV, where we felt the reso- 
nance contribution could be neglected. In this region, all 
the individual 817(0) error bars overlap, except for the 
Bochum result, which lies low. 

We also made a fit to data with E < 1250 keV, where 
the 1"*" resonance tail contributions had to be subtracted. 
We did this using the resonance parameters of |Junghans| 
' d]^003) (£:p=720 keV, Fp =35.7 keV and = 25.3 
meV) , adding in quadrature to data errors an error of 
20% of the resonance subtraction. In order to minimize 
the error induced by variations in energy-averaging be- 
tween experiments, we excluded data close to the reso- 
nance, from 490 to 805 keV, where the S-factor is strongly 
varying and the induced error is larger than 1.0 eV b. 
Ab ove the resonance , the d ata have smaller errors. Only 
the Filippone et al. ( 1983 ) and Weizmann group error 



bars overlap the UW-Seattle/TRIUMF error bars. 
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FIG. 9 (Color online) Si7(i5) vs. center-of-mass energy E, for 
E < 1250 keV. Data points are shown with total errors, in- 
cluding systematic errors. Dashed line: scaled Descouvemont 
(2004) curve with 817(0) = 20.8 eV b; solid line; including a 
fitted l"*" resonance shape. 



Figure [9] shows the best-fit Descouvemont (2004) (MN 
interaction) curve from the E < 475 keV fit (together 
with the 1"*" resonance shape determined in [Junghans 
et al. (2003), shown here for display purposes). Our fit 
results are shown in Table [VTTl The errors quoted include 
the inflation factors, calculated as described in the Errors 
Appendix. The main effect of including the inflation fac- 
tors is to increase the error on the combined result by the 
factor 1.7 for £: < 475 keV, and by 2.0 for £; < 1250 keV. 
Both the 817(0) central values and uncertainties from the 
combined fits for these two energy ranges agree well, the 
latter because the added statistical precision in the E < 
1250 keV fit is mostly offset by the larger inflation factor. 

We also did fits in which the low energy cutoff was 
varied from 375 to 475 keV and the high energy exclu- 
sion region was varied from 425-530 to 805-850 keV. The 
central value of 817(0) changed by at most 0.1 eV b. On 
this basis we assigned an additional systematic error of 
± 0.1 eV b to the results for each fit region. 

To estimate the theoretical uncertainty arising from 
our choice of the nuclear model, we also performed fits 



using the sh apes from other plausible models: [Descouve 
(12004') plus and minus the theoretical uncertainty 



mont 



shown in Fig. 8 of that paper; |Descouvemont and Baye| 
(|1994D; the CD-Bonn 2 000 calculation shown in Fig. 15 of 
Navratil et al. ( 2006b ) ; and four potential model calcula- 
tions fixed alternately to reproduce the ^Li -I- n scattering 
lengths, the best-fit ''Be -I- p scattering lengths, an d their 
upper and lower hmits (Davids and Typel 2003). The 



combined-fit resul ts for all these curv es, in cluding |De 
ouvemont ( 2004 ) , are shown in Table VIII 
We estimate the theoretical uncertainty on 817(0) from 
the spread of results in Table |VIII| ± 1.4 eV b for the 
E < 475 keV fits, and eV b from the E < 1250 keV 
fits (the smaller error estimate in the latter case reflects 



TABLE VII Experimental 817(0) values and (inflated) un- 
certainties in eV b , and /dof determined by fitting the |De-| 
{|2004|) MN calculation to data with E < 475 keV 



scouvemont 



and with E < 1250 keV, omitting data near the resonance in 
the latter case. 



Fit Range 


E < 475 keV 


E < 1250 keV 


Experiment 


Sl7(0) 


a 


x'/dof 


Sl7(0) 


a 


xVdof 


Baby 


20.2 


1.4" 


0.5/2 


20.6 


0.5" 


5.2/7 


Filippone 


19.4 


2.4 


4.7/6 


18.0 


2.2 


15.8/10 


Hammache 


19.3 


1.1 


4.8/6 


18.2 


1.0 


12.5/12 


Hass 








18.9 


1.0 


0/0 


Junghans BE3 


21.6 


0.5 


7.4/12 


21.5 


0.5 


12.3/17 


8trieder 


17.2 


1.7 


3.5/2 


17.1 


1.5 


5.1/6 


Mean 


20.8 


0.7 


9.1/4 


20.3 


0.7 


18.1/5 



"We include an additional 5% target damage error on the lowest 
3 points, consistent with the total error given in the text of |Baby] 
|et al.\ l [2003a[ l (M. Hass, private communication, 2009). 



TABLE VIII Experimental 817(0) values and (inflated) un- 
certainties in eV b, and determined by fitting nine calcula- 
tions to the data sets of Table |VII| The E < 475 keV fits have 
dof = 4 and the E < 1250 keV fits have dof=5. D04 is|D&l 



scouvemont 



(20041, DB94 is Descouvemont and Baye (19941, 



and NBC06 is Navratil et al. ( 2006b \ 



Fit Range 



E < 475 keV 



E < 1250 keV 



Model 


817(0) 


a 


x' 


Sl7(0) 


a 


x' 


D04 (central) 


20.8 


0.7 


9.1 


20.3 


0.7 


18.1 


D04 (upper) 


20.1 


0.7 


10.0 


19.7 


0.7 


18.5 


D04 (lower) 


21.5 


0.7 


8.1 


21.0 


0.7 


17.3 


DB94 


21.4 


0.7 


8.4 


21.5 


0.7 


16.7 


NBC06 


22.1 


0.7 


7.4 


21.8 


0.8 


18.5 


^Be-|-p (central) 


21.2 


0.7 


8.7 


20.2 


0.7 


19.7 


^Be-|-p (upper) 


19.4 


0.8 


11.7 


17.3 


0.7 


21.6 


^Be-|-p (lower) 


21.7 


0.7 


8.2 


21.0 


0.7 


19.4 


■^Li-l-n 


20.5 


0.7 


9.7 


19.1 


0.7 


20.9 



the exclusion of the poorer potential- model fits). We 
note that the estimated uncertainties are substantially 
larger than those given in Junghans et al. ( |2003 ) and in 
Descouvemont (2004 1. 

We expect the model dependence^ 



of the fit to be 
greater above the resonance because of the demon- 
strated dependence of the S-factor in this range on 



Recently |Yamaguchi et al.\ |2009[ | discussed a contribution of 
a possible higher energy (3.2 MeV) 2~ resonance to ''Be(p,7). 
They estimate its contribution by taking the transition strength 
to be a Weisskopf unit. As low-lying El transitions are typically 
strongly inhibited, this estimate is unlikely to be realistic. Our 
S-factor estimate is based on a fit to low-energy data that would 
be free from any significant influence of this distant resonance, 
regardless of such assumptions. 
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the less-constrained short-range part of the wave func- 



17.979 2* 



tions (|Cs6t6[ |1997[ |Descouvemont[ |2004[ [Jennings et al. 
1998b|. We Ijase our 817(0) recommendation on the E < 



475 keV fit, 

Si7(0) = 20.8 ± 0.7(expt) ± 1.4(theor) eV b. 



(49) 



This value is in agreement with, but substantially more 
precise than, the Solar Fusion I recommendation, 817(0) 

-4 

-2 



19tt eV b. 



X. THE SPECTRUM OF *B NEUTRINOS 

The ®B neutrino spectrum differs from an allowed 
shape primarily because the principal state populated in 
the decay is a broad resonance. A precise determination 
of the neutrino spectrum is important to the analyses of 
the *B neutrino data obtained by the Super-Kamiokande 
and 8N0 collaborations. Uncertainties in the spectrum 
are a source of systematic error in these experiments, po- 
tentially affecting conclusions about the hep flux, MSW 
spectral distortions, etc. The neutrino spectrum can be 
determined from laboratory measurements of /S"*" de- 
cay in which the decays of final-state ^Be resonances are 
observed. 

The ^B /3+ decay from the = 2+ ground state is 
followed by the emission of t wo a particles from excited 
2+ states of ^Be (see Fig. 10 ). Although the region of in- 
terest is dominated by a single state in ®Be with ^ 3 
MeV, the width of this resonance is quite large, F ~ 1.5 
MeV. Consequently the a spectrum yields a continuum, 
so that other 2+ states need to be considered. The a 



spectrum was first measured by|Farmer and Class ( 1960 1 , 



and later by Wilkinson and Alburger (1971). R-matrix 
analyses were presented by Barker ( 1989[) and Warburton 
( 1986D (but s ee the caveat of Bhattacharya and Adel ' 
berger| ( [2002| )). [Bahcall et Z] ( |1996p used the existing 
data to produce a recommended neutrino spectrum that 
was widely used in subsequent analyses of neutrino ex- 
periments. 



Ortiz et al. { 2000 ) claimed a discrepancy with previous 
determinations of th e a spectrum. Subsequent ly [Win 



ter et al. (2003) and Bhattacharya et al. ( [2006 ) studied 
the spectrum via experiments with very different system- 
atic uncertainties, finding excellent agreement with each 



other but disagreement with the claim of Ortiz et al. 
(|2000|) . It was reported (A. Garcia, private communica- 



tion, 2009)) that Ortiz et al.\ (2000) now recognize that 



they underestimated uncertainties related to the energy 
loss generated by carbon buildup in their targets, so that 
a claim of a disagreement with earlier measurements no 
longer should be made. We recommend using the a spec- 
trum of Winter et al. (2006) and the consistent and higher 



precision spectrum of jBhattacharya et al. (2006). These 
experiments do not suffer from the energy calibration 
problems that affected earlier experiments, as discussed 




^He+'He 



FIG. 10 Energy levels from the *B(/3+)*Be(2Q) decay chain. 



trino spectrum tabulated in Winter et al. (2006) . [The 



neutrino spectrum was not calculated by |Bhattacharya| 
et aL|(|2006D.l 



The positron spectrum can be deduced from the a 
spectrum in a similar fashion, and is useful as a test of 
data consistency. The meas urements of|Napolitano et aL\ 
{ 1987 ) have been shown by Winter et al. { 2006 ) to be in 



good agreement with the results from the a spectrum. 

Forbidden corrections are at the level of a few per- 
cent. Many measurements have been performed to de- 
termine needed matrix elements (Bowles and Garvey 



vey 



19781 De Braeckeleer et al. 1995 


McKeown et al.\ 1980 


Nathan et al. 1975| Paul et al.\ 


19771 Tribble and Gar- 



1974 1975). Radiative corrections are smaller at a 



frac tion o f one pe rcent and have been calculated by |Sir-| 
lin (1967) and by Batkin and Sundaresan (1995). Both 



sets of corrections are described by Winter et al. 



(20061, 



and incorporated into the spectrum given there. Bah- 



call] (1991) showed that red-shift distortions associated 
with the Sun's gravitational potential are insignificant, 
affecting the spectrum at the fractional level of ~ 10~^. 
Bacrania et al. ( 2007 ) have placed a 90% confidence-level 
bound on the branching ratio for ®B /3 decay to the 0+ 
ground state of ^Be (a second-forbidden transition) of 7.3 
X 10~^ (see Fig. 10), limiting uncertainties in the high 
energy portion of the ^B neutrino spectrum. 



by Bahcall et al. ( 1996 ) . Finally we recommend the neu- 



* The strength function and the neutrino and positron spectra are 
in electronic repositories available onUne through Phys. Rev. C. 
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XI. THE CNO CYCLES 

The need for two mechanisms to account for the stellar 
burning of hydrogen to helium was recognized in the pi- 
oneering work of Bethe and collaborators. The pp chain, 
which dominates energy production in low-mass main- 
sequence stars, can operate in metal-free stars, synthesiz- 
ing ^He from H, while creating equilibrium abundances of 
deuterium, '^He, and ^Be/''Li, the elements participating 
in intermediate steps of Fig. [2j 

Heavier main-sequence stars produce their energy 
dominantly through the CNO cycles, where reactions are 
characterized by larger Coulomb barriers. Hence, the en- 
ergy production rises more steeply with increasing tem- 
perature (ecNO T^^ compared to epp cx T"* at solar 
core temperature, as illustrated in Fig. U. The CNO cy- 
cle was proposed by Bethe and Weizsacker to account for 
the evolutionary tracks of massive stars. Unlike the pp- 
chain, the CNO bi-cycle of Fig. [2] requires pre-existing 
metals to process H into "'He. Thus the contribution to 
energy generation is directly proportional to the solar- 
core number abundance of the primordial metals. The 
CN-cycle, denoted by I in Fig. [2] is an important SSM 
neutrino source. It also accounts for about 1% of solar 
energy generation. The cycle conserves the number abun- 
dance, but alters the distribution of solar metals as it 
burns into equilibrium, eventually achieving equilibrium 
abundances proportional to the inverse of the respective 
rates. In the Sun this leads to the conversion of almost 
all of the core's primordial '^^C into '^^N. This change in 
the chemical composition alters the core's opacity and, 
at the 3% level, the heavy element mass fraction Z, SSM 
effects first explored by Bahcall and Ulrich (19881. 

The '^'*N(p,7) reaction - the slowest reaction in the CN 
cycle at low temperatures and thus the rate-controlling 
step - determines whether equilibrium has been achieved. 
The "'''N lifetime is shorter than the age of the Sun for 
temperatures > 1.33 x 10^ K. Therefore equilibrium for 
the CN cycle has been reached only for R <0.1Rq, cor- 
responding to the central 7% of the Sun by mass. Con- 
sequently, over a significant portion of the outer core, 
^^C has been converted to ^^N, but further reactions are 
inhibited by the '^''N(p,7) bottleneck. 



A. The reaction "N(p,7)^50 
1. Current status and results 



-CM 



[keV] 



Figure 11 shows the level structure of '^^O, relative to 
the threshold energy for ^'*N(p,7). 

Solar Fusion I gave 3.5l°'g keV b as the recommended 
total S-factor for the "'''N(p,7)^^0 reaction. This was 
based on the energy dependence dete r mined by [Schroder] 
et al.\ ( |1987[ ). In the [Schroder et al.\ ( [1987[ ) analysis the 
ground state transition accounted for half of the total 
S-factor at zero energy, primarily because of the con- 
tribution of a subthreshold resonance at E= —504 keV 



2187 



987 



7297 



259 



14N + P 



-506 



-1121 



-2117 



Ex [keV] 
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FIG. 11 The energy levels of '^^O and their relationship to 
the threshold energy for ^*N(p,7). 



(corresponding to the 6.79 MeV state in '^^O). How- 
ever, a reanalysis based on an R-matrix calculation by 
Angulo and Descouvemont (20011 indicated that the 
streng th of the ground state transition in [Schroder et al. 
([1987|, Sfi4(0)=1.55 keV b, had been significantly over- 



estimated, and should be reduced to 0.08 keV b. 

This finding prompted a series of new experiments 



using 


direct (Bemmerer et al. 


2006b Formicola et al. 


2004 


Imbriani et al.[ 2005| Lemut et al.[ |2006| |Marta| 


et al. 


2008 Runkle et al. 2005 


) and indirect approaches 


(Bertone et al. 2001 


2002 1 


Vlukhamedzhanov et al. 


2003 INelson et aL||2003 


ISchiirmann et a/.l|2008||Yamada| 


et al. 


2004 1 . The prompt-capture 7-radiation was mea- 



sured in experiments by the TUNL group (Runkle et al. 



2005 


) in a surface laboratory and by the LUNA group 


dFormicola et al. 


2004 [Imbriani et al. 


2005 


Marta et al. 


2008) in Gran Sasso. From these experiments - carried 



out with Ge detectors - the contributions of each transi- 
tion could be extracted. In an additional measurement by 



the LUNA Collaboration ( Bemmerer et al. 2006b Lemut 



et al. 2006) the total cross section was determined. 
These recent experiments cover an energy range from 70 
to 480 keV, still far from the solar Gamow window at Eq 
= 27 keV. Additional information is provided by experi- 
ments that probe the width of the subthreshold state at 
E — —506 keV by the Doppler shift attenuation method 
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(Bertone et al. 2001 Schiirmann et al. 20081 and by 



Coulomb excitation (Yamada et al. 2004). Asymptotic 
normalization coefficients (ANC) for the ground state 
and selected excited states were determined from trans- 
fer reaction measurements for ^''N('^He,d)^'^0 by ,Bertone| 
et al.\i2002\ and iMukhamedzhanov et aL|(|2003|). All ex- 



periments and subsequent analyses confirmed that the 
value for the ground-state contribution determined in 
the extrapolations of Schroder et al. { 1987) had been too 
hi gh. Current estimates of S^°|4(0) range from 0.0 8 keV 
b dAngulo and Descouvemont[ 200 1[ ) to 0.45 keV b ( |Run 



kle et a/.||2005p . Hence, the S-factor for "N(p,7)i^OTs 
now determined largely by the transition to the 6.79 MeV 
state. Minor contributions arise from transitions to the 
5.18, 5.24, 6.17, 6.86 and 7.28 MeV states in ^^O. 



2. R-matrix analysis and normalization 

We have performed an R-matrix fit to the three 
strong e st transitions using t he data of |Imbriani et al.\ 
(|2005D, [Marta et al.\ (|2008|), iRunkle et al.\ (|2005|), and 



Schroder et al. ( 1987) and the code of Descouvemont (De 



scouvemont and Baye 2010). In this way we obtain the 



most robust weighted mean. The recent direct experi- 



ments ( 


Bemmerer et al. 


2006b Formicola et al. 2004 


Imbriani et al. 2005; 'Lemut et al.[ |2006| |Marta et al. 


2008 


Runkle et al. 


2005 


) cover only a relatively narrow 



energy window. Thus, as no new information is avail- 
able for the higher lying resonances, a reliable extrap- 
olation to zero energy requires the high-energy data of 
Schroder et al. (1987). However, systematic differences 



are apparent in the data sets of Imbriani et al. { 2005 1 , 



Runkle et al. (2005), and Schroder et al. (1987). In order 



to minimize systematic uncertainties, all data sets were 
renormalized to the weighted mean of the strength of the 
259 keV resonance in "N(p,7)i50. Table 



IX 



summarizes 



the resonance 



the available absolute determinations of 
strength with a weighted mean of W7259 — 13.1 ± 0.6 
meV. The uncertainty was obtained by calculating the 
error on the weighted mean, excluding the common sys- 
tematic uncertainty on the stopping power of protons in 
nitrogen ( Ziegler et al. 2008 ) . The latter was summed in 



quadrature with the weighted mean error to obtain the 
final uncertainty. 



In Schroder et al. (1987) the data were normalized to 
an absolute cross section determination at E = 760 keV, 
a{E = 760 keV) = 620 ± 80 nb. This value is an adopted 
mean based on several experimental methods, while the 
measurement relative to W7259 gives a{E = 760 keV) — 
609 nb (iSchroder et a/.l|1987|). Thus, based on the differ- 



ence between the value for ^7259 used by [Schroder et al.\ 



( |1987[ ), 14 meV (Becker et al.\ ,1982 ), and the new de- 
termination, 13.1 ± 0.6 meV, a precise renormalization 
of (j{E = 760 keV) can be made, relative to this reso- 
nance. One finds a{E = 760 kcV) = 570 nb. Moreover, 
we no te that the energy dependence of [Schroder et al\ 
(1987) was corrected for summing contributions, as dis- 



TABLE IX Summary of the published values for (^7259, along 
with their estimated statistical, systematic, and total uncer- 
tainties. All quantities are in units of meV. The last row gives 
the recommended value. 



Becker et al. (19821° 



Runkle et al. 



0J'y259 Stat. syst. total 



(20051 



14 
13.5 
12.9 

Bemmerer et al. (|2006bl) 12.8 



Imbriani et al. 



(20051 



0.4 
0.3 



1.0 

1.2 1.2 
0.8 0.9 
0.5 0.6 



recommended value 



13.1 



0.6 



"used in [Schroder et aL| ( [l987[ 



cussed by Imbriani et al. (2005). The renormalizations 



for [Runkle et al. (2005) and Imbriani et al. (2005) are 3% 
and 2%, respectively. 

The ANCs for the ground, 6.79 MeV, and 6.17 MeV 
states as well as of the 6.79 MeV state are important 
parameters in the R-matrix analysis determining S(0). 
Parameter values determined in the analysis will refiect 
the quality of the input data. Thus the R-matrix results 
can be validated by comparing these values with those 
determined independently by transfer reactions and other 
indirect measurements (see Table [X|) . 



3. Transition to the ground state and 6.79 MeV in ^^0 

The transitions to the ground and 6.79 MeV states in 
^^O are connected through the reduced proton width of 
the —0.506 MeV subthreshold state. This width can also 
be expressed in terms of the subthreshold state ANC via 
the Whittaker function at the R-matrix radius a that 



appears in Eq. (3.60) of [Descouvemont and Baye (2010) 
(see references therein). Both transitions are discussed 
together here. 

Transition to the 6. 19 MeV state: The reaction mech- 
anism for the transition to the 6.79 MeV state appears 
rather simple, primarily an external capture process 
whose magnitude is determined by the value of the ANC. 
Hence Sf jf (0) is dominated by the external capture pro- 
cess. In_the_^resentaM^ the data of [Runkle et al.\ 



(2005), Imbriani et al. (2005), and Schroder et al. (1987) 



are included after renormalization, as described above. 
As the recent low-energy data do not strongly constrain 
the R-matrix radius, high-energy data are needed. The 
resulting SiJ2{E) fails to reproduce the high-energy data 
for radii 5.5 fm < a < 6.5 fm, as in Fig. 4 of [Angulo and[ 
Descouvemontj (2001). A better fit can be obtained by 
choosing smaller radii. However, this choice also impacts 
fits for the ground state transition, which favor larger 
radii. Consequently, we have not used the transition to 
the 6.79 MeV state to determine the R-matrix radius in 
this way. Instead, R-matrix fits were done 



i) taking all renormalized data (Imbriani et al. 
IRunkle et al. 



2005 Schroder et al. 1987) into ac 



2005 



37 



TABLE X Published ANC values and T-y for the 6.79 MeV transition. All ANC values are given in the coupling scheme of 
Angulo and Descouvemont (20011. The recommended values in the last row were obtained as a weighted mean considering 



as weights the experimental errors only. Finally, the recommended uncertainty was obtained by summing in quadrature the 
weighted mean error and an average theoretical uncertainty. The latter is according to information provided by the authors. 
As existing measurements of r.^(6.79 MeV) are discrepant, no recommended value is given. 



C„ 



,3/2 



(fn 



-1/2 



.^6.79 



(fm-1/2) C, 



6.171/2 



(fm- 



-1/2-1 i) 



C 



6.173/2 



(fn 



-1/2-1 



r.,(6.79) (eV) 



Mukhamedzhanov et al. { 2003 1 



Bertone et al. 



(20021 



2001 



Bertone et 

Yamada et al. (2004!_ 
Schiirmann et al. (2008) 



recommended value 



7.4 ±0.4 
7.9 ±0.9 



7.4 ±0.5 



4.9 ±0.5 
4.6 ±0.5 



4.8 ±0.5 



0.47 ±0.03 
0.45 ±0.05 



0.47 ±0.03 



0.53 ±0.03 
0.51 ±0.06 



0.53 ±0.04 



0.41 



+0.34C 
-0.13 

> 0.85 



"channel spin 1 = 3/2 
''channel spin 1 = 1/2 

'^the quoted uncertainty represents a 90% confidence limit 



count; 

ii) limiting the data sets to £' < 1.2 MeV; and 

iii) same as i), but introducing an unidentified = 
5/2- pole at £; = 6 MeV. 

In each case the ANC values and the radii were deter- 
mined. The results for the three cases are 

i) C6.79 = 4.61 ± 0.02 fm-1/2 for a = 4.14 fm and 
S6f|(0)=l.ll keV b. This solution has the lowest 



X but was rejected for the reasons given above. 



ii) C 



6.79 



4.65 ± 0.02 fmi/2 for 
= 1.15 keV b. 



iii) C6 79 = 4.69 ± 0.02 fm^/^ for 
Sf -^Ko) = 1T8 keV b. 



4.6 fm and 



= 5.4 fm and 



The latter two fits are in very good agreement with Run- 
kle et al. ( |2005 ) and about 5% lower than Imbriani et al. 
( |2005 ). An three fits are shown in Fig. [l2 

In summary, the dominant systematic uncertainty for 
^114(0) arises from the interpretation of the high-energy 
data. This uncertainty is estimated from cases i) to iii) 
to be about 4%. One could speculate that the deviation 
of the higher energy data from the R-matrix fit is due to 
broad unidentified structures in this transition (Fig. 12 1. 



We recommend Sfjf{0) = 1.18 ± 0.05 keV b. The er- 
ror includes both systematic and statistical uncertainties, 
though the former are much larger. 

The weighted mean of the ANC for the 6.79 MeV state 
from indirect measurements, C6.79 = 4.8 ± 0.5 fm-^/^ 
(Table |x| , is in excellent agreement with the results of 
the R-matrix analysis. 

Ground state transition: Three data sets (Imbriani 
\et al, 20(55| [Runkle d] |2005[ [Schroder et al[|1987[ ). 



normalized to a;7259 as discussed above, were used in 



the gr ound-state analysis. The results from |Marta et al. 
( 2008 ) - three data points with high precision above the 



259 keV resonance and essentially free from summing ef- 
fects " are relative to the yield of the transition to the 6.79 



J3 

> 



LU 
CO 




FIG. 12 (Color online) R-matrix fits to t he ^'^N(p,7)^^0 
6.79 MeV transition together with the data of |Schroder et al\ 
(19871 (open squares), Imbriani et al. (20051 (open triangles) , 
and Runkle et al. ( 2005 1 (open circles). The cases i, ii, and 
iii (see text) are represented by the dotted green, dash-dotted 
blue, and dashed red lines, respectively. The black line is a 
calculation similar to iii), but without the unidentified J'^ = 
5/2" pole at E=& MeV, comparable to fits in past work. 



MeV state. These data were normalized to the weighted 



mean of the renormalized S-factor (see Sec. 


XI.A.2 


) from 


Schroder et al. 


(1987 


1, 


Runkle et al. 


(2005), and|Imbriani 



et al. (2005) in the energy region 311 keV < < 360 



keV. 

The R-matrix fit was based on the same poles as in 
Angulo and Descouvemont (2001 1 with starting parame- 
ters as given in Ajzenberg-Selove ( 1991 ). The sensitivity 



to radius was tested for a broad range of ANC values, 
6 fm-^/^ < Cgs3/2 < 9 fm^^/^. The minimum was 
obtained for a = 5.6 ± 0.1 fm. Thus, we selected a = 
5.5 fm as an appropriate average for the ground and 6.79 
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MeV states, employing this value for all subsequent R- 



et al. (20051, 


Imbriani et al. 


2005), 


Marta et al. { 


2008 


1, 


and 'Mukhamedzhanov et al. 


(12003). The reduced width 



for the subthreshold state was fixed through Ce.rg (see 
above) to 7^=0.37 MeV. The narrow resonances at 0.987 
MeV {Tp ^ 



3.6 keV, see Fig. 13) and 2.191 MeV (J'^ % 



5/2^, Tp = 10 keV) are not relevant for Sf\4(0) and 
thus were excluded from the fit. In order to optimize the 
fit ofF-resonance, contributions to from points near the 
2.191 MeV and 0.259 MeV (Tp - 1 keV) resonances were 
omitted. As slopes are steep and counting rates peak 
near the resonances, the inclusion of near-resonance data 
forces the fit in arbitrary ways. The region excluded de- 
pends on resonance width and on target thickness, which 
can spread the effects of a resonance over a larger en- 
ergy interval. We omitted data in the interval between 
E]i — 20T and £'i^-f 1.5A, where A is the target thickness. 
Target thickness effects are especially prominent in the 
data of Schroder et al. (19871, representing the integral 
over the target thickness of ^ 30 keV. 

In the fit the decreases with increasing ANC, reach- 
ing a minimum at C, 



W 0.1 




,3/2 ~ 11 fm a value out- 

side th e ra nges determine d by [M ukhamedzhanov et al.\ 
• ■ ■ Tjl 



(2003) and 



Bertone et al. 



(20021 



FIG. 13 R-matrix fit to the ^'^ N(p,7)"0 ground st ate transi- 
tion. The filled circles are from Marta et al. (20081. All other 
data are labeled as in Fig. |12[ 



Most recent treatments of "'^^N(p,7)-'^^0 direct measure- 
ments have failed to address issues connected with the 
total radiative width. 



At the 9 fm up- 
per bound for Cgs3/2, we obtain Sf\4(0) = 0.29 keV, 
while at the 6 fm'^/^ lower bound, Sfi4(0)== 0.24 keV 
b. These fits do not include the possibility of a small 

contribution from Cg3i/2, interfering with the 259 keV 4. Transition to tlie 6.17 MeV state 
resonance. We expand the uncertainty to account for 
such a possibility, recommending Sf\4(0) — 0.27 ± 0.05 
keV b with r^(int) = 1.1 eV. The latter value is the inter- 
nal part of the —0.504 MeV subthreshold state radiative 
width (at E = 0), a fit parameter in the R-matrix cal- 
culation. The total radiative width, which can be com- 
pared to experimental values obtained from, e.g., life- 
time measurements, is derived following the approach of 



This transition was analyzed with the poles given by 
Angulo and Descouvemont (2001 ) except that we also al- 



Holt et al. (1978) and Barker and Kajino (1991), giving 
r^(6.79) = ]r7int)i/2 ± r^(ch)^/^|^ where the relative 
sign of the two amplitudes is unknown. The channel (ex- 
ternal) radiative width r^(ch) = 0.57 eV can be directly 
calculated from the adopted value of Cgg3/2 . If the minus 
sign is chosen in the relationship for r-y(6.79), one obtains 
a lifetime in excess of 4 fs, in disagreement with [Bertone] 
et al.\ ( |2001[ ) and [Schiirmann et al.\ ( |2008[ ). If the plus 



lowed for an external capture contribution (channel spin 
I = 3/2), improving the fit substantially. The primary 
uncertainty in predicting Si'll{Q) arises from the choice 
of the poles, i.e., more poles at higher energies and their 
interference pattern, respectively, could be included in 
the fit. However, a full study of all possible minor contri- 
butions is far beyond the scope of the present work and 
would be hampered by the lack of precise data. The 
best fit yields SljJiO) = 0.13 keV b with 

-1/2 



6.171/2 



-1/2 



sign is chosen, a lifetime shorter than 0.2 fs is obtained. 
Such a lifetime is presently beyond the reach of Doppler 
shift lifetime measurements, but still in agreement with 
Schiirmann et al. (2008). However, the Coulomb excita- 
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0.43 ± 0.02 fm-i/2 and C6.173/2 = 0.49 ± 0.02 fm 
Th ese ANCs are in good agr e ement with those deduced 
by Mukhamedzhanov et al. (2003) and Bertone et al. 
(2002) (see Table |Xp. Previous results without the con- 
tribution from chann el spin 3/2 external c apture led to 



Sf ll(O) = 0.08 keV b d lmbriani et al. 2005|) and 0.04 keV 



b (|Runkle et aL[|2005[). Thus, we have adopted SfU{0) 



tion work of Yamada et al. { 2004 ) gives a lower limit of 
0.4 fs, apparently ruling out such a short lifetime. We 
conclude that the current experimental situation is un- 
satisfactory and calls for further work. Lifetimes larger 
than 0.4 fs require 0^^3/2 < 6 fm~^/^, again in disagree- 



= 0.13 ± 0.06 keV b where the error refiects the uncer- 
tainty in the R- matrix input as well as the spread of this 



value in the literature ( 


Angulo and Descouvemont 


2001 


Imbriani et al. 2005 Nelson et al.[ |2003| |Runkle et al. 


2005 


). In Nelson et al. 


(2003 


) a Ml contribution was in- 



ment with Bertone et al. (2002) and Mukhamedzhanov 



ferred from an analyzing power experiment. The fit only 
extends to £^ ~ 327 keV and trends above the data for 



et al. (2003). The somewhat larger range in Cgs3/2 used 
m tlie present analysis, compared to the uncertainty rec- 
ommended in Table pCl takes account of this dilemma. 



higher energies. Runkle et al. (2005 1 showed that there is 
no significant difference in S^ ll (0) results from including 
the A/1 contribution specified by Nelson et al. (2003). 
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TABLE XI 8114(0) and the fractional uncertainty ASii4(0) for the different transitions. Note that tr(5.24)— >-0 includes 
contributions from the transition tr->-6.86->5.24 and tr-5>7.28->5.24 with 8114(0) = 0.037 ± 0.011 and 0.019 ± 0.006 keV b, 



respectively (from Schroder et al. (1987| with a 30% uncertainty). The contribution of tr(7.28)— >-0 observed by |Schroder et al. 
(|1987[) is negligible. 



transition 



8114(0) (keVb) A8ii4(0) 



reference 



tr-^0 

tr->-6.79 

tr^6.17 

tr^5.18 

tr(5.24)-^0'' 



0.27 ±0.05 
1.18 ±0.05 
0.13 ±0.06 
0.010 ±0.003 
0.070 ±0.021 



19% 
4% 
38% 
30% 
30% 



present 
present 
present 



"value from the analysis of the secondary transition 
''uncertainty from the R-matrix analysis only 
'^from normalization to CJ7259 



Imbriani et al. 




20051 


Imbriani et al. 




20051 



R-Matrix sum 


1.66 ±0.08' 


5% 


additional systematic uncertainty'^ 




5% 


total 


1.66 ±0.12 


7% 



5. Total Si 14(0) and conclusions 



We have obtaine d S^°'^4(0) from the d ata se ts of Im- 
briani et aL\ pOOSl, IMarta et al.\ (|2008|), and [Schroder 



et al. (1987), normalized to the 259 keV resonance, and 
supported by an R-matrix analysis that defines the ex- 
trapolation to astrophysical energies. The R-matrix anal- 
ysis focused on the systematic uncertainties associated 
with fitting and extrapolating the data, and made use 



of indirect measurements (Bertone et al. 2001 



Mukhamedzhanov et al.\ |2003| [Schiirmann et al 



2002 



2008 



Yamada et al. 2004 1 to constrain parameters in the fit- 
ting. Systematic uncertainties in this analysis dominate 
the errors: statistical uncertainties have minor conse- 
quences for the resulting S'°f4(0). The R-matrix radius 
a is a key parameter, fixed in the present analysis to the 



best-choice value of 5.5 fm (Sec. X1.A.3 1. The extrapola- 
tion for the strongest transition to the 6.79 MeV state is 
robust within 4%, while the extrapolations for transitions 
to the ground and 6.17 MeV states are less constrained. 
The transitions to the 5.18, 5.24, 6.86, and 7.28 MeV 
states combine to contribute 0.08 keV b to S'°i4(0), ^ 
5% of the total. These contributions were obtained from 



literature ( Imbriani et al. 2005 Schroder et al.\ 19871, 



scaled to the weighted mean of L0'y259- The errors on the 
individual transitions were enlarged to a more realistic 
uncertainty of 30%. Note that some of the weak tran- 
sitions often have been neglected in past work. Finally, 
an additional systematic error of 5% due to the normal- 
ization of CL17259 (see Table IX) is included. Table XI 
summarizes the various contributions. 

We find, after summing all contributions, S'i°]'4(0) = 
1.66 ± 0.12 keV b. The S-factor fits derived in the present 



study are shown in F igs. [T2l[T3l and 



T4| together with the 



renormalized data of |Imbriani et al. 


(2005 




Marta et al. 


(2008) 


, Runkle et al.\ ( 


2005 


, and 
ilts for 


Schroder et al. 


(1987 


1. 


Figure 


15 compares our rest] 


the total S^°|'4(£;) with 



the data from Lemut et al. (2006) and Bemmerer et al. 




0.6 0.£ 
E [MeV] 

FIG. 14 (Color onhne) R-matrix fit to the "N(p,7)^''''0 6.17 
MeV transition. The symbols are the same as in Fig. |12| The 
dotted green line corresponds to the present analysis. The 
solid black, das hed red, and dash-dot t ed black lines are th e 
R-matrix fits of Imbriani et al. (2005), Runkle et al. (20051, 



and Nelson et al. 



( 2003 1 , respectively. 



(2006b). Below E — 108 keV the gas-target results and 
the R-matrix fit are not inconsistent, given uncertainties; 
at higher energies, E ^ 200 keV, the average deviation is 
~ 8%. These data are an absolute determination of the 
S-factor and thus do not depend on the normalization of 

W7259- 

S*°*4(i?) below E ^ 130 keV can be approximated to 
better than 1% by a second order polynomial 

sr*4(o) 



*°14'(0) 



1.66 keV b 
-0.0033 b 

4.4 X 10"^ b/keV. (50) 
The absolute scale of this energy dependence has an un- 



S*i' 

sri4"(o) 
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FIG. 15 Comparison of the S5°i4 obtained from the present R- 
matrix fit and gas target data. Note that the gas target data 
are corrected for electron screening (see Table 2 in Bemmerer 
et al. (2006b I ) according to calculations of Assenbaum et al. 



certainty of ± 7%. Recently, a coupled channel analysis 
of the data for ^^ N(p,7)^^0 has been reported (Grinevl 
ciute et al. 2008 1 which gives S*' 



^*°l*4(0) 



1.68 keV b, in 

excellent agreement with the results presented here. 

Further work on ^'*N(p,7)-'^^0 is needed. A better 
understanding of the reaction mechanism governing the 
transition to the 6.79 MeV state at high energies would 
help reduce systematic uncertainties. Moreover, addi- 
tional experimental and theoretical work on the transi- 
tion to the 6.17 MeV state is needed, as the existing 
database is lacking. A new determination of for the 
6.79 MeV state with an alternative method would be de- 
sirable to constrain the R-matrix fit and to resolve slight 
discrepancies in existing data. Elastic scattering exper- 
iments could give an additional constraint. Finally, a 
high-precision measurement of 1UJ7259 with significant im- 
provements in the accuracy of stopping power data would 
reduce the systematic uncertainty in the normalization.^ 



^ Note added in proof: A new R-matrix analysis of ^*N(p,7)^^0 
reaction appeared ( [Azuma et al.\ \20W\ after submission of the 
present work. This analysis, which served as a validity test for 
the AZURE code, yielded S*°j'4(0) = 1.81 keV b, 9% larger than 
the central value recommended here. No unc e rtaint y was pro- 
vided. The differences between [Azuma et aL| ( |2010[ l and Solar 
Fusion II are connected with the 6.79 MeV transition. In the 
present work (i) a normalization procedure is employed to ad- 
dress needed corrections in the high-energy data and (ii) a back- 
ground pole is introduced to achieve a better representation of 
that data. Without such adjustments, the procedure of |Azurna| 
I et all\ l |2010| produces a fit that underestimates the high-energy 
data and consequently yields a larger Sj'^J'(O). Nevertheless, the 
present and [Azuma et al.\ { 2010| l results are consistent if one 
assigns a reasonable uncertainty to the latter. 



B. Other CNO-cycle reactions 

While the ^^N(p,7)^^0 reaction controls the cycling 
rate and the energy production by CN reactions at so- 
lar temperatures, other reactions in the cycle determine 
the extent to which the reaction flow moves out of the 
CN cycle toward heavier metals, oxygen in particular. 
These trends in turn affect the opacity evolution and 
temperature profiles as a function of solar age. There 
has been significant recent progress in determining the 
rates of many of these other reactions. The reader is 
referred to Solar Fusion I for summaries of other reac- 
tions for which there has not been new work reported 



gulo et al. 


( 


19991 (th 


Wiescher et al. 


pOlO) 



1. ^^C(p,7)"N 

In the starting phase of the CN cycle, before it has 
reache d its equilibrium, this reactio n controls the buildup 
of ^^N ( Haxton and Serenelli |2008 1 . A recent study using 



the ANC method by [Burtebaev et al.\ (|2008|) yields a 
reaction rate consistent with that of Angulo et al. { 1999 1, 



the rate recommended here. 



2. ^^'N{p,aY'^C 



^^N(p, a)^^C reaction competes with 
, a parallel study of the two is highly 



As the 
i5n(p, 7)160 

desirable. In Solar Fusion I, a weighted average of 
^115(0) = 67. 5 lb 4.0 MeV b was r ecom mended using 
the results of iRedder et al.\ { 1982 ) and Zyskind and 

-^° N(p, a ) -^^C r eaction has 
(|2007D, using the 



Parker (1979). Recently the 



been measured by La Cognata et al 



indirect Trojan Horse Method (TH method) (see Sec. 
XII). The new data have been analyzed along with 
''^N(p, 7)1^0, using a common R-matrix approach. The 
TH method allows one to extend the explored energy 
range down to about 20 keV, without the complication 
of electron screening enhancements that enter for direct 
measurements. Thus the TH measurements provide 
complementary information that can be helpful in 
checking the overall consistency of S-f actor fits. |La| 
Cognata et d] ( [2007| determined S?i5(0) = 68±11 MeV 
b from TH measurements. New R-matrix fits to direct 



data of [Redder et a/.| ( |1982[ ) by |La Cognata et aL| ( |2009[ ) 
yielded S?i5(0) = 73 ± 5 and 74 ± 9 MeV b, depending 



on the resp ective energy ranges fit (see La Cognata 
TTai] ( [2009 1 for details), and 8^15(0) = 70 ± 13 MeV b 



for the indirect TH method data of 



La Cognata et aL\ 



( [2007[ ). An R-matrix fit by [Barker] (|2008af , which did 
not include the TH method results, gave S"i5(0) = 80 
MeV b. We recommend the value Sf i^(0) = 73 ± 5 



MeV b obtained by La Cognata et al. ( [2009 ) by fitting 
direct data as the new best value for the ^^N(p, a^^C 
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FIG. 16 Summary of the available measurements of 8115(0), 
showing values as originally reported on the dates indicated. 
The shaded band corresponds to the NACRE compilation 



(Angulo et al. 19991. From La Cognata et al. (20091, by 



pernussion. 



reaction (see Table XII). It is consistent with the two 
direct measurements, the indirect TH method data, 
and the R-matrix fit by Barker. A summary given by 
La Cognata et al. (2009) of S"i5(0) determinations is 
shown in Fig. |1^ In Table |XII| the deriva tives shown 
are those reported by Zyskind and Parker (1979), and 



may therefore not be completely consistent with the 
R-matrix energy dependence calculated by |La Cognata 
erall(|2009|. 



3. ^5N(p,7)1«0 

The ^^N(p,7)^^0 reaction provides the path to form 
^^O in stellar hydrogen burning^, thus transforming the 
CN cycle into the CNO bi-cycle and CNO tri-cycle. In 
stellar environments, the reaction proceeds at very low 
energies, where it is dominated by resonant capture to the 
ground state through the first two interfering = 1^ 
s-wave resonances at En = 312 and 964 keV. In addition 
there is some direct capture to the ground state. Direct 



measurements have been reported by Hebbard ( 1960 ) for 
proton ene rgies down to 220 keV and by [Rolfs and Rod-| 
ney (1974) down to proton energies of 155 keV. These 
measurements disagree significantly below 300 keV. In 
order to fit their low-energy data, [Rolfs an d Rodney] 
(1974) included the interference of the two 1^ resonant 



capture amplitudes with the nonresonant (direct) com- 
ponent to the ground state of ^^O calculated in the hard- 
sphere approximation. The absolute normalization of the 
direct term is entirely determined by the ANC of the 



® Most of the ^^O found in the Sun originates not from hydrogen 
burning in the Sun itself, but instead from the ashes of helium 
burning in earlier stars. 
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FIG. 17 (Color online) S(0) for the i^N(p,7)"0 reaction. 



Data from Hebba rd ( 1960 |) (blu e circles, limited to E ^ 210 
keV), Rolfs and Rodney, ( |1974 1 (green triangles), and Bem- 
merer et al. (2009) (red squares). Error bars reflect statistical 
and systematic uncertainties summed in quadrature. Dashed 
line, previous R-matrix fit and shaded area, its quoted 17% 
uncertainty, from iMukhamed zhanov et al. ( 2008 1. Dotted 
line: previous extrapolation by Angulo et al. (19991. Figure 



from Bemmerer et al. (20091, by permission 



bound state for -^^N-l-p — » ^^O. The spectroscopic factor 



adopted by Rolfs and Rodney (1974) corresponds to an 



ANC almost an order of magnitude l arger than the one 
determined from ^^N(^He,d)^^0 by Mukhamedzhanov 



et al. (2008) 



A new 
two-level, 



analysis of 
two-channel 



the direct data using the 
R-matrix was presented by 
Mukhamedzhanov et al. (2008). The contribution from 
the a — ^^C channel was also taken into account. The 
determined astrophysical factor Si]^5(0) = 36 ± 6 keV 
b is about a factor of two lower than the previously 
accepted value S7 



(21 



64 ± 6 keV b from Rolfs 



and Rodney (1974). Hebbard (1960) reported 



32 ifc 6 keV b at 23.44 keV, which was converted by 
Mukhamedzhanov et ai\ (|2008| to S7i5(0) = 29.8 ± 5.4 



keV b using the polynomial extrapolation given by Heb- 
bard. Mukhamedzhanov et al. (20081 conclude that for 



every 2200 ± 300 cycles of the main CN cycle, one CN 
catalyst is lost due to this reaction, rather than 880 cycles 
recommended by [Rolfs and Rodney ( 1974 ) and 1000 cy- 



cles re commended by the NACRE compilations (Angulo 
et al. 1 1999 J. Their result coincides with the R-matrix 



analysis by [Barker (2008b), which yielded a leak rate of 
1/2300. Barker's analysis was completed before the ANC 
data were available and shows a larger spread of S values. 

New measurements of this reaction at LUNA by [Bem 
merer et al. (2009) yielded cross sections with improved 



precision for energies between 90 to 230 keV (Fig. 17). 
The extent of the agreement between the new LUNA 
data and the Hebbard data point to a possible uniden- 



Rolfs and Rodney ( 


19741. The value S7i 


M = 


keV b obtained by IV 


ukhamedzhanov et al. 


(20081 



36 ± 6 
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regarded as an interim recommendation pending an up- 
dated analysis taking full account of new data (e.g., com- 
pletion of the analyses for recent LUNA and Notre Dame 
experiments). Further measurements at higher energies 
are also desirable in order to constrain the R-matrix fits. 



et al. (2007). The large uncertainty (~ 50%) makes a 
new round of measurements of the direct capture cross 
section desirable.'' 



7. i»0(p,a)i''N 



4. ^«0(p,7)^^F 

The cross section is dominated by direct capture to the 
ground and first excited states of ^ ''F. Because the latter 
is weakly bound, its S-factor rises rapidly at low energies 
and the ground-state transition plays a minor role. Cal- 
culations of the direct capture process by Rolfs (19731, 



Morlock et al. (1997), Baye et al. (1998), and Baye and 



Brainis ( 2000j ) give a quantitativ e account of the en ergy 
dependence of both transitions. Baye et al. ( 1998 ) cal- 



culate S7i6(0) with two choices for the nuclear force, ob- 
taining S7ig(0) = 10.2 and 11.0 keV b when normalized 



to the data of |Rolfs| ( |1973[ ) and |Morlock et al.\ ( |1997[ ) . The 
value adopted here is S/;^g(0) = 10.6 ± 0.8 keV b and the 
derivative is S7i fi(0) = —0.054 b. A recent reevaluation 



by Iliadis et al. using both R-matrix theory and 

a potential model yielded reaction rates at temp eratures 
> lO^'K that are consistent with these values (Angulo 
et al. 1999), but with a lower assigned uncertainty. 



5. i^O(p,a)^*N 

The i^O(p,a)i**N reaction closes branch II of the CNO 
bi-cycle. The reaction rate at solar energies is dominated 
by a subthreshold resonance at Eji — —3.1 keV and a 
resonance at Eji — 65.1 keV. Several recent experiments 
have clarified the strength and location of a 2" resonance 
at 183.3 keV that plays a significant role at the higher 
temperatures characteristic of novae and asymptotic gi- 



ant branch stars ( Chafa et al. 


2005 


2007 


Moazen et al. 


2007 




Chafa et al. ( 2007 ) find a low-energ 


y cross section 



gulo et aL| ( |1999[ ), refiecting a re-evaluation of the proton 
width of the subthreshold resonance. No calculated value 
for S"i7(0) has been published. 



The ^^0-|-p interaction represents a branching point in 
the CNO cycle: the ^^0(p, a)^^N reaction leads to a re- 
cycling of CN catalytic material, while ^^0(p, 7)^^F may 
lead to a loss of this material, depending on the fate of 
the produced ^^F. Nine resonances below 1 MeV influence 
the astrophysical rate for ^^0( p, Q!)^^N, with those at 20, 
144, and 656 keV dominating ( [Angulo e^"al^ |1999[ ). The 
presence of strong resonances in the astrophysical regime 
makes extraction of a value for S"]^g(0) inappropriate. 

The strength of the 20-keV resonance had been known 
only from spectroscopic measurements performed by 



Champagne and Pitt ( 1986 ) through the transfer reac- 
tion ^''OCHe, d)^^F and throug h the direct capture re - 
action 



0(P>7)^^F measured by |Wiescher et a^.| ( |1980 |. 
The cross section at 20 keV is a factor ~ 10^^ smaller 
than the one at 70 keV owing to the Coulomb barrier 
penetration factor. This makes a direct measurement 
of the cross section impossible with present-day nuclear 
physics facilities. Furthermore the spin and parity of the 
8.084 MeV level in ^^F (corresponding to a 90 keV reso- 
nance in the ^*0(p, a)^^N cross section) was not known. 
In order to reduce the nuclear uncertainties affecting the 
reaction rate, which |La Cognata et al. ( 2008) estimated 
at about an order of magnitude, a new round of mea- 
surements has been made with the TH method by |La| 



Cognata et aL| ( |2008[|2009[|2010[ ). The deduced strength 
of the 20 keV resonance u;7 = 8.3;2:6 ^ 10"^^ eVelimi- 
nates much of the broad range given by NACRE (Angulo 

19 ev 



et al. 



19991, W7 = 6_5 X 10 ^ eV, and decreases the 
uncertainty of the reaction rate by about a factor 8.5 (La 
Cognata et aLl i2008, ,20091 120101). In addition, the spin 



;3/2-) and strength of the 90-keV resonance, which was 
seen in the work of Lorentz-Wirzba et al.' (1979 ) 

(2008, 2009 



determined. The La Cognata et al. 



and Lorentz-Wirzba et al. ( 1979 ) strengths agree. 



were 



2010) 



6. ^^0(p,7)^«F 

The cross section shows a number of resonances in the 
range relevant to the hot CNO cycle in novae. Effort has 
been recently devoted by Chafa et al. (2005 2007) and 



Fox et al. ( 2004 ) to measuring the resonance parameters 
inboth^'0(p77)iSF and i^O(p, a)"N. While the higher- 
lying resonances are not directly relevant to solar CNO 
processing, they do have a significant influence in modern 
interpretations of the work of Rolfs ( 1973 1 , who measured 



the direct capture cross section that dominates at solar 



energies. |Fox et al.\ ( |2005[ ) and [C hafa et al.\\2001\ both 
concluded that significant corr ection s are require d. The 
recommended S7i7(0) in Table XII is taken from Chafa 



XII. INDIRECT METHODS AND THEIR VALIDATION 

Three classes of experiments contribute to our under- 
standing of solar fusion reactions, direct cross section 
measurements, indirect methods, and ancillary nuclear 
structure techniques for determining the properties of 
resonances (energies, 7 and particle widths, and spins 
and parities). Indirect methods involve the use of nu- 
clear reactions related to, but not identical to, the solar 



Note in proof: The direct capture cross section was recently ex- 
tracted fr om new measuremen ts between lab energies of 193 and 
519 keV (INewton et a«.||2010||. 
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TABLE XII Summary of updates to S-values and derivatives for CNO reactions. 



Reaction 



Cycle 



S(0) 
keV b 



S'(0) 
b 



S"(0) 
keV"^ b 



References 



^^C(p,7)"N 



14 



N(p,7)''0 



1.34 ± 0.21 
7.6 ± 1.0 
7.0 ± 1.5 

1.66 ± 0.12 



2.6x10"- 
-7.83x10" 



-3.3x10"^ 



8.3x10 
7.29x10 

4.4x10" 



-5 
7" 



Recommended: Solar Fusion I 
Recommended: Solar Fusion I 
NACRE: 



Angulo et al. 



(19991 



Recommended: this paper 



^^N(p,Qo)'^C I 
^5N(p, 7)1^0 II 



(7.3 ± 0.5) X 10"* 
36 ±6 
64 ±6 
29.8 ±5.4 



351 



11 



Recommended: this paper 



Mukhamedzhanov et al. 



Rolfs and Rodney 



Hebbard ( 1960 1 



( 1974 1 



(20081 



^«0(p,7)^'F 


II 


10.6 ±0.8 


-0.054 


Recommended: this paper 




II 




Resonances 


Cliafa et al. 


(2007) 




III 


6.2 ±3.1 


1.6x10"^ 


-3.4x10"'^ Cha.ia et al. 


(2007) 



'0(p,a)^^N 



III 



Resonances 



See text 



^0(p,7)"F IV 



15.7± 2.1 



3.4x10" 



-2.4x10" 



Recommended: Solar Fusion I 



reactions under study, as tools to probe properties of the 
solar reactions. References have been made in this review 
to three indirect methods, asymptotic normalization co- 
efficients, Coulomb dissociation, and the Trojan horse 
method. As the connection between the indirect observ- 
able and the solar reaction of interest must be established 
through reaction theory, such methods entail a greater 
degree of model dependence, impacting systematic un- 
certainties. But indirect methods also have many virtues: 
they can be applied when direct measurements are diffi- 
cult or impossible, have systematic uncertainties that are 
different from those of direct measurements, and provide 
supplementary information that can constrain R-matrix 
and other models used in the extrapolation of data from 
direct measurements. The role of indirect measurements 
in validating and constraining models is apparent from 



asymptotic normalization coefficient (ANC) is defined as 
the amplitude of this component (apart from an over- 



the discussions, for example, of Sec. XI. A 



A. The asymptotic normalization coefficient method 

The asymptotic normalization coefficient method con- 
strains S(0) by exploiting the peripheral nature of many 
radiative capture reactions in nuclear astrophysics. Be- 
cause of Coulomb and/or centrifugal barriers, most (p,7) 
and (a, 7) reactions are peripheral at solar energies. The 
cross section for a nonresonant radiative capture reac- 
tion A(p, 7)B at zero relative energy depends only on the 
long-distance behavior of the p+A wave function (and 
on the overlap of that extended wave function with B). 
The detailed short-range behavior of the scattering state 
p±A or bound state B, governed by the strong interac- 
tion and nuclear length scales, are not relevant to the 
reaction mechanism. The bound-state wave function at 
long distances will contain a component corresponding 
to two separated clusters, p and A, with the cluster rel- 
ative radial motion given by a Whittaker function. The 
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. A distinct ANC will govern the nonresonant 



capture into each final state, i.e., the ground or bound 
excited states of B. Therefore, if one can identify another 
nuclear reaction that includes the vertex A ± p B and 
is sensitive only to the tail of the radial overlap function, 
the needed ANC can be determined from that reaction. 
This measurement in a different system then determines 
the radiative capture cross secti on at zero relative energy 
( [Mukhamedzhanov et al. 2001), up to small corrections 
determined by the scattering wave function and the po- 



tential in the continuum ( Capel and Nunes 2006 Typel 



and Baur[ |2005|). While the method is limited to S(0), 



providing a data point below the Gamow peak, this often 
complements the data from direct measurements, which 
are frequently limited to energies above the Gamow peak. 

In most applications, the ANC is deduced from trans- 
fer reactions. The extraction relies on the distorted wave 
Born approximation (DWBA) and the direct proportion- 
ality between the transfer cross section and the square of 
the ANC. Provided that the transfer reaction is com- 
pletely peripheral and the measured angular distribu- 
tions are well described within the single-step DWBA, 
the ANC can be extracted. The main source of uncer- 
tainty comes from the optical model description, typically 
> 10% for reactions above the Coulomb barrier. For this 
reason, it is often important to also measure the elas- 
tic channel of the corresponding transfer reaction over 
a wide angular range, to help constrain optical model 
parameters. Investigations of effects beyond the single- 
step DWBA arising from target excitation suggest that 
deformed targets with strong couplings to low-lying ex- 
cited states are not good candidates for the ANC method 
(Azhari et al. 2001). Some of the applications of the 
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method involve loosely bound nuclei, opening up the pos- 
sibility of multi-step processes through continuum states 
as viable alternatives to the direct reaction mechanism. 
So far there has only been one reaction for which the 
magnitude of this effect has been evaluated; in this case 
it was found to be negligible ( |Moro et al. 20031, but a 
more systematic study should be done. 

In Solar Fusion I the ^^0(p,7)^^F reaction was identi- 
fied as a good test for the method. As a consequence, the 
i'^0(^He,d)^'^F reaction was measured at 30 MeV. The 
angular distributions of the ground state and the first 
excited state were well described within the DWBA and 
the inferred S factors agreed with the radiative capture 
data to better than 9% ( Gagliardi et al. 1999 ) . 



Finally, it should be mentioned that breakup reac- 
tions B + T ^ A + -p + T can also be used to extract 



There have been many subsequent applications of this 
method, mostly involving peripheral transfer reactions 
on intermediate mass targets. Here we focus on those 
relevant to validating the method for solar fusion re- 
actions. Two transfer reactions, ^''B(''Be,®B)^Be and 
^4N(^Be,^B)iy , were u sed to extract the ANC for Si7(0) 
(Azhari et al. 1999alb|. For both targets, the peripheral 



nature of the transfer reactions were checked carefully 
by evaluating the sensitivity of the extracted ANC to 
the single particle parameters of the binding potential in 
the DWBA analysis. Similar analyses have been done 
by inv oking a radial cutoff in the distorted wave calcu- 
lation ( Fernandez et al.\ '2000^ iMukhamedzhanov et al. 



1997 ). In Tabacaru et al. (2006) a joint analysis was per- 
formed, yielding Si7(0)=18.0 ± 1.9 eV b, which can be 
compared to the best value from direct measurements, 
20.8 ± 0.7 ± 1.4 eV b. In additio n, the low-energy reac 



tion7Be(d£i)*B at Euh=7-5 McV ('Liu et aL[|T996l|Ogata 



et al. 2003 1 was studied, but difficulties were encoun- 



tered in the analysis. The (d,n) reaction model depends 
on the poorly constrained exit-channel neutron optical 
potential. In addition, the use of low energies, neces- 
sary to satisfy the peripherality condition given the low 
Z of the deuteron, leads to significant compound nuclear 
contributions, introducing additional uncertainties. 

This review includes several illustrations of the use 
of ANC determinations to validate R-matrix descrip- 
tions of direct reaction data. In Sec. IXI.AI the 
example of the subthreshold-state (6.79 MeV) con- 
tribution to ^''N(p,7)^^0 is described in some de- 
tail: the ANC determined from the R-matrix fit is 
in good agreement with that extracted by |Bertone| 



et al. ( [2002^ and ^Mukhamcdzhanov et al. ( 2003 ) from 



N(^He;d)^0. Analogous work using ^'^"Npfl e^dr^O to 



XI.B 



IT 

study ^^N(p,7)^^0 is discussed in Sec. 

As ANCs can be related to spectroscopic factors, the 
latter can also be used to parameterize cross sections. 
However, spectroscopic factors have an additional depen- 
dence on the single-particle bound state orbitals assumed 
in their extraction. Consequently radiative capture re- 
actions parameterized through ANCs and spectroscopic 
factors have somewhat different uncertainties. Further 



and Bertone et al. (20021 



discussion can be found in Mukhamedzhanov et al. (2001 ) 



ANCs when they meet the peripherality condition (Tra- 
che et al. . 2004) . However a detailed study of the uncer- 



tainties involved in the reaction theory has not yet been 
completed. 



B. The Coulomb dissociation method 

Coulomb dissociation (CD), originally proposed as a 
method for extracting information on astrophysical fu- 
sion cross sections by Rebel, was developed theoretically 

1986j). The process oc- 



shortly thereafter (Baur et al 



curs when a beam of fast projectiles interacts with a 
heavy target such as Pb. An energetic virtual photon 
from the target can then dissociate the projectile, lib- 
erating a nucleon or a particle. To the extent that the 
experimentalist can exploit the kinematics of this process 
to enhance the contributions from the long-distance ex- 
change of single photons, this process can then be related 
by detailed balance to the corresponding radiative cap- 
ture reaction. But several effects complicate this simple 
picture. Whereas nonresonant radiative captures gen- 
erally proceed almost exclusively by El transitions, the 
strong E2 field in CD can be important. Moreover, the 
simple radiative capture/CD correspondence is compli- 
cated by multiple photon exchange and by the strong 
interaction, which can lead to nuclear diffraction dissoci- 
ation and Coulomb-nuclear interference. Strong interac- 
tion effects can be reduced by restricting measurements 
to small angles, where long-range electromagnetic transi- 
tions dominate nuclear interactions. Multiple photon ex- 
change (also known as post-acceleration) can be reduced 
by increasing the beam energy, shortening the time the 
projectile spends in the target's field. 

In Solar Fusion I a proposal was made to test the va- 
lidity of the CD method quantitatively through compari- 
son with a corresponding radiative capture measurement. 
The radiative capture reaction was to have suitable prop- 
erties, including a low Q value, a nonresonant El reaction 
mechanism, reactants with similar mass/charge ratios, 
and a final nuclear state with relatively simple structure. 
Although no perfect reaction was identified, ^Be(p, 7)^B 
appears to be a good choice. Several new measurements 
were made, and a great deal of theoretical effort was in- 
vested in their interpretation and in extracting the S fac- 
tor. This work is summarized in Sec. IIX.CI and will not 
be discussed further here, except to repeat the conclu- 
sion that, while in several cases agreement between the 
CD method and direct measurements has been demon- 
strated at the 10-20% level, remaining uncertainties in 
the magnitude of S(0), in independently determining the 
shape of S{E), and in the theory argue that the inclu- 
sion of CD data in the current S17 evaluation would be 
premature. 

Efforts also have been made to validate the CD method 
for the ^'*C(n,7)-'^^C reaction. Although this reaction is 
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not directly relevant to solar fusion, the radiative cap- 
ture rate is now known to a precision of ~ 10% (Reifarth 
et al.\\2008l The corresponding CD of JfConf^fPbhas 



that is not strongly suppressed by Coulomb barriers or 



recently been remeasured at RIKEN (Nakamura et al. 
|2b09| ) . Reaction models predict that the ^'^C breakup has 
an insignificant nuclear contribution and is dominated by 
El transitions, provided the analysis is limited to events 
in which the ^^C center-of-mass scattering angle and the 
relative energy of the breakup fragments are small. In- 



dependent^^ data ( Esbensen 2009 Sum- 



mers and Nunes 20081 find that the neutron capture 



cross section extracted from CD agrees very well with 
the direct measurement and has comparable precision. 
This appears to be a favorable case for the theoretical 
treatment due to the dominant nonresonant El reaction 
mechanism, small E2 and nuclear contributions, and rel- 
ative simplicity of ^^C, which can be described reason- 
ably in a single-particle ^'*C-|-n potential model. While 
the agreement in this case is promising, some caution 
is warranted because the radiative capture measurement 
has not been confirmed by an independent measurement. 

The ANC and CD methods are both well suited to 
measurements with low intensity radioactive beams be- 
cause the transfer reaction and CD cross sections are 
much larger than the corresponding radiative capture re- 
actions. Moreover, they are both applicable to radiative 
capture reactions. 



C. The Trojan Horse method 



The Troja n Horse (TH) method ( |Baur[|1986[|Spitaleri 
et al. 20041 is an indirect technique to determine the 
astrophysical S factor for rearrangement reactions. It 
allows inference of the cross section of the binary process 



X + b + B 



(51) 



at astrophysical energies through measurement of the TH 
reaction 



a + A 



y 



B. 



(52) 



The measurement is done with quasi-free kinematics, in 
which a TH a having a strong x + y cluster structure is 
accelerated to energies above the Coulomb barrier. After 
penetrating the Coulomb barrier, the nucleus a breaks 
up, leaving x to interact with the target A while the pro- 
jectile fragment y flies away. From the measured cross 
section of reaction ( 52 ), the energy dependence of the bi- 
nary subprocess (51) is determined. While the reaction 



( 52 1 can occur in a variety of ways, the TH reaction mech- 



anism should dominate in a restricted region of three- 
body phase space in which the momentum transfer to the 
spectator nucleus y is small, i.e., quasi-elastic scattering 
conditions apply. Since the transferred particle x in the 
TH reaction ( |52| is virtual, its energy and momentum 
are not related by the on-shell equation E^ = p^/{2mx)- 
The main advantage of the TH method is that the 
low-energy cross sections can be deduced from a reaction 



strongly altered by electron screening ( Assenbaum et al. 



1987 Spitaleri et al. 20011. The TH cross section can 



be used to determine the energy dependence of the bare 
nuclear S factor for the binary process (51) down to zero 



relative kinetic energy of x and A. The absolute value 
of S{E), however, must be determined by normalizing 
to direct measurements at higher energies. To ensure 
quasi-free kinematics one should measure the momentum 
distribution of the spectator fragment y and the angular 
distributions of the fragments of the binary sub-reaction 
to check for contributions from non-TH mechanisms. As 
a check on distortions due to final state interactions, the 
momentum distribution of the spectator can be measured 
and compared with that of the spectator in the free TH 
nucleus (Pizzone et al. 20091. Final state distortions 



can b e treated in DWBA calculations ( La Cognata et al. 
20101). 



The uncertainty of the S{E) extracted from the TH 
method includes contributions from statistics, uncertain- 
ties due to the need to normalize the TH data, finite 
experimental energy resolution, and backgrounds due to 
other reaction mechanisms. The first successful test of 
the TH method was conducte d for the ''Li(p, a)'*He re- 
action (Lattuada et al. 2001). The extracted S(0)= 55 
± 6 keV b includes an uncertainty of 10% from the nor- 
malization of the TH data to the direct data ( [Engstler 



et al.\ [1992 ) and 5.5% from other sources, mainly statis- 
tics. In addition, in Sec. |XI.B.2| we compare results 
for TH and direct determinations of the cross section 
for ^^N(p,a)-'^^C. Although promising, the TH method 
requires further validation by experiment, and its signif- 
icant dependence on reaction theory calls for more in- 
vestigation of the approximations by which TH reactions 
are related to their astrophysical analogs. 

The TH method also provides an important test of 
electron screening potentials, which can be obtained from 
comparisons of direct and TH cross sections. 



D. Summary 

The three indirect techniques discussed here provide 
alternatives to direct measurements of astrophysically 
important reaction rates. In some cases they provide 
the only practical means for determining stellar reaction 
rates. While their connection to solar reactions requires 
an additional level of reaction theory, experimental tests 
of their validity have often yielded agreement with di- 
rect measurements within 10-20%. Significant progress 
has been made since Solar Fusion I in benchmarking in- 
direct techniques. Indirect methods are best applied to 
cases where there is a supporting body of experimental 
data that can be used to constrain the needed nuclear 
model input, such as optical potentials and effective in- 
teractions. 

In actual practice, the distinction between direct and 
indirect methods is not sharp, but rather a matter of de- 
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gree. While a measurement may probe a stellar reaction 
directly, it often does so at a different energy or in a dif- 
ferent screening environment. Thus direct methods also 
depend on reaction theory, to extrapolate data to stellar 
energies or, in cases like S33 where data in the Gamow 
peak have been obtained, to correct for the effects of 
screening in terrestrial targets. Still, the connection to 
stellar physics is typically much closer. Models play a 
less important role, and increasingly the needed model- 
ing can be done microscopically, as direct measurements 
involve light nuclei. 

For this reason we maintain a distinction between di- 
rect and indirect methods in this review, basing our rec- 
ommendations on results from the former. However, in- 
direct methods have had a significant impact on our anal- 
ysis: they have been used in this review to constrain R- 
matrix fits to direct data and to check the consistency 
of conclusions based on analyses and modeling of direct 
data. 

We recommend extending the benchmarking of indi- 
rect methods against direct methods over a wider range 
of reactions, as more data would be useful in quantifying 
the uncertainties in such techniques. 



XIII. FUTURE FACILITIES AND CURRENT 
CAPABILITIES 

We noted in the introduction to this review the cru- 
cial role nuclear astrophysics experiments have played in 
the development of a quantitative SSM and in motivating 
solar neutrino experiments. We outlined the important 
goals that remain in this field - tests of weak interactions 
and of solar properties that make use of high precision 
solar neutrino measurements, hclioseismology mappings 
of c(r), and detailed solar modeling. There are also a 
host of related problems - Big Bang nucleosynthesis, red- 
giant evolution, the evolution of supernova progenitors, 
and a variety of transient explosive phenomena in as- 
trophysics - where a quantitative understanding of the 
nuclear physics is essential. This chapter deals with the 
experimental facilities that have allowed progress in this 
field, and discusses the instrumental developments that 
will be important if we are to continue a similar rate of 
progress over the next decade. 

The measurements that support the development of 
a quantitative theory of main-sequence stellar evolution 
primarily involve low energy proton- and a-capture reac- 
tions that traditionally have been studied with small ac- 
celerators. The machines must be able to provide proton 
or a beams of sufficient intensity to allow cross section 
measurements near the very low energies of the Gamow 
peak. 

Because low energy charged-particle reaction cross sec- 
tions are small, experiments must be designed for sig- 
nal rates much smaller than background rates associ- 
ated with cosmic rays, the natural radioactivity of the 
laboratory environment, and the induced activity aris- 



ing from beam interactions with target impurities. The 
ambient background can be roughly divided into muons 
and neutrons associated with cosmic rays, and 7 rays and 
neutrons from natural radioactivity (uranium, thorium, 
potassium, and radon from surrounding geology). Today 
most charged-particle reaction measurements for nuclear 
astrophysics are being performed at above ground facil- 
ities, with various techniques then employed to mitigate 
backgrounds. The common technique is passive shielding 
around the detection region. Typically a layered combi- 
nation of lead, copper, and polyethylene is used to re- 
duce 7 and neutron backgrounds within detectors with 
relatively small volumes. But additional strategies are 
available to further reduce backgrounds and thus allow 
measurements at energies nearer those relevant for astro- 
physics, including 

1. use of more sophisticated detector setups with both 
passive and active shielding and with triggers to aid 
in event identification; 

2. measurements in inverse kinematics using recoil 
separators in facilities above ground; and 

3. measurements with direct kinematics using acceler- 
ators that are sufficiently deep underground to sup- 
press penetrating cosmic-ray muons and the neu- 
trons and other secondary activities they induce. 

Passive shielding, active shielding, and coincidence 
gating techniques can enhance event identification and 
significantly reduce backgrounds in above-ground labo- 
ratory environments. As most resonance levels of astro- 



physi cal interest decay via 7-cascades (Rowland et at. 
2002[ ) 77-coincidence techniques can be used to signifi- 
cantly reduce the single-7 background. Q-value gating 
techniques, where only events in coincidence with the 
summing peak of the radiative capture reaction are ac- 
cepted (ICouture et al.\ [20081), can allow one to extend 



measurements to lower energies, but at the cost of a de- 
creased overall counting efficiency due to the coincidence 
requirement. 

Alternative techniques have been developed to reduce 
backgrounds without such losses in detection efficiency. 
Two ideas that have demonstrated their promise are mea- 
surements in inverse kinematics - one detects the reaction 
recoil particles rather than the light particles or 7s of the 
reaction - and measurements in underground environ- 
ments. Below we describe past and current experience 
with these two techniques as well as the future facilities, 
in progress or planned, that would allow these techniques 
to be further advanced. 



A. Inverse kinematics measurements using recoil separators 

In an inverse-kinematics experiment a heavy ion in- 
duces (p,7) or (a, 7) reactions when it interacts in a hy- 
drogen or helium gas target. The projectiles and reac- 
tion products move within a narrow cone in the forward 
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direction. A recoil separator is used to reject the pri- 
mary beam while focusing the reaction products for de- 
tection. The charged recoils can be detected with higher 
efficiency than the 7s produced in conventional proton- 
or a-beam experiments. By detecting the 7s in coinci- 
dence with the reaction products, dramatic reductions 
in backgrounds can be achieved. Existing recoil sepa- 
rator facilities fior nuclear astrophysics experiments in- 



clude DRAGON at ISAC in TRIUMF (jHutcheon et al. 



A clean separation is difficult for recoils with large en- 
ergy spreads, making low-energy solar fusion reactions 
particularly challenging. Recoil separators are therefore 
more typically used for higher energies characteristic of 
helium- or explosive hydrogen-burning reactions. For ex- 
ample, the recoil-separator measurements of S34 at the 
ERNA facility in Bochum were limited to data above a 
center-of-mass energy of 700 keV ( Di Leva et al. 2009[ ) . 



2003 ) , the Daresbury separator at HRIBF in Oak Ridge 



Below this energy the angular divergence of the recoils 
exceeds the angular acceptance of the separator, ± 25 



(Fitzgerald et al. 2005), ERNA at the DTL in Bochum mrad (Di Leva ei aL 2008). 



dRogalla et aL[|2003^nd t he RMS at KUTL in Kyushu, 
Japan (Sagara et al. 2005). 



Recoil separators are not useful for (Q!,n) reactions be- 
cause separator acceptance angles are too small, given 
the momentum transfer in this process. 

Recoil separators present several experimental chal- 



lenges (Rogalla et al. 2003), particularly for the low en- 



ergies important in solar fusion cross section measure- 
ments. At such energies, the energy spread and the an- 
gular aperture are, for most solar fusion reactions, larger 
than the acceptance of any of the recoil separators cited 
above. 

The following conditions must be fulfilled in experi- 
ments on absolute cross sections: 

• the transmission of the recoils must be exactly 
known and should ideally be 100%; 

• the charge-state distribution of the recoil products 
must be known or the reaction must be studied for 



20081; 



all charge states produced (Di Leva et al. 
and 

the interaction region must be well defined. 



Therefore, experiments coming on-line in the near future 
are all planning to use compact high-density gas-jet tar- 
gets instead of extended windowless gas targets. 

Recoils of solar fusion reactions typically have rela- 
tively large emission angles and large energy spreads, 
both of which increase with decreasing reaction energies 
E, when E < Q. The angular distribution of recoils fol- 
lowing emission of capture 7-rays of energy E-y is char- 
acterized by an emission cone half-angle of 



= arctan 



E-y 
P 



(53) 



where p is the momentum of the beam (c = 1). The total 
energy spread AE of the recoil accompanying 7 emission 
is 



AE 

If 



AE.y 

p 



(54) 



Furthermore a large spatial separation between the re- 
action products and the beam is required, as the pri- 
mary beam intensity is typically many orders of magni- 
tude larger than that of the recoiling reaction products. 



Two dedicated next-generation separators for low- 
energy nuclear astrophysics studies with stable ion beams 
will soon come on line, the St. George facility at 



Notre Dame's Nuclear Science Laboratory ( Couder et al. 



2008[ ) and the ERNA separator at the CIRCE facility in 

Caserta, Italy. 

[2003 ). 



The latter is based on a redesign of the 
Bochum ERNA separator (IRogalla et al. 



Both 

separators feature large acceptances in angle and energy 
and will be equipped with high density gas-jet targets to 
ensure well defined interaction regions. Figure [l8| shows 
the layout of the St. George recoil separator. The design 
is optimized for low-energy radiative a-capture reactions 
important to stellar helium burning. It has a large angu- 
lar acceptance of ± 40 mrad, an energy accepta nce of ib 
7.5%, and a mass resolving power M/ AM ^ 100 (Couder 



et al. 2008). 



B. Underground facilities 

In all direct-kinematics capture-reaction measure- 
ments using 7 or neutron spectroscopy, whether per- 
formed above ground or underground, sources of envi- 
ronmental radioactivity must be controlled. Background 
sources include radioactivity from intrusions and impuri- 
ties in the rock and from construction materials, as well 
as sources intrinsic to targets and detectors. External 
sources can be reduced by careful shielding of the tar- 
get and the detector environment. In addition, beam- 
induced backgrounds (e.g., backgrounds from activation 
of impurities in the target) must be controlled through 
careful ion beam optics and choice of vacuum component 
materials. Active shielding techniques and complex event 
identification can also help. 

In surface facilities, however, the most difficult back- 
grounds are frequently those associated with cosmic rays. 
This background can be removed by exploiting the nat- 
ural shielding provided by the rock overburden in under- 
ground sites. The improvements possible with this strat- 
egy have been demonstrated by the 50 keV LUNA I and 
400 keV LUNA II programs at Gran Sasso. The labora- 
tory's depth, ~ 3.0 km.w.e. (kilometers of water equiva- 
lent, flat-site equivalent (Mei and Hime 2006)), reduces 



the fluxes of muons and secondary neutrons, relative to 
surface values, by factors of 10^ and lO'^, respectively. 



Consequently, the LUNA I collaboration (Bonetti et al. 



1999 ) was able to map S33 throughout the Gamow peak: 
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FIG. 18 Layout of the St. George recoil separator. 



a counting rate of one event per month was achieved at 
the lowest energy, E = 16 keV, with an uncertainty of 20 
fb or 2 X 10^'^^ cm^. Other critical pp chain and CNO 
cycle cross sections were made at energies far lower than 



previously possible (Formicola et al. 


2003 


Greife et al. 


1994 


Imbriani et al. 


2005 


Junker et al. 


19981. 



The successes of LUNA have inspired plans for the new 
underground facilities we discuss in this section. Figure 
[T9| shows a schematic of the present LUNA II set-up in 
Gran Sasso, which consists of a commercial 400 kV accel- 
erator, a windowless gas target, and a solid target line. 

Nuclear astrophysics has rather modest depth require- 
ments. The hadronic cosmic-ray component is quickly at- 
tenuated, leaving penetrating high-energy muons as the 
dominant source of background at depth. These muons 
interact in the rock to produce neutrons and a continuous 
spectrum of high energy 7s. Thus the main requirement 
is an overburden sufficient to reduce muon-associated ac- 
tivities to a level well below natural background levels 
associated with activities in the laboratory's rock and 
concrete wall s. The neutron fluxes in Gran Sasso, ~ 4 x 
10~^/cm^/s (Bemmerer et al. 2005 [Laubenstein et al. 



2004), and in Spain's underg round laboratory Canf ranc, 
JTm±Q.A4) X lO-VcmVs ( [Carmona et a/.t|2004[ ), are 
almost entirely due to local radioactivity. Taking these 
deep-laboratory values as typical of the environmental 
background component, one can determine the depth 
necessary to reduce the cosmic-ray-associated neutron 
contribution to 1% of the total. The simulations of |Mei| 
and Hime (2006) yield ~ 1.5 km.w.e. (flat site equiva- 



lent). 

Similar results are found for the 7-ray flux. The LUNA 
^^N(p,7) counting goal was 10~^ counts/keV/hr. The 
cosmic-ray muon-induced rate at 1.5 km.w.e. would 
be approximately an order of magnitude lower ( )Haxton| 



et al. 2007). As almost all deep physics laboratories 
now operating are at depths in excess of 1.5 km.w.e, one 
concludes that many locations are suitable for nuclear 
astrophysics - at least until order-of-magnitude reduc- 
tions in the laboratory environmental neutron and 7-ray 
background are made. 

Based on the success of the LUNA collaboration, sev- 
eral underground accelerator facilities are now being pro- 
posed. Table |XIII| shows the parameters of these facili- 
ties. The plans reflect design improvements from fifteen 
years of experience with LUNA. 

The present LUNA facility is small and limited to the 
measurement of proton- and a-capture reactions below 
400 keV, with typical beam currents between 100 and 
200 /iA. The available beam current has limited the sta- 
tistical accuracy of data taken at the lowest energies. In 
addition, many reactions have complex resonance struc- 
tures that must be adequately mapped, to provide the in- 
formation needed to extrapolate cross sections to Gamow 
energies. This requires measurements over a broader en- 
ergy range than is currently available at LUNA. There- 
fore, the LUNA collaboration has submitted a letter of 
intent for the installation of a higher energy accelerator 
that would allow the LUNA program to grow beyond so- 
lar fusion physics. This upgrade proposal is currently 
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FIG. 19 A schematic of the present LUNA 400 keV set-up. 



under review ( Prati et al. 2008 1 . 



Three initiatives for new underground accelerator fa- 
cilities are also under discussion: 

• ELENA is a proposed facility for the Boulby salt 
mine in the UK, a site that has environmental neu- 
tron backgrounds less than half those of Gran Sasso 
[(L72 ± 0.61 (st at.) ± 0.38 (syst )) x 1 0-^ /cm^/s 
above 0.5 MeV (Carmona et al.\ 2004)] and 7-ray 



backgrounds that are 5-30 times lower than Gran 
Sasso values, for < 3 MeV ('Aliotta:'2009'). This 
reflects the low U and Th concentrations ni salt. As 
the site is approximately at the same depth as Gran 
Sasso (2.8 vs. 3.1 km.w.e., taking proper account 



of the topography ( Mei and Hinie , 2006) ), full ad- 
vantage can be taken of the reduced environmental 
background. 

• CUNA is a 3 MeV accelerator facility that has 
been proposed for Spain's Canfranc Laboratory, lo- 
cated in an abandoned train tunnel in the Pyrenees 



• DIANA, Dakota Ion Accelerators for Nuclear As- 
trophysics, would be the nuclear astrophysics facil- 
ity for DUSEL (Deep Underground Science and En- 
gineering Laboratory), a laboratory being planned 
in the abandoned Homestake gold mine. South 



Dakota (DIANA Collaboration 2009). 



mountains (Bettini 2009) 



As in the case of the proposed LUNA upgrade, these 
facilities would be capable of mapping cross sections over 
broad energy ranges with fixed configurations for target 
and detector. 

We discuss DIANA is more detail, as an example of 
the improvements that would be possible in next gener- 
ation nuclear astrophysics facilities. DIANA's proposed 
site is the 4850-foot level of Homestake, the same level 
where Davis operated his chlorine detector. The design 
combines a low-energy 400 kV high-intensity accelerator, 
a high-energy accelerator with a maximum voltage of 3 
MV, and flexibly configured target stations and detector 
systems. Both accelerators will be coupled to a shared 
target station, in order to reduce uncertainties that would 
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FIG. 20 Proposed layout of the DIANA facility. 



arise when cross sections are measured at different facil- 
ities, with different targets and detector configurations. 
The accelerators will have a substantial overlap in their 
energy ranges due to the design of the ion source on the 
high- voltage platform of the low energy accelerator. This 
will reduce uncertainties in combining data sets. The pro- 



posed beam current of several mA is at least one order 
of magnitude higher than any presently available. This 
enhances counting rates, but also requires increased at- 
tention to beam-induced backgrounds as well as targets 
capable of handling the power. Figure [20| shows DIANA's 
conceptual design. 
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Appendix: Treating Uncertainties 

A. Introduction 

This section describes our method for dealing with dis- 
crepant data sets that may occur, for example, when de- 
riving recommended S(0) values from experimental mea- 
surements of nuclear reaction cross sections. 

While the conventional minimization method is ade- 
quate for analysing data sets that are in good agreement, 
there is no rigorous method for dealing with discrepant 
data sets and their underlying unidentified systematics. 
But reasonable procedures exist. In Solar Fusion II we 
adopt the Scale Factor method, here called the i nflation 
factor method (IFM), that is used by the PDG (Amsler 
et al. 2008). In this method, the fit errors from a con 



ventional x minimization are inflated by a factor that 
depends on \/x^Ji', where v is the number of degrees of 
freedom. This method is well known, widely used, and 
straightforward to apply. 
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TABLE XIII Attributes of proposed second-generation underground facilities for nuclear astrophysics. 





LUNA 


DIANA 


ELENA 


CUNA 


FACILITY 


Tiphoratorv TTTiHerPTOiinr] 


Dakota Ton Accelerators 


Exp6riniGiit3il Low-EiiGrgy 


Canfranc Nuclear 




for Nuclear Astrophysics 


for Nuclear Astrophysics 


Nuclear Astrophysics 


Astrophysics Facility 


Location 


Gran Sasso, Italy 


Homcstakc Mine, USA 


Boulby Mine, UK 


Canfranc, Spain 


Rock type 


hard limestone 


metamorphic rock 


salt 


hard limestone 


Depth 


3.1 


4.3 


2.8 


~ 2.0 


(km.w.e, flat site) 










Low energy 


50-400 keV 


50-400 kV HV olatform 


none 


none 


accelerator 


0.5-1.0 mA 

RF ion source (pjlle"*") 


> 10mA 

ECR ion sources, 
single, multiply charged 






High energy 


0.4-3.5 MeV 


0.35-3.0 MeV 


3.0 MeV accelerator 


up to 5.0 MeV 


accelerator 


electrostatic 


electrostatic 


electrostatic 


electrostatic 




up to 0.3 mA 


up to 10 mA 


0.5 mA 






ECR ion source 


ECR ion sources 
single, multiply charged 


ECR ion source 





While the IFM is the only one discussed in the PDG 
Introduction, alternatives exist. We discuss some exam- 
ples at the end of this Appendix. 



B. The inflation factor method 

The IFM addresses systematic uncertainties when 
combining results from different and possibly discrepant 
data sets. The method inflates errors in proportion to 
the quoted errors originally given by the experimenters. 

Discrepant data may be defined by the P- value of the 
fit, where P = P{xl^,v) is the probability of obtaining 
a value at least as large as the observed value. The 
inflation factor is conventionally chosen to be \fx^Jv and 
is commonly applied in cases where 'X^ jv > 1. We use an 
alternative inflation factor \/x^lx^{P = 0-5) to account 
for the fact that, for small v and non-discrepant data, 
the expected value of is smaller than unity. For large 
!/, the two scaling factors are equivalent. 

The IFM scales all experimental errors by the same 
fractional amount, resiilting in equal internal and exter- 
nal errors on the mean. Because one generally cannot 
identify a specific mechanism accounting for discrepant 
data, this procedure (like all other procedures) has no 
rigorous mathematical justification. However qualitative 
arguments support its reasonableness. As the method 
maintains the relative precision of discrepant data sets, 
it apportions a larger absolute fraction of the identified 
systematic error to the less precise data sets. This is 
consistent with naive expectations that a large, uniden- 
tified systematic error is more likely to "hide" within a 
low-precision data set than within a high-precision one, 
given the advantages a high-precision data set offers an 



experimentalist who does "due-diligence" cross checks to 
identify systematic errors. The IFM is generally con- 
sidered the most appropriate procedure in the absence 
of information that would support alternatives, such as 
omitting certain data, or increasing errors on some data 
but not others. 

We employ error inflation whenever x^ > X^{P — 0-5), 
and no error scaling otherwise. With this general rule, 
errors are inflated a bit even when x^ is only slightly 
in excess of x^{P = 0.5), despite the lack of compelling 
evidence of discrepancy in such a case. This procedure 
yields a continuous formula and avoids the introduction 
of an arbitrary threshold for inflation. 

In extreme cases one may obtain errors that are 
deemed too small. For example, when analyzing data 
containing a few results with small errors and a larger 
number of results with large errors, the large-error data 
will reduce the error on the mean by increasing i^, even 
though they may have little effect on the central value. 
In such a case, we agree with the PDG's recommendation 
that, to mitigate this problem, data be excluded which 
have an error larger than some (arbitrary) limit, specif- 
ically 3SVN, where N is the number of measurements 
and 6 is the unsealed error on the mean. However, ap- 
plying this exclusion criterion may not be adequate to 
resolve this difficulty in all cases. 

While the IFM makes no assumptions about the rea- 
sons for discrepant data, in actual applications it may be 
apparent that not all data sets are equally reliable. In 
such cases judgment is necessary, and data selection is 
appropriate. Data should be discarded if the error anal- 
ysis is poorly documented or inadequate. Data may be 
discarded if the procedure used to generate them involves 
questionable assumptions, or if corrections were not made 
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for effects now known to be important. Data errors may 
be modified (e.g. increased) if such new information is 
available. 



C. Application of the inflation factor method 

The following is based on the discussion in the Intro- 
duction of the PDG compilation of Amsler et al. (20081: 



1. In general, statistical and systematic data errors 
should be identified and specified separately. Sys- 
tematic errors should be subdivided into varying 
(random) and common-mode (normalization) er- 
rors. For a single data set, normally the statistical 
and varying systematic errors should be combined 
in quadrature and used as data errors in a min- 
imization to determine unknown parameters. The 
resulting fit error (s) should be multiplied by the 
inflation factor (see below). The common- mode er- 
ror is then folded in quadrature with the inflated fit 
error to determine the overall normalization error. 

For multiple data sets, the systematic errors should 
be examined to determine if they are independent 
among the different data sets. Parameters deter- 
mined from multiple, independent data sets may be 
combined in a separate minimization in which 
each parameter value is characterized by its to- 
tal error determined by combining statistical and 
systematic (normalization) errors in quadrature. 
Again, this fit error should be multiplied by the 
inflation factor. If the systematic errors in differ- 
ent data sets are correlated, then this correlation 
must be taken into account in the fitting. A conve- 
nient method for handling correlations is described 
in the 2008 PDG compilation. 

2. Whenever > X^{P = 0-5) the fit errors should 
be increased by the multiplicative inflation factor 
\/x^/x^{P = 0-5), where x^{P = 0.5) is the x^ 
corresponding to a P value of 0.5 for v degrees of 
freedom. The x^ and v should be stated, along with 
the inflation factor when it is larger than unity. 
Large reported inflation factors serve to alert the 
reader to potential problems. 



3. Data with uncertainties larger than 3v7V(5, where 
N is the number of measurements and 5 is the (un- 
sealed) error on the mean should be excluded. One 
should be aware of possible error underestimation 
in certain cases as mentioned above. The resolution 
of such situations may require additional judgment. 



D. Other methods 

Other error analysis methods follow somewhat differ- 
ent strategies. The cost function methods used in CO- 
DATA analyses (Cohen and Taylor 19871 are designed 



to reduce the x^ by selective re-weighting of data; i.e. by 
increasing the errors nonuniformly on the data, in such 
a manner as to minimize the "cost", i.e. the error on the 
mean. Alternatively, D'Agostini (1994) has advocated 



a procedure for fitting multiple data sets in which one 
minimizes the sum of a data x^ and a normalization x^- 
One method that has been applied to the analysis of so- 
lar fusion cross section is t hat of Cyburt ( 2004, ) (see also 
Cyburt and Davids (2008)). This approach introduces a 
"discrepancy error", a disc, that is added in quadrature 
with the normalization errors of individual experiments 
when fitting mixed data sets. Effectively this procedure 
distributes the unexplained discrepancy equally over the 
data sets, regardless of their stated accuracy, in contrast 
to the PDG procedure, which assigns the discrepancy in 
way that preserves the relative stated accuracy of data 
sets. The Cyburt (2004) method leads, in cases where 



there is excess dispersion, to increased de-weighting of 
the more precise data points, compared to the IFM. In 
addition, the contribution of (Xdisc to the error of the mean 
does not decrease as the number of measurements N in- 



creases. 



The Cyburt ( 2004 ) and IFM methods refiect two limits 



in how one apportions an unexplained discrepancy among 
data sets: one could construct other models that inter- 
polate between these two limits (equal vs. proportionate 
allocation of the discrepancy error). The argument for 
the IFM procedure has been stated previously: it is eas- 
ier to miss a large systematic error within a low-quality 
data set than within a high-quality one. In addition, it 
avoids a situation where archival data of poor quality, 
containing an unidentified systematic error, unduly im- 
pact the weight that would otherwise be accorded a new 
experiment of exceptional quality - thereby inappropri- 
ately diluting the impact of the best results. Alternatives 
to the IFM tend to produce roughly equivalent results un- 
less the discrepancies among data sets are large. We are 
fortunate in this paper to be dealing with discrepancies 
that are modest. 
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